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Summary  of  Results 
(keyed  according  to  list  of  publications) 


1.  This  paper  represents  the  first  operation  of  a  molecular  beam  maser  in 
the  submiilirneter  spectral  region.  The  extremely  high  resolution  of  this 
device  makes  it  an  especially  powerful  spectroscopic  tool  and  presents  the 
possibility  of  the  construction  of  a  primary  frequency  standard  in  the 
submillimeter. 

2.  This  paper  reports  the  first  pressure  broadening  study  on  HgS.  Because 
of  a  paucity  (none)  of  transitions  at  conventional  microwave  frequencies,  this 
important  species  had  not  been  previously  studied.  More  importantly  our 
theoretical  calculations  (Anderson  theory)  show  that  the  theory  w as  in  fact 
generally  valid  (  to  10%)  and  that  previous  inconsistencies  (^  100”)  were 
due  to  inaccurate  infrared  data  rather  than  theoretical  shortcomings. 

3.  This  paper  reports  the  first  observation  of  symmetric  top  forbidden 
transitions  in  the  NMMW  spectral  region.  The  results  lead  to  a  substantially 
improved  knowledge  of  the  energy  level  structure  of  PII„ 

O 

4.  This  paper  demonstrates  the  power  of  millimeter  and  submillimeter 
spectroscopy  as  a  diagnostic  probe  of  molecular  lasers.  Several  conventional 


wisdoms  about  the  HCN  laser  were  shown  to  be  untrue.  This  project  and  its 
spin  offs  will  become  a  major  part  of  our  effort  in  the  next  contract  period. 

5.  This  is  a  short  note  that  relates  centrifugal  distortions  to  changes  in 
molecular  constants  and  is  a  spin  off  of  paper  1. 

o 

6.  This  paper  reports  the  results  of  a  high  temperature  1500  C)  study  of 
SiO  and  represents  the  first  laboratory  observation  of  one  of  the  interstellar 
masers. 

7.  This  is  an  invited  review  of  the  applications  and  techniques  of  millimeter 
and  submillimeter  spectroscopy. 

8.  Molecular  laser  plasmas  are  rich  sources  of  excited  vibrational  states. 

This  paper  reports  the  use  of  the  HCN  laser  plasma  as  a  source  for  these 
excited  state  studies.  An  interesting  symbiosis  between  this  project  and  that 
reported  in  4)  exists.  So  little  is  known  about  the  excited  stales  of  molecules 

in  the  millimeter  and  submillimeter  that  one  ordinarily  must  do  the  spectroscopic 
project  first  so  that  you  know  what  you  are  diagnosing. 

9.  This  paper  represents  an  extension  of  the  work  reported  in  3  to  Asll^  • 


10.  Because  of  its  extreme  centrifugal  distortion,  water  is  the  most  difficult 
'2  molecule  to  model.  This  paper  reports  a  combined  analysis  of  our 
millimeter  and  submillimeter  microwave  data,  new  FIR  Fourier  transform 
data,  and  energy  levels  derived  from  IR  combination  differences  to  upgrade 
the  energy  levels  of  AFGL  spectral  table. 

11.  This  paper  reports  the  results  of  an  extensive  project  to  model  the  HNC 

o 

spectrum  from  0-300  GHz.  The  primary  result  is  a  model  that  predicts 
~  5000  transitions  in  this  region  -a  1  MHz.  I  understand  that  this  work  is  of 
some  interest  to  BMD  for  their  modeling  of  the  "disturbed"  atmosphere. 

12.  This  work  has  been  the  subject  of  two  papers  deFivored  at  the  Columbus 
Meeting  and  is  being  finalized  for  inclusion  in  a  thesis  and  journal  publication. 
Again,  I  understand  that  this  work  (especially  H09)  relates  to  BMD  interests. 

Cl 

13.  This  is  a  spectroscopic  project  on  the  millimeter  and  submillimeter 
spectrum  of  the  extremely  reactive  and  unstable  mixed  diatomic  halides. 

14.  This  is  a  paper  that  reports  an  interesting  and  potentially  very  valuable 
approach  to  the  theoretical  modeling  of  the  absorptions  of  hot  water  vapor. 

15.  and  16.  These  papers  represent  the  last  papers  from  Professor  Gordy's 


With  the  exception  of  12.  (which  has  shared  support  with  NASA)  the 
work  reported  here  was  supported  by  no  other  federal  agency  than  ARQ. 
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Extension  of  high  resolution  beam  maser  spectroscopy  into  the 
submillimetre  wave  region1 

R.  Michael  Garvey2  and  Frank  C.  De  Lucia 

Department  t>1  Physics.  Duke  I  hiiccrsi/v.  Durham.  \(  2770b,  l  !.S.A. 

Received  December  20.  1976 

High  resolution  molecular  beam  maser  spectroscopy  has  been  extended  into  the  submilli¬ 
metre  wave  region  of  the  spectrum.  The  maser  spectrometer  consists  of  a  multiple  beam  source, 
stacked  quadrupole  state  selector,  and  Fabrv  Perot  cavity.  The  submillimetre  ware  stimulating 
signal  is  produced  by  harmonic  generation  from  a  millimetre-wave  klystron,  directly  referenced 
to  WWVB.  The  maser  signal  is  detected  in  a  cryogenic  photodetector.  Linew  idths  50  100  times 
narrower  than  the  Doppler  width  and  good  signal-to-noise  ratios  were  obtained.  This  system 
has  been  used  to  resolve  the  small  deuterium  hyperfine  splittings  of  both  the  1 ,0  -*  l0i  transi¬ 
tion  of  D;0  at  317  GH/  and  the  •/=!-■»  0  transition  of  ND,  at  309  GHz.  The  relevance  of 
the  D.O  quadrupole  interaction  data  to  the  electric  field  gradient  in  water  is  discussed.  The 
spectroscopic  constants  in  kilohertz  for  D.O  in  the  states  I  and  l„,  are,  respectively,  ( cQqDo 
=  34.81  +  0.32,  Co  =  —2.661,  and  U’Q‘Ij)d  =  -29.60  ±  0.55,  C„  =  -2.532.  The  constants 
for  NO,  are,  in  kilohertz:  (e@i/j)N  =  816.2  +  1.1,  Cs  =  2.0  +  1.5,  U’0‘h)o  =  19.9  ±  0.8, 
and  rD  =  -  2.2  ±  0.8. 


La  spectroscopie  a  haute  resolution  par  maser  et  faisccau  nioleculaire  a  etc  etendue  a  la 
region  vies  longueurs  d'onde  inframillimetriques.  Le  speetrometre  maser  est  constitue  par  unc 
source  it  multiples  faisceaux,  un  selecteur  d'etats  quadrupolaire  et  une  cavite  Fabry  Perot. 
Le  signal  stimulant  les  ondes  inframillimetriques  est  produit,  par  generation  d’harmoniques, 
it  partir  d’un  klystron  fonctionnant  dans  la  region  millimetrique,  avec  reference  dirccte  aux 
signaux  WWVB.  Le  signal  maser  est  detecte  dans  un  photodetecteur  cryogenique.  On  a  obtenu 
des  raies  ayant  vies  largeurs  de  50  a  100  fois  infericures  a  la  largeur  Doppler  et  de  bonnes 
valeurs  du  rapport  signal-bruit.  Le  systeme  a  ete  utilise  pour  resoudre  les  faibles  separations  de 
structure  hyperfine  du  deuterium  dans  les  deux  transitions  l,0->  I o i  de  D.O.  a  317GHz,  et 
J  =  l  ->  Ovle  NO  ,,  iv  309  GHz.  La  povtee  vies  donnees  concernant  1‘ interaction  quadrupolaire  de 
D.O  sur  la  question  du  gradient  du  champ  electriquc  dans  I'eau  est  discutee.  Les  valeurs  cn 
kilohertz  des  constantcs  sped roscopiq ues  de  D.O.  pour  les  etats  I  et  I  sont  respeetivement 
(cCh/jlo  =  34.81  ±  0.32,  C„  =  -2.661  et  U’Q</j)u  =  -29.60  ±  0.55.  C„  =  —  2.532.  Les 
constantcs  pour  ND,  sont.  en  kilohertz:  ieQt/j) s  =  816.2  +  1.1.  CN  =  2.0  +  1.5.  (cQ‘ljhs  = 
19.9  +  0.8  et  CD  =  -2.2  +  0.8. 

ITraduit  par  le  journal] 

Can.  J.  Plus..  55.  111511977)  1  1  1 


I.  Introduction 

Molecular  beam  maser  spectroscopy  was  first 
reported  by  Gordon  (1955)  in  the  centimetre 
microwave  region  and  the  high  resolution 
capabilities  of  this  technique  have  now  been 
used  to  study  numerous  molecular  species  in  this 
spectral  region.  A  number  of  proposals  have 
been  made  to  extend  beam  maser  spectroscopy 
into  the  submillimetre  wavelength  region  and 
beam  absorption  spectroscopy  has  been  reported 
in  this  region  by  van  Dijk  and  Dymanus  (1970). 
These  proposals  have  been  motivated  by  the 
general  desire  to  extend  the  capabilities  of  these 
devices  as  well  as  by  the  need  for  studies  of 

'This  work  supported  by  the  U.S.  Army  Research 
Office,  Grant  DAAG29-77-G-0007. 

1  Present  address:  National  Bureau  of  Standards,  325 
Broadway,  Boulder,  CO  80302,  U.S.A. 


specific  transitions.  For  many  of  the  fundamental 
and  interesting  molecules,  these  submillimetre 
transitions  are  either  the  only  ones  available  or 
are  more  advantageous  to  study  than  the  higher 
./  transitions  which  fall  at  longer  wavelengths. 
Low  J  transitions  exhibit  the  largest  hyperline 
splittings  and  are  least  affected  by  centrifugal 
distortion.  In  this  paper  we  report  the  extension 
of  high  resolution  molecular  beam  maser  spec¬ 
troscopy  into  the  submillimetre  region  and  the 
study  of  the  J  =  I  ->  0  transition  of  ND,  at 
309  GHz  and  the  l,0  -»  ln,  transition  of  6,0 
at  317  GHz. 

II.  F.vpcrinicntal 

Numerous  articles  which  describe  centimetre 
wave  molecular  beam  maser  spectrometers  have 
appeared  in  the  literature  (Gordon  vt  «/.  1955: 
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Fig.  I .  Submillimctre  beam  maser  spectrometer. 


Thaddeus  et  al.  1964).  In  addition,  molecular 
beam  masers  for  88  GHz  (Marcuse  1961).  72 
GHz  (Krupnov  el  al.  1964).  and  various  fre¬ 
quencies  between  50  and  200  GHz  (De  Lucia 
and  Gordy  1969:  De  Lucia  and  Cederberg  1971 ; 
Garvey  and  De  Lucia  1974:  Garvey  el  al.  1976) 
have  been  described.  Because  of  the  difference  in 
wavelength  between  the  more  common  centi¬ 
metre  wave  devices  and  the  one  described  here, 
a  simple  scaling  down  of  dimensions  is  not  a 
satisfactory  design  approach.  The  operation  of  a 
molecular  beam  maser  may  be  divided  as  follows: 
( /)  beam  formation,  (2)  state  selection.  (2)  micro- 
wave  cavity.  (4)  generation  and  detection  of  sub¬ 
millimetre-wave  energy,  (5)  frequency  measure¬ 
ment.  Figure  I  shows  a  diagram  of  the  experi¬ 
mental  apparatus.  Multiple  molecular  beams  are 
formed  in  the  beam  source,  pass  through  a 
stacked  quadrupole  electrostatic  state  selector, 
and  interact  with  the  microwave  field  of  a  tunable 
Fabry  Perot  cavity.  The  microwave  field  in  this 
cavity  is  produced  by  harmonic  generation  from 
a  millimetre  wave  klystron.  The  maser  output 
signal  is  detected  by  a  helium  temperature  (1.6  K) 
InSb  photoconducting  detector. 

A.  Beam  Formation 

Molecular  beam  masers  have  ordinarily  used 
effusive  beam  sources.  Kantrowitz  and  Grey 


(1951)  proposed  the  use  of  supersonic  jets  to 
obtain  beams  with  greater  forward  directivity 
and  intensity.  Until  recently  these  sources  have 
not  been  used  in  beam  masers.  These  supersonic 
nozzles  produce  increased  forward  beam  inten¬ 
sities  compared  to  effusive  sources,  but  for  beam 
spectroscopy  the  greatest  gains  result  from  the 
cooling  effects  of  the  supersonic  expansion 
(Dyke  el  al.  1972).  Rotational  cooling  results  in 
partition  functions  which  are  substantially 
smaller  than  can  be  obtained  in  a  bulk  gas  sample 
before  condensation.  The  increase  in  population 
which  results  is  largest  for  low  J  transitions.  This 
is  fortuitous  because  of  the  previously  discussed 
spectroscopic  advantages  of  low  J  transitions. 

Because  of  the  large  frontal  area  of  Fabry- 
Perot  cavities,  it  is  advantageous  to  employ 
multiple  beams  to  increase  the  filling  factor  of 
the  cavity.  We  have  used  a  seven  by  nine  array 
of  nozzles,  each  beam  having  its  ow  n  electrostatic 
quadrupole  state  selector.  Collimation  of  each 
beam  is  accomplished  by  an  arrangement  of 
cryogenically  cooled  channels.  It  has  been  found 
that  the  complex  and  precise  nozzle  geometries 
originally  used  in  studies  of  supersonic  jet 
phenomena  may  be  replaced  by  small  holes  in  a 
thin  metal  sheet.  To  take  full  advantage  of  the 
nozzle  source,  it  is  desirable  to  expand  the 
molecules  into  the  maser  vacuum  chamber  from 
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as  high  a  source  pressure  as  is  feasible.  Ulti¬ 
mately.  the  gas  load  presented  to  the  collimation 
structure  limits  the  maximum  source  pressure 
that  may  be  used.  With  holes  of  0.34  mm  diam¬ 
eter.  approximately  10  Torr  source  pressure  "as 
found  to  be  appropriate  in  our  experiment.  Al¬ 
though  this  design  is  an  order  of  magnitude  from 
the  source  pressures  used  in  single  beam  systems, 
we  have  realized  gains  of  approximately  five  over 
a  similar  effusive  source. 

tt.  Stale  Selection 

Population  inversion  may  produce  intensity 
gains  of  the  order  of  k  T  hv  for  transitions  of 
frequency  v.  Although  this  factor  is  inversely 
proportional  to  the  frequency,  gains  of  20  may 
be  realized  at  300  GHz.  The  process  of  state 
selection  with  inhomogeneous  electric  fields  has 
been  the  subject  of  extensive  discussion  (Beck'r 
and  Ange\\  l%3;  Mednikov  and  Parygin  1963). 
A  number  of  different  state  selectors  have  been 
successfully  employed,  but  the  multiple  beam 
arrangement  allowed  by  the  submillimetre  wave 
Fabry  Perot  cavity  places  constraints  upon  the 
state  selector  geometries  which  can  be  used 
efficiently.  Figure  2  shows  an  end  view  (beam 
flow  perpendicular  to  figure)  of  the  state  selector 
used  for  this  experiment  (De  Lucia  and  Gordy 
1 969).  The  poles  of  the  quadrupolc  are  fabricated 
from  copper  rods  of  3  16-in.  diameter,  spaced 
3  8-in.  on  centers.  The  structure  is  1 2  in.  long  and 
is  cryogcnically  pumped  from  within  by  cooled 
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copper  tins.  Arcing  within  the  state  selector 
places  an  upper  limit  on  the  applied  voltage  of 
about  35  kV. 

C.  Cavity 

At  millimetre  and  submillinictre  wavelengths 
the  closed  cylindrical  cavities  used  as  resonant 
structures  for  centimetre  wave  masers  become 
significantly  less  advantageous  than  quasi- 
optical  Fabry  Perot  cavities.  Because  the  mole¬ 
cular  time-of-llight  in  the  microwave  field  is  the 
ultimate  limit  to  maser  resolution,  the  diameter 
of  the  microwave  mode  in  the  cavity  should  be 
maximized.  For  l  abry  Perot  cavities  made  with 
the  usual  spherical  mirrors,  large  mode  diam¬ 
eters  are  difficult  to  obtain  at  submillimetre  wave¬ 
lengths  and  the  modes  themselves  arc  sensitive  to 
imperfections  in  the  long  focal  length  mirrors 
which  are  required. 

We  have  developed  a  Fabry  Perot  mirror, 
shown  in  Fig.  3.  which  overcomes  the  submilli¬ 
metre-wave  mode  size  problem.  The  center  of  the 
mirror  is  flat:  the  edges  are  curved  to  confine  the 
microwave  mode  in  a  manner  analogous  to  that 
of  conventional  spherical  mirrors.  This  mirror 
may  be  cut  in  a  stepwise  fashion  and  polished  to 
the  final  curve  as  illustrated  in  Fig.  3.  The  cavity 
is  formed  from  this  curved  mirror  and  an  op¬ 
posing  plane  mirror,  both  fabricated  from  cop¬ 
per.  Diffraction  losses  are  insignificant.  Q' s  ap¬ 
proaching  the  limit  imposed  by  copper  losses  are 
attainable,  and  high  order  modes  are  effectively 
suppressed.  This  new  mirror  design  has  the 
added  advantage  that  the  mode  diameter  is  only 
a  very  slow  function  of  wavelength  (essentially 
filling  thecav  itv  to  the  curv  ed  edges)  and  is  there¬ 
by  usable  over  a  wide  wavelength  range.  Power 
is  coupled  into  and  out  of  the  cavity  through 
small  holes  in  the  plane  mirror.  Tuning  is  ac¬ 
complished  by  translation  of  one  of  the  mirrors. 

For  the  cavity  used  in  this  work  (the  mirror  of 
Fig.  3  and  a  plane  mirror  separated  by  about 
20  cm)  a  mode  of  approximately  10  cm  diameter 
results.  The  Q  is  approximately  proportional  to 
v1  *’  and  at  I  mm.  Q  %  500000. 


I  Hi.  I  ahry  l*erol  mirror  fabrication. 
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/).  Generation  and  Detection  of  Subntilliinetre 
Radiation 

Harmonic  generation  (King  and  Gordy  1953) 
is  a  proven  technique  for  production  of  modest 
amounts  of  submillimetre  radiation  and  is  par¬ 
ticularly  attractive  for  beam  maser  spectroscopy; 
relaxation  rates  in  a  molecular  beam  are  slow 
and  saturation  occurs  at  relatively  low  power 
levels.  With  the  high  Q  cavity  used,  we  have 
found  it  possible  to  generate  enough  power  to 
saturate  the  ll0  -»  l()l  transition  of  D.O  at  317 
GHz  from  a  300  mW  53  GHz  klystron  (OKI- 
55V 11). 

Since  the  observed  linewidth  depends  upon 
spectral  purity  of  the  stimulating  source,  it  is 
important  that  this  source  be  as  stable  as  pos¬ 
sible.  Free  running  klystron  oscillators  arc  too 
unstable  and  it  is  necessary  to  lock  the  signal 
source  to  a  stable  crystal  oscillator.  A  diagram  of 
this  stabilization  system  is  included  in  Fig.  I. 
The  Frequency  Engineering  Labs  1 33A  Synchro¬ 
nizer  contains  a  variable  crystal  oscillator  and 
multiplier  amplifier  chain  which  produces  a  48 
MHz  frequency  comb  extending  to  several  giga¬ 
hertz.  A  variable  sideband  is  added  and  the  comb 
with  sidebands  is  amplified  in  a  traveling  wave 
tube  and  applied  to  a  crystal  mixer  (1N26) 
coupled  to  the  klystron.  One  of  the  resulting  beat 
frequencies  is  down-converted  to  27  MHz  in  a 
modified  Clarke  Instruments  receiver.  Tnis  27 
MHz  signal  is  returned  to  the  synchronizer  where 
it  is  phase  compared  to  a  crystal  reference  oscil¬ 
lator.  An  error  signal  is  generated  and  used  to 
stabilize  the  klystron.  The  klystron  is  swept  in 
frequency  by  means  of  the  variable  sideband 
oscillator. 

The  very  low  power  levels  at  which  saturation 
of  maser  levels  occur  make  superheterodyne 
linear  detectors,  which  arc  almost  universally 
used  in  centimetre  wave  masers,  more  attractive 
than  square-law  detectors.  However,  at  sub¬ 
millimetre  wavelengths,  local  oscillators  and 
mixers  for  superheterodyne  detection  arc  either 
unavailable  or  prohibitively  expensive.  We  have 
found  that  the  square-law  InSb  photodetector 
cooled  to  1.6  K  has  adequate  sensitivity  for  these 
experiments  in  the  submillimetre  region  (Hel- 
minger  et  at.  1971).  This  results  partly  from  the 
increased  power  output  and  saturation  threshold 
of  submillimctrc  masers  compared  to  centimetre 


masers  and  partly  from  the  increased  sensitivity 
of  the  cryogenic  detector  compared  to  the  room 
temperature  detectors  more  commonly  employed. 
This  detector  is  broadbanded,  requires  no  local 
oscillator,  and  is  extremely  reliable;  it  has  proven 
to  be  very  satisfactory  for  our  experiments.  The 
output  signals  arc  recovered  by  means  of  phase 
sensitive  detection  (at  3  kHz)  of  the  signal  pro¬ 
duced  as  a  result  of  the  demodulation  by  the 
maser  line  of  a  small  1.5  kHz  frequency  modu¬ 
lation  of  the  klystron. 

E.  Frequency  Measurement 

A  frequency  measurement  system,  separate 
from  the  frequency  stabilization  system,  is  used 
to  avoid  errors  that  might  be  introduced  either 
by  inexact  loop  locks  or  by  instabilities  or  in¬ 
accuracies  in  the  several  oscillators.  For  the  fre¬ 
quency  measurement  system,  a  5  MHz  crystal 
oscillator  (General  Radio  III5-B),  phase  com¬ 
pared  to  WWVB  (HPI 1 7A).  generates  a  standard 
reference  signal  which  is  amplified  and  multiplied 
to  provide  a  60  MHz  comb  extending  to  several 
GHz.  This  comb  is  mixed  with  a  sampled  portion 
of  the  stabilized  klystron  signal  and  the  beat  fre¬ 
quency  vb  is  amplified  and  mixed  with  a  counted 
reference  oscillator  (Gertch  FM-3)  in  a  com¬ 
munications  receiver.  As  vb  is  swept  (as  the 
stabilized  klystron  is  swept)  the  resulting  zero 
beat  is  recorded  simultaneously  with  the  spectra 
on  a  dual  channel  chart  recorder.  Ten  kilohertz 
sidebands  may  be  added  to  the  reference  oscilla¬ 
tor  signal  to  provide  a  lattice  of  zero  beats  w  htch 
serves  as  a  frequency  dispersion  calibration  on 
the  chart  paper. 

111.  Theory 

The  theory  of  molecular  hyperfine  structure 
has  been  extensively  discussed  in  the  literature 
(Thaddeus  et  at.  1964:  Garvey  ct  at.  1976:  Cook 
and  De  Lucia  1971;  Cederberg  1972).  Powerful 
techniques  have  been  developed  which  make  the 
calculation  of  transition  frequencies  and  intensi¬ 
ties  straightforward.  We  use  the  results  and  nota¬ 
tion  of  Garvey  et  at.  (1976)  and  Cook  and  De 
Lucia  (1971). 

The  two  identical  deuterium  nuclei  of  D,0 
make  the  only  contributions  to  the  nuclear 
hyperfine  structure.  The  Hamiltonian  for  the 
relevant  interactions  is: 


of  =  £  '  F(0,)<?(0,)  +  /(D,)  C(/),)  i  +  /(/),) f>(/>,0,)/(02) 

i  I  ® 


[i] 
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where  the  terms  included  represent  the  electric  quadrupole,  n  gnetic  dipole,  and  spin-spin  inter¬ 
actions,  respectively.  This  problem  is  identical  in  form  to  the  D2S  investigation  reported  by  De  Lucia 
and  Cederberg  (1971)  and  the  matrix  elements  for  the  D20  hyperfine  interactions  may  be  found  in 
that  work. 

Analysis  of  the  hyperfine  structure  of  ND2  is  considerably  more  complex,  partially  because  of  the 
increased  number  of  interactions  but  more  significantly  because  of  the  three  identical  particle 
problem.  Application  of  the  techniques  of  Cederberg  (1972)  circumvents  ambiguities  which  may 
occur  in  such  an  analysis.  The  appropriate  Hamiltonian  for  the  spectrum  observed  is: 

[2]  ye  =  ^K(N)g(N)  +  I(N)-C(N)-J  +  X  V{D)-Q(Di)  +  HD,yCiD,)-J 

+  HN)-D(ND,)-i(D,)  1  +  '  X  [/(D,)-D(D,D,)-/<D/)] 
J  1  <-i=  i 


where  the  terms  included  represent  the  nitrogen 
electric  quadrupole  and  magnetic  dipole  inter¬ 
actions,  the  deuterium  electric  quadrupole  and 
magnetic  dipole  interactions,  the  nitrogen-deu¬ 
terium  spin-spin  interaction,  and  the  deuterium- 
deuterium  spin-spin  interaction,  respectively. 
The  matrix  elements  of  these  terms  are  included 
in  the  Appendix  of  this  paper.  The  expressions 
which  result  lend  themselves  to  machine  compu¬ 
tation  and  the  energy  levels  and  transition  fre¬ 
quencies  are  obtained  by  numerical  diagonaliza- 
tion  of  the  resulting  Hamiltonian  matrix. 

IV.  Results 

The  experimental  beam  maser  spectrum  of  the 
1I0  -+  l0,  transition  of  D,0  is  shown  in  Fig.  4 
along  with  a  theoretical  spectrum  calculated  from 
the  hyperfine  constants  of  Table  I.  Table  2  show  s 
both  the  observed  and  calculated  transition  fre¬ 
quencies.  The  spin-spin  interaction  constant 
may  be  calculated  from  the  molecular  geometry. 
Because  of  the  small  size  of  the  nuclear  g-factor 
for  deuterium  compared  to  that  of  hydrogen,  the 
spin  rotation  constants  C,l0  and  C,01  can  be 
calculated  from  previous  beam  maser  work  on 
H20  and  HDO  (Bluyssen  el  al.  \961a,b)  with 
higher  accuracy  than  they  can  be  determined 
from  the  D20  spectrum.  The  deuterium  quadru¬ 
pole  constants  of  Table  I  were  calculated  from 
the  observed  spectrum.  Comparison  of  the  experi¬ 
mental  and  theoretical  data  shows  agreement 
well  within  the  experimental  error. 

The  large  nitrogen  quadrupole  interaction  in 
NDj  produces  a  triplet  which  may  be  resolved  in 
a  conventional  absorption  spectrometer,  but  the 
much  smaller  splittings  which  result  from  the 
deuterium  quadrupole  interaction  are  obscured 
by  the  Doppler  broadening.  Figures  5-7  show 


the  resolution  of  this  structure  obtained  in  this 
experiment.  As  in  the  case  of  D.O  .the  spin-spin 
interactions  may  be  calculated  from  the  molecu¬ 
lar  geometry.  The  remaining  hyperfine  constants 
have  been  evaluated  from  the  beam  maser  spec¬ 
trum  and  appear  in  Table  3.  Table  4  show  s  a  com- 


Fig.  4.  Experimental  and  theoretical  spectra  of  the 
I  io  -*  I o i  transition  in  D20. 
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Tabu  1.  Hyperfine  constants  of 
0.0  for  l,D  -*  1„, 


I 


10 


1 


01 


Constant 


U‘Q‘Ij)  o 

CD 


<11  £>....  Il> 


U‘Q(Ij)o 

Co 

<11  DJ:  II) 


Value  (kHz) 


34.81  ±0.32“ 
-2.661'’ 
-0.163“ 

—  29 . 60  ±  0 . 5S 
-2.532'’ 
0.325“ 


•Uncertainties  represent  one  standard  deviation  of  the 
experimental  data. 

^Calculated  from  Bluyssen  ct  al.  (I967«). 

'Calculated  from  molecular  geometry. 


Tabu  2.  Hyperfine  spectrum  of  13.0  for  1  I0, 


Measured" 

Calculated 

Difference 

F 

(kHz) 

(kHz) 

(kHz) 

0  ->  1 

1  ->  0 

-46. 6±  1.8 

-46.9 

0.3 

2  — *  2 

-  16.3+  1 .4 

-  16.5 

0.2 

2  -  1 

23 . 3  ±  1 . 2 

24.4 

-1.1 

1  ->  2 

41.5+1.2 

41.5 

0.0 

1  -  1 

82. 4±  1 .2 

82.1 

0.3 

“Measured  relative  to  316  799  845  kHz.  Uncertainties  represent  one 
standard  deviation  of  the  experimental  data. 


Fie;.  5.  Experimental  and  theoretical  spectra  of  the 
J  =  I  -»  0 ,  FN  —  I  transition  in  NO.,. 


Fic;.  6.  Experimental  and  theoretical  spectra  of  the 
J  =  I  — »  0,  An  =  2  transition  in  ND,. 


parison  between  the  experimental  hyperfine  fre¬ 
quencies  and  those  calculated  from  the  constants 
of  Table  3. 

Although  there  has  been  considerable  experi¬ 
mental  and  theoretical  work  on  the  hyperline 
spectra  of  the  isotopic  species  of  ammonia, 
numerous  uncertainties  and  inconsistencies  re¬ 
main.  These  problems  and  the  relation  of  new 
experimental  work  to  them  have  been  discussed 
recently  by  Garvey  ct  at.  (1976).  The  submilli¬ 
metre  capabilities  of  the  spectrometer  discussed 
in  this  paper  allow  the  observation  of  the  7  =  I 
-v  0  transition,  which  is  less  complex  than  the 
spectra  which  fall  at  longer  wavelengths. 

V.  Discussion 

The  deuterium  quadrupole  interaction  in  the 
water  molecule  may  be  employed  to  investigate 
the  electric  field  gradient  tensor  at  the  deuterium 
site.  Specifically, 

^  (,J  =  27  +  3  ?  7(7  +T)  1 " 

where  c/j  is  the  spectral  quadrupole  constant,  g 
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Fig.  7.  Experimental  and  theoretical  spectra  of  the 
7  =  1  -»  0,  Fn  =  0  transition  in  NDj. 


Tabu:  3.  Hyperlinc  constants  of 
7=1-,  Oof  ND., 


Constant 

Value  (kHz) 

U’&Ij)  N 

816. 2±  1.1“ 

Cn 

2.0±1.5 

19.9  ±0.8 

(  I> 

-2. 2  +  0. 8 

<11  D..(DD)  1 1  > 

-0. 13* 

<I1|D„(ND)  ll> 

-0. 15" 

‘Uncertainties  represent  one  standard  dela¬ 
tion  of  the  experimental  data. 

'‘Calculated  from  molecular  geometry. 


Tabu  4.  Hyperlinc  spectrum  of  ND.,  for  7  =  I  -»  0 


Upper 

state 

Observed" 

(kHz) 

Calculated 

(kHz) 

Difference 

(kHz) 

Is  =  0 

2037  7  +  0.8 

2037.7 

0.0 

235 . 3  ±  2 . 9 

235.4 

-0.1 

II 

£ 

214. 3±  1.9 

213.2 

1  1 

181 . 7  ±  3 . 0 

182.2 

-0.5 

- 1002.9+  1.3 

-  1002.9 

0.0 

F*  =  1 

-  1032. 4±  1 .4 

- 1032.5 

-0. 1 

-  1060.8+  1 .4 

- 1060.7 

0. 1 

‘Measured  relative  to  .MW  <88  kll/.  I  nccrtaintics  represent  one 
standard  deviation  of  the  experimental  data. 


refers  to  the  molecular  principal  inertial  axes. 
PK  is  the  rotational  angular  momentum  along  g, 
and  are  components  of  the  electric  field  gra¬ 
dient  tensor.  Combination  of  spectroscopic  data 
from  HDO  and  D.O  allow  complete  solution  of 
V  due  to  the  rotation  of  the  molecular  inertial 
axes  from  one  isotopic  form  of  water  to  the  other. 
As  a  result  of  Laplace's  equation  and  molecular 
symmetry,  the  tensor  may  be  described  by  values 
of  two  of  its  components  and  by  the  location  of 
its  principal  axes  in  the  molecular  plane.  Posener 
(I960)  first  performed  such  an  evaluation  of  the 
electric  field  gradient  tensor  in  water.  Using 
beam  maser  techniques.  Bluyssen  el  al.  ( \9bla.h) 
were  able  to  significantly  improve  this  analysis. 
As  was  pointed  out  by  Bluyssen  cl  al.  ( 1967a,A), 
Posener  erroneously  attributed  the  wrong  direc¬ 
tion  to  the  angle  of  rotation  from  the  O  I)  bond 
to  the  principal  axis  system  of  the  tensor  V.  The 
direction  of  the  angle  of  rotation  has  more  re¬ 
cently  been  incorrectly  reported  in  the  literature 
(De  Lucia  cl  al.  1974)  because  of  an  erroneous 
communication  to  the  authors.  Posener 's  data, 
as  well  as  that  of  Bluyssen  cl  al.  and  the  present 
experiment,  indicate  that  the  principal  axes  of 
the  quadrupole  tensor  arc  rotated  from  the  O  D 
bond  direction  by  approximately  I  clockwise 
about  Posener's  positive  £  axis. 

Figure  8  shows  the  correct  geometry  for  the 
location  of  the  quadrupole  tensor  relative  to  the 
direction  of  the  molecular  bond,  and  Table  5 
gives  the  components  of  V  in  its  own  principal 
axis  system.  These  results  are  calculated  from  all 
available  data  except  that  from  the  744  and  7,4 
states  which  Verhoeven  cl  al.  ( 1968)  have  shown 
to  be  contaminated  by  centrifugal  distortion 
effects.  The  uncertainties  shown  in  Table  5  in¬ 
clude  contributions  from  both  experimental  un¬ 
certainties  and  model  errors.  The  model  errors 


Flo.  8.  Rotation  of  the  quadrupole  tensor  from  the 
017  bond  direction.  The  positive  £,  r  axes  point  into 
the  paper  for  D,.  Different  sign  conventions  and  (or) 
axis  assignments  would  lie  required  for  Z>2. 
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Tabu  5.  Components  of  the  quadru- 
pole  coupling  tensor  in  its  principal 
axis  system 


D.O 

D.S 

V„ 

307.5  ±0.6 

149.0 

y„ 

-  132.6+0.6 

-59  8 

-  174.9  +  0.6 

-89.2 

(1 

-1  16+3' 

-1  .35 

arise  principally  from  neglect  of  centrifugal  dis¬ 
tortion  and  isotope  effects.  These  errors  were 
not  included  in  the  earlier  calculation  of  Bluyssen 
et  at.  (I967cj,/>)  because  their  smaller  data  set 
did  not  provide  adequate  redundancy  for  their 
evaluation. 

Also  included  in  Table  5  are  similar  results  for 
D,S  (De  Lucia  and  Cederberg  (1971).  It  is  inter¬ 
esting  and  reassuring  to  note  the  similarity  both 
in  the  rotation  of  the  principal  axis  of  the  V 
tensor  from  the  bond  direction  and  the  deviation 
from  cylindrical  symmetry  for  both  molecules. 
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Appendix 

The  matrix  elements  of  the  various  Hamil¬ 
tonian  terms  for  the  nuclear  hyperline  inter¬ 
actions  in  NDj  may  be  conveniently  evaluated  in 
the  coupling  scheme: 

J  +  IN  =  Fn 

[A.I]  /(£>,)  +  I(D2)  +  l(Ds)  =  ID 

F»  +  In  =  F 

In  the  equations  to  follow  the  states  are  desig¬ 
nated  by  spccilic  cases  of  the  general  form 
<A7/N/rN/D/rA/|.  For  brevity,  the  Hamiltonian  is 
retained  in  Cartesian  form. 


[A. 2]  ^01 1  FnIdFM  ‘  V(N)- Q(\)  01 1  FNlDFM^  =  ( —  1  )f  N  j  \  \  ^  [  (1 5l2)(eQqj)N 

[A. 3]  <01 1  Fn/dFJW|/(!V)-C(N)- J|0l  1  FNI„FM)  =  (-  l)F"  1 1  j  ^  j  (6)(Cj)n 

The  relationship  between  C'j  and  C  is  considered  by  C'ook  and  De  Lucia  (1971) 

[A. 4]  (o\\Fn  I„'FM  ‘  X  Q(Di)  V(Dl)  Oil  FN/„FM^  = 
for  /„'  =  /„  =  3: 


( _  I)2,n+  i  +f[(2FN'  +  I )(2Fn  +  I)]1  ’ 


I  fN'  F  N  2()  I 
1 3  3  F  f  \Fn 


'  \  2(  (9V 
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or  lo  =  l D  =  l: 


. . '((2/y  +  1X2F„  +  »]"  f"  If,  V  )  1  2? ,rQ‘h>" 


”?  Y  f}{f.  Y  i}^"*** 


or  la  =  1.  If)  -  3: 

(—  l)2fN  + 1  +  ,[(2Fn'  +  IK2Fn  +  I)]"! 

[A. 5]  /oil  FN'/DFM  £  /(Di)-C(Di)' ./  Oil  FN/DFAf^  = 

x  [(2Fn'  +  1H2Fn  +  1H2/d  +  1)(/D  +  D/d]'  2  J  ^  Jr)  |  FN  ^  6)K  j)” 

[A. 6]  ^011  FN'/pFA/  l(N)  D(NDi)  HDi)  Oil  FN/DFM^  =  (-l)'s''"  "" 

r  .  n 

Ifv  fn  i) 

x  15s  2  <1  l|D_..(ND)|l  l> 


F  i  |  l*  1 

x  [3(2Fn'  +  1)(2Fn  +  l)(2/„  +  1)(/D  +  1  )/[>]'  2  |  y*  /.  (  1  1 


<  1 1  |F)_.|  1 1  >  is  the  specialized  form  of  <F  J\D:;\J  J>,  the  spin-spin  interaction  term,  discussed  by 
De  Lucia  and  Cederberg  (1971). 


[A. 7]  /oil  FN7D  FAf  i  £  l(Di)  D(DiDj)  I(DJ)  01 1  Fn/dFaA  - 

\  *•  ij=)  ' 


i*  i 


for  /D'  =  /D  =  3: 


(-11"-*' »[(2FC  +  l*2F»  +  l)],  !  Y  f}{fn  Y  ;}|I8sT4K1I|0„(DD,|M> 


1/2 /Fn'  2  \  )  1  Fn  U(_6)<i||Dr:(DD)|ll> 


|l  I  FI  }Fn  I  21 


for  Id  =  A>  = 

( —  ,  )2Fn  +  f  + 1  [(2F  N#  +  1)(2Fn  +  D] 
for  la  =  I  >  A)  =  3- 

'[(2F,+1X2F»+1)]“!{Y  Y  YY  Y’  5}l  — 3s®)<l'IB..(DD)l"> 
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Abstract— Microwave  measurements  of  the  self  broadening  parameters  of  four  pure  rotational  transitions  of 
H.S  have  been  carried  out  in  the  millimeter  and  submillimeter  wavelength  region.  The  resulting  parameters 
are  (in  MHzlorr  '):  1,,„-  1„„  7.18  ±0.5;  2:.„-2, ,.  6.78±0.5;  2,., -2«.  9.102:0.5:  7.76 i 0.5.  In 

addition.  Anderson  theory  calculations  have  been  carried  out  for  these  transitions  and  are  found  to  be  in  good 
agreement 

INTRODUCTION 

Because  H-S  is  a  light  asymmetric  rotor  without  transitions  in  the  conventional  microwave 
region,  we  have  used  high  resolution  millimeter  and  submillimeter  microwave  techniques'"  to 
measure  the  self-broadening  parameters  of  four  pure-rotational  transitions  of  H;S  in  the  region 
between  150  and  4(H)GHz.  To  the  best  of  our  knowledge,  these  results  are  the  first  experimental 
measurements  of  the  pressure  broadening  parameters  of  H;S.  We  have  also  carried  out  Anderson 
theory"’’  calculations  and  find  them  to  be  in  good  agreement  with  our  experimental  values. 

These  results  are  of  both  theoretical  and  practical  importance.  fLO  is  the  only  other  light 
asymmetric  species  for  which  experimental  line-broadening  data  have  been  reported  and  some  of 
these  data  are  in  substantial  conflict  with  theoretical  predictions.',’  H:S  is  both  a  minor 
constituent  of  the  atmosphere  and  one  of  the  principal  molecular  components  of  the  interstellar 
medium. 

EXPERIMENTAL  TECHNIQUE  AND  RESULTS 

We  have  discussed  in  previous  publications  the  techniques  for  producing  and  detecting 
millimeter  and  submillimeter  microwave  energy."’  Our  pressure-broadening  measurement 
technique  follows  closely  the  double  modulation  scheme  of  Rusk.14’  Figure  I  shows  an  outline  of 
our  apparatus.  Millimeter  wave  energy  in  the  50-60  GHz  region  was  produced  by  an  OKI  55V 11 
klystron.  A  small  amount  of  this  energy  was  compared  with  an  accurate  local  frequency  standard 
and  the  rest  coupled  onto  a  crystal  harmonic  generator.  The  resultant  harmonic  energy  in  the 
150-400  GHz  region  was  transmitted  by  quasi-optical  techniques  through  the  sample  cell  and 
detected  by  a  l.6°K  InSb  photodetector.  The  pressure-broadened  line  acted  as  a  discriminator  for 
the  small  10  kHz  FM  modulation  that  was  applied  to  the  tube.  The  resulting  10  kHz  signal  was 
amplified  by  the  tuned  amplifiers  of  a  PAR  HR-8  lock-in  and  displayed  on  an  oscilloscope.  A 
cathetometer  and  oil  manometer  were  used  to  monitor  pressure.  Although  H3S  is  well  behaved 
and  reaches  equilibrium  pressure  rapidly,  measurements  were  made  both  in  series  of  increasing 
and  decreasing  pressure  and  on  different  days  to  avoid  possible  systematic  effects.  Because 
rotational  transitions  of  light  asymmetric  rotors  are  such  strong  absorbers  in  the  millimeter  and 
submillimeter  spectral  region,  we  used  absorption  cell  lengths  as  short  as  10  cm.  Even  so.  it  was 
necessary  to  include  a  correction  for  the  large  absorption  by'" 

aj  +  ^(«<>/)3-  •  •  U) 

tThis  work  was  supported  by  the  Army  Research  Office.  Grant  Number  DAHC04  74  (J  0034. 
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Table  I  Hydrogen  sulfide  line  broadening  parameters 


Trar.Jl'lon 

Lliv* 

:  r*qu*ncy 
(GHz) 

KxperLu:--n";l  5 

Broa-tenini 

(MHe/Torr) 

•>  If  bn  Jenin  i 
Parameter 
(MHz/T  --) 

ArvJer. .  n 

The  ry  . -Jue 
UvIHr/Tcrr) 

1 1 , 0  *  . 

ItJS 

7.BO 

7.  18  t  C.  5 

d.S8 

■  X  1,1 

c. 

’ 

Old 

’  -  ' 

U.  ■  A  *  U.  ^ 

’ 

_ 

7.  70 ,0.5 

DISCUSSIONS 

Until  this  work.  H?0  was  the  only  light  asymmetric  rotor  for  which  comparison  could  he 
made  of  experimental  and  theoretical  pressure-broadening  parameters.  Unfortunately,  only  the 
microwave  measurements  of  the  3i.a-22.0(4>  and  6, .a- 5s.,1"’  transitions  are  in  good  agreement 
with  the  theoretical  predictions.  Several  i.r.  measurements  all  gave  results  that  are  about  50 % 
higher  than  the  values  of  Benedict  and  Kaplan.'”  The  good  agreement  between  our  microwave 
results  and  theory,  coupled  with  the  fact  that  the  apparent  discrepancies  in  H?0  are  all  based  on 
lower  resolution  i.r.  measurements,  would  indicate  that  Anderson’s  theory  is  capable  of  giving 
good  results  for  pure  rotational  transitions  of  light  asymmetric  species.  Recent  microwave  work 
by  Nere"21  on  the  somewhat  heavier  H;CO  leads  to  the  same  conclusion. 

Tejwani  and  Yenung"”  have  recently  carried  out  Anderson's  theory  calculations  on  H;S  with 
and  without  the  inclusion  of  higher  order  interactions.  Their  results  are  in  general  agreement 
with  ours  (~10%)  although  the  differences  are  substantially  greater  than  the  differences 
between  our  H:0  calculations  and  those  of  Benedict  and  Kaplan. 
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"Forbidden"  Rotational  Spectra  of  Symmetric-Top  Molecules: 

PH3  and  PD3X 

David  A.  Hei.ms  and  Walter  Gokdy 

Department  of  Pliysiis,  Duke  l ' niversity ,  Durham,  Xortlt  Carolina  27706 

A  millimeter-wave  spectrometer  having  a  sensitivity  of  4  X  ME10  cm  1  in  the  2-mnt  region 
has  been  constructed  for  observation  of  extremely  weak  millimeter-wave  spectra  of  gases.  It 
has  been  used  to  measure  J  —>J,K  =  0  —  .1  transitions  in  I’li.i  and  J  — >  J,  K  =  0 «—  3  as 
well  as  K  =  ±1  ±  4  transitions  in  I’D;.  The  and  t'n  spectral  constants  (in  MHz)  are: 

for  I’M-,  Bo  =  133  480.15  ±0.12  and  Co  =  1 1 7  488.85  ±0.16;  for  IT).,,  Ba  =  69  471 .10  ±  0.03 
and  Co  =  58  974.37  ±0.05.  The  effective  ground-state  values  obtained  for  the  bond  angle 
and  bond  length  are:  for  ITT.  ra  (a)  =  t.420,  and  ao(°)  =  93.34..;  for  I'D.,,  r0  (A)  =  1 .41 76 
and  a o(°)  =  93.35.,.  The  corresponding  zero-point-average  values  were  calculated  to  be:  for 
1'Hi,  r,  (A)  =  1.4269,  ±  0.0002  anti  a,(°)  =  93.228,;  for  PD3,  r,  (A)  =  1.4226;,  ±  0.0001 
anti  tt:(a)  =  93.256,  ±  0.(M)4.  For  both  species,  the  equilibrium  values  are  r.  (A)  =  1.4115, 

±  0.0006  and  n,C)  =  93.32s  ±  0.02. 

1.  INTRODUCTION' 

I  he  pure  rotational  transitions  normally  observed  for  symmetric-top  molecules 
conform  to  the  selection  rules  A./  =  ±1,  \K  =  0.  As  a  result  of  their  axial  symmetry, 
there  is  no  off-axis  component  of  the  permanent  dipole  moment  to  couple  a  radiation 
field  with  the  rotation  about  the  symmetry  axis.  Hence,  resonance  radiation  does  not 
normally  induce  changes  in  the  K  quantum  number  which  measures  the  angular  momen¬ 
tum  about  this  axis.  Consequently,  the  moment  of  inertia  for  rotation  about  the  sym¬ 
metry  axis,  and  the  associated  centrifugal  distortion  constants,  cannot  be  obtained  from 
the  normal  rotational  spectra. 

As  a  result  of  theoretical  developments  and  notable  advancements  in  the  experimental 
techniques  for  spectral  detection,  it  has  become  possible  to  detect  rotational  transitions 
in  the  ground  vibrational  state  corresponding  to  changes  in  the  energv  of  rotation  about 
the  symmetry  axis  of  symmetric-top  molecules,  i.e.,  changes  in  the  K  quantum  number, 
l  or  molecules  with  C3,  symmetry,  such  as  PH,  and  PI)3  which  are  reported  here,  the 
observable  transitions’ — forbidden  in  the  normal  theory  correspond  to  K  changes 
of  ±3. 

I  he  possibility  of  A/\  =  ±3  changes  in  molecules  of  (  a,  symmetry  was  apparently 
first  recognized  in  19(>7  by  Hanson  (/).  The  basic  theory  for  such  transitions  was  later 
contributed  by  Watson  (2).  These  transitions  were  first  detected  as  pure  rotational 
spectra  in  1973  by  Chu  and  Oka  (.?),  who  observed  XI  =  U,  K  =  ±1  < —  =F2  transitions 
in  PH:,  occurring  in  the  43.7  to  47.4  GHz  region.  In  1974  they  reported  (4)  similar 

1  This  work  was  supported  by  the  Army  Research  Office,  Grant  PAAG29-77  -G-0007. 
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measurements  of  A J  =  0,  A  =  ±1  <—  ^2  transitions  for  I’D;)  in  the  29.9  to  .50.8  GHz 
region,  and  for  AsH3  in  the  21.2  to  21.8  GHz  region.  Chu  and  Oka  were  helped  in  their 
measurements  by  the  line  frequency  predictions  based  on  the  previous  observation  of 
A (A  —  l)  =  ±3  vibration-rotation  transitions  of  I’ll.,  and  AsH3  by  Maki,  Sams,  and 
Olson  (5).  Likewise,  our  detection  (6)  of  the  higher-frequency  A  =  0  <—  3  forbidden 
transition  of  PH3  and  PD3  was  assisted  by  a  preliminary  prediction  of  their  frequencies 
from  the  results  of  Maki  el  al.  as  well  as  those  of  Chu  and  Oka. 

The  need  for  observation  of  these  higher  A  transitions  described  in  this  paper  was 
indicated  by  Chu  and  Oka  ( 4 ),  who  reported  that  their  spectrometer  was  not  sufficiently 
sensitive  at  the  higher  frequencies  to  detect  them.  Measurements  of  these  transitions 
confirm  the  assignments  of  the  frequencies  measured  by  Chu  and  Oka  and  provide 
increased  accuracy  in  the  values  of  the  spectral  constants,  particularly  those  influenced 
by  the  distortion  caused  by  rotation  about  the  symmetry  axis.  For  observation  of  weak 
transitions  such  as  these  in  the  shorter  millimeter-wave  region,  we  designed  the  excep¬ 
tionally  sensitive  millimeter-wave  spectrometer  described  below. 

The  0-branch  A  A'  =  ±3  transitions  represent  only  one  class  of  forbidden  transitions 
in  symmetric  molecules  included  in  the  general  theory  of  Watson  el  al.  (2,  7)  and  others 
(A,  9).  Oka  (10)  has  recently  written  an  excellent  review  of  the  experimental  and  theo¬ 
retical  developments  of  this  relatively  new  type  of  spectra.  It  should  be  mentioned  that 
forbidden  transitions  of  symmetric-top  molecules  in  vibrationally  excited  states  were 
theoretically  predicted  many  years  earlier  by  Mizushima  and  Venkateswarlu  (11), 
before  sufficiently  sensitive  spectrometers  were  available  for  their  detection.  In  the 
ground  vibrational  state  the  ofl-axis  dipole  moment  which  couples  the  radiation  field  to 
the  rotation  is  induced  by  centrifugal  distortion.  The  selection  rules  are  determined  by 
the  molecular  symmetry.  In  the  molecules  studied  here,  the  rotation  about  the  b  axis 
of  inertia  induces  a  small  dipole  moment  perpendicular  to  the  symmetry  axis  which 
leads  to  a  slight  admixture  of  the  A  levels  and  a  breakdown  of  the  AA  =  0  selection 
rule.  The  threefold  symmetry  about  the  c  axis  limits  this  admixture  to  A  levels  differing 
by  3.  It  is  of  interest  that  centrifugal  distortion  due  to  rotation  about  the  symmetry 
axis  can  likewise  generate  a  small  dipole  component  along  the  symmetry  axis.  This  com¬ 
ponent  is  of  no  particular  consequence  in  molecules  such  as  PHj,  which  already  has  a 
permanent  dipole  moment  along  the  symmetry  axis,  but  it  can  give  rise  to  forbidden 
XI  =  ±1  transitions  in  molecules  such  as  CH4,  which  has  no  permanent  dipole  moment . 
In  fact,  XJ  =  1  rotational  transitions  of  CH4,  SiH4,  and  GeH4  have  been  observed  by 
Rosenberg  cl  al.  (/J  15)  in  the  far-infrared  region.  Forbidden  transitions  of  CH4  in  the 
ground  vibrational  state  have  also  been  detected  in  the  microwave  region  by  Holt, 
Gerry,  and  Ozier  (16).  These  and  other  forbidden  sped ra  of  symmetric-top  molecules 
are  treated  in  Oka’s  review  (10).  He  also  reports  that  he  and  his  associates  have  recently 
detected  AA  =  ±3  transitions  of  \'H3  with  an  infrared  laser  (17). 


1.  EXPERIMENTAL  DETAILS 

The  microwave  spectrometer  used  in  these  measurements  has:  for  the  millimeter- wave 
source,  a  klystron-driven  harmonic  generator;  for  the  absorption  cell,  a  high-(t  Fabry - 
Perot  transmission  cavity  (semiconfocal  arrangement);  and  for  millimeter-wave 
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detection,  an  indium  antimonide  detector  operated  at  1.6  K..  The  power  from  the 
harmonic  generator  is  injected  into  the  cavity  by  way  of  a  small,  rectangular  aperture 
(0.075  X  0.034  in.)  in  the  center  of  the  tlat,  circular  mirror.  The  power  level  in  the 
cavity  is  monitored  by  radiation  from  another  small,  rectangular  aperture  in  the  same 
plate,  the  output  from  which  is  channeled  to  the  indium  antimonide  detector.  The  Q  of 
the  resonant  cavity  mode  as  indicated  by  the  >'o  Ae,  the  resonant  frequency  divided  by 
the  half-width  of  the  Lorentzian  mode  shape,  is  approximately  600  000  at  140  GHz. 

In  operation,  the  resonant  cavity  is  tilled  w  ith  the  absorbing  sample  of  gas  to  a  pressure 
of  about  15  to  25 /rm.  The  klvstron  is  tuned  so  that  the  proper  harmonic  coincides  with 
the  predicted  absorption  frequency,  and  the  cavity  is  adjusted  to  resonate  at.  that 
frequency.  The  klystron  harmonic  is  then  locked  to  the  peak  of  the  cavity  resonance, 
and  the  cavity  is  tuned  through  the  molecular  absorption  line.  The  drop  in  the  moni¬ 
tored  power  resulting  from  the  change  in  the  loaded  (J  of  the  cavity  corresponds  to  the 
detected  spectral  line. 

To  achieve  control  of  the  klystron  frequency  by  the  cavity,  a  2.6  kHz  sine  wave  is 
superimposed  on  the  klystron  reflector  to  sweep  the  harmonic  back  and  forth  across  the 
center  of  the  cavity  resonance.  The  signal  produced  at  the  detector  (with  the  loop 
opened)  changes  phase  by  180°  as  the  klystron  is  tuned  through  the  cavity  resonance. 
This  detector  signal  is  processed  by  a  phase-lock-in  amplifier  which  generates  a  correction 
voltage  with  the  proper  phase-and-time  response  to  lock  the  klystron  harmonic  to  the 
peak  of  the  cavity  resonance.  Thus,  as  the  cavity  is  tuned  by  movement  of  the  curved 
mirror  in  and  out,  the  klystron  harmonic  and  the  cavity  resonance  track  together.  In 
practice,  the  amplitude  of  the  2.6  kHz  stabilizing  signal  is  reduced  to  about  40  mV,  or 
less,  to  keep  most  of  the  harmonic  power  within  the  envelope  of  the  cavity-response 
curve. 


I  k;.  I.  Diagram  of  spec ammeter. 
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ROTATIONAL  ENERGY  LEVELS  OF  PD3 

KiO.  2.  Energy  level  diagram  of  the  first  five  A"  states  for  ./  -  10  of  l’H:i  showing  the  energy  shifts 
caused  hy  the  .1C'  perturbation,  the  A."  -  2  splitting,  the  level  symmetry  under  particle  exchange),  and 
the  nuclear  spin  degeneracy.  The  transition  indicated  l>y  (i  i  is  that  measured  by  I'liu  and  Oka  /  Those 
indicated  by  (ii)  and  fin)  were  measured  in  this  work. 


The  inside  of  the  cavity  is  fitted  with  Stark-tu adulation  plates  which  are  positioned 
well  out  of  the  microwave  field  to  avoid  lowering  the  {].  A  square-wave  generator  is  used 
to  produce  Stark  modulation  of  the  absorption  line  at  KKI  kHz.  As  the  clock-driven 
cavity  tunes  through  the  absorption,  I  he  detector  signal  is  processed  by  a  phase-lock  in 
amplifier  referenced  to  the  100-kHz  signal.  The  first -derivative  line-shape  signal  from 
the  amplifier  is  displayed  on  the  chart  recorder  (see  Fig.  1).  After  exploratory  scans  lo 
establish  the  shape  of  the  line,  the  peak  of  the  cavity  response  is  tuned  to  the  center  of 
the  absorption-line  contour,  and  the  corresponding  klystron  frequency  is  measured  with 
a  frequency  standard  monitored  by  station  \Y\YYB.  Repeated  measurements  are  made 
for  each  direction  of  approach  to  the  line  center.  The  accuracy  achieved  for  the  line- 
frequency  measurements  is  estimated  to  be  1  pari  in  10s.  The  sensitivity  of  our  spectrom¬ 
eter  (  =  4  X  10  10  cm'1  at  1-14  (ill/,  for  a  time  constant  of  1  sec)  enables  us  to  measure 
the  weak  AA."  =  ±4  transitions  with  this  high  precision.  The  factors  limiting  the 
sensitivity  appear  to  be  the  degree  of  klystron  stability,  the  power  available  from  the 
harmonic  generator,  and  the  extraneous  vibrations  of  the  cavity  end  plates. 


3.  THEORETICAL  HACKUROTM) 


As  a  result  of  centrifugal  distortion,  the  rotational  wavefunction  tpj  k  of  a  molecule 
having  (',<„  symmetry  mixes  slightly  with  the  wavefunction  ipj.Kt*  in  the  ground  vibra¬ 
tional  state.  The  term  in  the  Hamiltonian  responsible  for  the  interaction  may  be  ex- 
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pressed  (/.V)  as 

x’  =  (A*r„,f  4)[(J*t  +  J*-)Jt  +  MJ\  +  •/'-)],  0) 

where  T„,:  is  a  centrifugal  distortion  constant  ami  J  ■,  J . ,  and  J.  are  rotational 
operators. 

In  the  absence  of  the  perturbing  .1C'  term,  tile  rotational  levels  are  all  doubly  de¬ 
generate  for  each  value  of  A  (A  (I),  l'igure  2  illustrates  the  effect  of  the  |iert urb.it ion 
on  the  J  =  10  energy  levels  in  I’Da.  The  energy  levels  corres|xmding  to  K  *  2  arc- 
raised  or  lowered,  but  they  remain  degenerate  with  symmetry  li.  The  A  =  3  level  is, 
however,  split  into  .h  and  .1;  com|xinents.  Ihe  +  and  —  signs  indicate  the  parity  of 
the  total  wave-function.  The  symmetry  assignments  .li,  At,  and  E  indicate  the  sym¬ 
metry  of  the  rovibrational  wavefunction  (excluding  spin)  for  an  interchange  of  two 
hydrogens.  The  integer  preceding  the  symmetry  assignment  (.1 1,  .1*,  or  li)  indicates  the 
nuclear  stat  ist  ical  weight  of  the  state.  Ihe  symmet  ry  select  ion  rules  for  ((type  A  A  =  ±3 
transitions  (1*))  are  ±  +-*  ±,  along  with  . 1 1  <— ► . 1 1 .  At,  or  li*-*  li. 

Kach  uniK-rturbed  energy  level  of  a  symmetric-top  molecule  is  normally  doubly 
dcgcncr.it c  for  each  nonzero  value  of  A.  It  is  necessary  to  form  linear  combinations  of 
fj  .  K  and  4>j.  -  k  in  order  to  obtain  wavefunct  ions  which  have  the  projicr  symmetry  for 
an  exchange  of  two  of  the  hydrogens.  The  new  wave-functions  for  PH.i,  j.k  and  tp  j.k, 
are  defined  in  terms  of  t  lie  wavefunct  ions  <pj . .  k  and  i pj.  k  and  t  lie  nuclear  wavefunct  ions 
<«fM,  ada,  etc.  (where  a  and  J  stand  for  spin  up  and  spin  down,  res|>ectively)  as  follows 
(JO) :  for  A  5^  0  or  a  multiple  of  .1, 

i'j.K  =  (1.  <>!)0.mk (t*««  +  ofitxt-  h +  atope*"*1*) 

±t.r.  K(liaa  +  <*fiac  +  aa/ic  1  3)]  (2) 

and 


1 1'  j.k  =  (1, 
for  A  =  0, 


K  (afifi  +  fh<tl3(~’  h  a  +  dtiaC*"1'1*) 

±  i.,  K  (itfifi  +  ffetfc  5"'i,  +  >¥<><■  4’,K  ’)]  i  (3) 


I  poo  -f-  aflat  -f-  aafl  | 


(4a) 


and  for  A"  =  3 


<P '  ./.x 


aaa 


(1  (>!)(l h.  a  =fc  t h.-n) 


Iddd 

flan  -f-  afia  -f-  aafl 
a(i(i  r 


(4b) 


The  similar  wave-functions  for  PI):|  are  more  complicated  because  of  the  spin  value 
/  =  1  for  I). 

The  formulas  giving  the  energy  levels  may  be  expressed  (J/,  for  A  7*  3  as: 


f- 
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These  A  =  3  levels  contain  the  term  ±  hJ(J  +  l)[-/(2  +  1)  —  2][./ (J  +  1)  —  0], 
first  derived  by  Nielsen  and  Dennison  (id)  in  their  explanation  of  the  anomalous  effects 
in  A  =  3  lines  observed  by  early  microwave  spectroscopists.  The  total  splitting  of  the 
A  =  3  levels  is 


Ac  =  [_E~(J,  3)  -  E f  (J ,  3)1/ h 

r  1  1 

=  2 a, - -/(-/+  1 )[-/  (-/  +  1)  -  2][y(y  +  1)  -  0],  («) 

L  /•(-/, «)  -  /•'(/, 3) J 

The  part  of  Eq.  (9)  that  is  enclosed  in  brackets  is  not  a  constant  because  of  the  presence 
of  the  I;(J,  0)  —  E(J,  3)  expression  in  the  denominator.  The  denominator  F(J,  0) 
—  /■  (./,  3)  is  a  power  series  in  J  (J  +  1) : 

T  9 1)K  Djk  ~  9 Ujkk 

1-(J,  0)  -  /•'(./,  3)  =  9(/>  -  O  1  + - —  + - - - J(J  +  1) 

L  B  —  C  B  —  C 


1/ JJK  +  OLjjkK 

- - - JHJ  +  1)-  - 

/>’  -  C 


I-JJJK 

— - J\J+  l)3 

B  -C  J 


(W) 


It  is  not  sufficient  to  retain  only  the  first  term,  9 (B  —  C).  We  found  that  a  simultaneous, 
least -squares  fit  of  the  available  data  for  normal  transitions  and  forbidden  transitions 
for  the  molecular  constants  did  not  converge  properly  unless  the  full  expression  was  used. 
The  line  strengths  of  the  forbidden  SK  =  ±3  transitions  (1,  2)  are  given  by 

.S'  =  1(0", )'-,„(./  =F  A')(./  T  K  -  \)(J  T-  K  -  2) 

X  {J  ±  A  +  1)(./  ±  A  +  2)(J  ±  A  +  3)(2J  +  1),  J(J  +  1)  (11) 

where 

(6'\U  =  0”x  +  h* Trrr;n:.  llh(B  -  O]  (12) 

is  the  effective  dipole  moment  inducing  the  transition.  The  term, 


is  an  induced  dipole  component  which  results  from  admixture  of  the  ground  vibrational 
state  with  excited  vibrational  modes.  This  cont ribution,  first  recognized  by  Watson  (2), 
is  called  vibrational  intensity  borrowing  bv  Oka  (At).  The  last  term  results  from  cen- 
trifugallv  induced  mixing  of  rotational  levels  of  the  ground  vibrational  state.  This 
contribution  to  the  forbidden  transitions,  first  derived  by  Hansen  (I),  is  called  rotational 
intensity  borrowing  by  Oka  (10).  In  this  expression  r„„  is  the  centrifugal  distortion 
constant  in  ,tC'  of  Eq.  (1),  is  the  permanent  dijrole  moment,  and  C  and  B  are  the 
spectral  constants  with  reference  to  the  symmetry  axis  and  the  perpendicular  axis, 
respectively. 

for  the  present  molecules,  the  contribution  of  the  vibrational  intensity  borrowing  is 
small  in  comparison  with  that  of  the  rotational  intensity  borrowing.  For  I’Hj,  Chu  and 
( )ka  (4)  calculated  that  0"x  =  1.6  X  19  "  D  and  that  the  eentrifugally  induced  moment, 
the  second  term  of  Eq.  (12),  is  8.3  X  10  6  D.  The  relative  contributions  of  these  mo- 
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nunts  to  the  forbidden  line  strengths,  Eq.  (12),  is  (1.6 '8.3)*  =  0.037.  Thus,  for  l’Il3  in 
the  ground  vibrational  state,  the  contribution  of  the  vibrational  mixing  to  the  intensities 
of  the  AA  =  ±3  transitions  is  negligible  in  comparison  with  that  from  rotational 
mixing.  Chu  and  Oka  also  estimate  that  the  centrifugal  rotational  intensity  borrowing 
is  the  dominant  term  for  I’Dj. 


4.  RESULTS 

Twelve  A'  =  ()<— 3  transitions  of  1  ‘  H  3  as  well  as  ten  A  =  ±  1  < —  ±4  and  live  A 
=  0<—  3  transitions  in  l’|)3  were  measured.  The  observed  frequencies  are  given  in  Tables 
1  and  It.  Displacement  in  positions  of  lines  of  l’H3  caused  by  the  Xielsen  Dennison 
splitting  (23)  of  the  A  =  3  levels  (It >3  Mil/  for  3  -  10  in  I’ll  d  increased  the  difficulty 
of  their  detection.  This  splitting  does  not  occur  in  the  transitions  observed  by  Chu  and 
t  >ka  (VI.  but  it  has  been  directly  observed  for  I’ll  bv  Davis  cl  ill.  (JVl  in  a  molecular 
beam  resonance  e.x|>eriment.  Location  of  the  A  =  <) * —  3  transitions  would  have  lieen 
even  more  difficult  were  it  not  for  our  access  to  their  molecular  lieam  measurements. 

The  AA  =  ±3  transitions  exhibit  a  rapid  change  in  the  observed  intensity  with 
increasing  value  of 3.  figure  3  illustrates  the  calculated  intensity  profile  (see  Eq.  (14)). 
This  variation  is  due  partly  to  the  increase  in  the  rotat tonally  induced  coupling  moment 
with  3  and  partly  to  the  variation  in  the  population  of  3  states.  At  the  temperature  of 
the  observation,  3(H)  K,  the  intensity  reaches  a  maximum  for  3  =  10.  The  solid  vertical 
bars  in  Fig.  3  correspond  to  transitions  that  we  have  measured.  The  dashed  vertical  bars 
representing  calculated  lines  have  been  included  for  better  indication  of  the  trend, 
f  igure  4  shows  actual  recordings  of  two  lines. 

The  peak  absorption  coefficients  listed  in  Tables  I  and  II  were  calculated  from  the 
formula, 


«(cm  ') 


1.147  X  10 


(  hen ' 

,2i>  (1 - 

\  IkT 


lii'ii  \ 


X  exp ( —  l'.j .k,  k'T)\_(/}-C) 'g s g i .V[ (3 ,  A  =fc  3)  «—  (3,  A)]],  (14) 


where  the  line-strength  factor  .V  is  given  by  Eq.  (11),  T  =  3(H)  K,  g;  is  the  reduced 
nuclear  statistical  weight  factor  (25),  g*  is  the  statistical  weight  factor  for  each  A  level, 
and  Ac  =  20  MHz  is  the  line  breadth  at  a  pressure  of  1  Torr.  This  value  for  the  Ac  is 
the  same  as  that  used  by  Chu  and  Oka  (V)  for  l’H3.  Other  factors  included  in  the 
numerical  constant  are  described  in  Ref.  (25). 

We  had  little  difficulty  in  detecting  the  3  =  3,  A3  =  0,  A  =  0<—  3  transition  which 
has  a  predicted  intensit  v  of  only  7  X  10  10  cm  From  comparisons  of  the  signal  to  noise- 
ratios  for  the  four  weaker  lines  of  I’Hj  and  l’l)3  with  the  calculated  line  intensities  given 
in  Tables  I  and  II,  we  estimate  the  sensitivity  of  our  spectrometer  to  be  approximately 
4  X  10  10  cm  l. 

A  combined  computer  analysis  was  made  of  the  available  frequencies  for  different 
Ivpes  of  transitions  measured  by  us  and  by  others  as  listed  in  Tables  I  and  II.  The 
analxsis  consisted  of  a  least -squares  titling  of  10  molecular  constants  to  the  observed 
frequencies  of  l’ll3  and  0  to  those  of  l’l)3.  The  program  calculates  the  difference  in  the 
energy  levels  l>\  use  of  Eqs.  (5)  (S)  with  assumed  values  of  the  molecular  constants. 
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I  n;.  3.  Calculated  intensities  for  the  K  0  * —  3  lines  observed  in  the  present  work  (solid  lines)  and 

for  some  not  observed  (dotted  lines). 

Transition  frequencies  are  calculated  and  compared  to  the  measured  transition  fre¬ 
quencies  to  generate  new  estimates  of  the  molecular  constants.  This  iterative  process  is 
repeated  many  times  by  the  computer  until  convergence  is  attained.  The  degree  of  fitting 
of  the  complete  frequency  set  is  shown  by  column  d  of  Tables  1  and  II.  The  molecular 
constants  for  I'll.,  and  I’D,,  obtained  from  the  combined  analyses  are  listed  in  Table  III. 
The  constants  li,i  and  r r,x;  for  IT)3  are  rather  well  determined  since  we  were  able  to 
measure  both  the  K  =  dbl  <—  ±4  and  0<—  d  transitions.  These  same  two  constants  for 
1’H.i  are  not  nearly  so  well  determined,  because  of  their  high  statistical  correlation  and 
the  lack  of  data  from  the  K  =  ±1*— ±4  transitions.  The  large  uncertainties  in  the 
relative  values  of  //„  and  the  Txir-.  do  not  appreciably  limit  the  accuracy  of  the  values 
of  the  other  constants  for  PHj. 

There  are  only  two  independent  structural  parameters  in  PHg  or  PDj,  the  bond  angle 
and  bond  length.  With  the  two  measured  moments,  /;,  =  It  8ir/4u  and  /,  =  /;8jt2C’o, 
an  effect ive  zero-jwint  vibrational  structure  can  be  calculated  for  both  species.  The 
effective  structural  parameters  thus  derived  are  listed  in  Table  IV.  They  are  limited  in 


PHj  K  =0*—  3  TRANSITIONS 


I  ic. .  4.  kciordiiij;*  of  die  .1  4  and  J  -  <>,  K  0  <—  d,  lines  for  I’ll,  at  300  K. 


accuracy  by  t he-  uncorreded  inertial  defect.  A  second  structure,  the  zero-point  average 
structure  ( 4 )  can  also  be  calculated  for  each  species  from  the  observed  /hi  and  C'D  values. 
This  structure  was  calculated  by  Chu  and  Oka  for  I'H3.  We  have  recalculated  the  zero- 
point  average  values  for  the  structural  parameters  with  the  /hi  and  C„  values  derived 
from  a  computer  analysis  of  the  measured  ground-state  rotational  transitions,  as 
described  above.  The  resulting  parameters,  listed  in  Table  IV,  are  essentially  the  same 
as  those  obtained  by  Chu  and  Oka;  the  method  used  for  calculation  of  the  average 

TABLE  IV 
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structural  values  is  the  one  they  describe.  Discussion  of  the  significance  or  meaning  of 
the  various  tyj)es  of  structures  and  methods  for  calculations  of  each  tyjK*  is  given  by 
(’■only  and  Cook  (25). 
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Millimeter  and  submillimeter  microwave  spectroscopy  is  shown  to  be  a  powerful  technique  for  the  study  of 
the  fundamental  processes  of  molecular  lasers.  The  calculation  of  absolute  chemical  compositions, 
excitation  parameters,  and  the  observation  of  collisional  deexcitation  are  demonstrated  for  the  HCN  FIR 
laser. 


PACS  numbers:  42.55.Hq,  33.20.Ea,  52.70.Gw 

Although  a  number  of  techniques  have  been  described 
for  the  study  of  the  fundamental  processes  of  molecular 
lasers, 1-5  much  remains  to  be  understood  about  their 
chemistry,  excitation  mechanisms,  and  decay  pathways 
This  letter  reports  the  results  of  experiments  that 
demonstrate  the  unique  capabilities  of  millimeter  and 
submillimeter  microwave  spectroscopy  in  this  field 

In  this  spectral  region,  molecular  laser  plasmas  are 
transparent  except  at  molecular  resonances  and  emit 
insignificant  noise.  This  latter  property  is  directly 
related  to  the  well-known  difficulty  of  producing  radia¬ 
tion  in  the  submillimeter  far-infrared  region  of  the 
spectrum  Equally  important,  the  noise  modulation  of 
the  probing  signal  due  to  plasma  fluctuations  decreases 
rapidly  as  the  signal  frequency  is  increased  from  the 
centimeter -wave  region  into  the  millimeter  and  sub- 


a,\Vork  supported  by  the  t'.S.  Army  Research  Office,  Grant 
No.  DAAG29-77-G-0007. 


millimeter  region.  Rotational  line  strengths  also  in¬ 
crease  as  approximately  1  A1  and  their  transition 
moments  are  orders  of  magnitude  stronger  than  those 
of  the  vibrational  transitions  studied  in  the  infrared 
Because  of  the  widely  tunable  spectral  coverage  of  our 
technique,  studies  of  both  nonradiative  states  and 
species  as  well  as  those  that  participate  directly  in  the 
lasing  process  are  possible  The  sharpness  of  the 
resonances  (Ac  -  1  MHz)  is  such  that  the  different  states 
and  species  are  well  resolved  and  can  be  studied  in¬ 
dependently  Of  equal  importance  is  the  direct  relation 
between  the  experimental  observables  (frequencies, 
absorption  coefficients,  linewidths,  and  lifetimes)  and 
the  dynamics  of  the  molecular  laser  systems. 

Our  experiment  is  shown  in  Fig.  1.  :  Microwave  ener¬ 
gy  in  the  30— 60-GHz  range  is  produced  by  a  reflex 
klystron.  A  small  portion  of  this  signal  is  compared 
with  a  frequency  standard  and  the  rest  matched  onto  a 
crystal  harmonic  generator.  The  harmonic  output  of 
this  generator  (typically  3rd  -  12th  harmonic)  is  trans- 
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Three  studies  of  the  HCN  FIR  molecular  laser  which 


FIG.  X.  Spectrometer  for  the  study  of  laser  plasmas.  illustrate  the  kinds  of  information  that  can  be  acquired 

are  described  below.  CH,CN  was  used  as  the  fuel. 


mitted  through  the  laser  tube  by  quasioptical  techniques 
and  detected  by  a  1,4  K  InSb  photodetector.  The  sensi¬ 
tivity  of  this  technique  is  such  that  states  whose  popula¬ 
tions  are  ~  10'  cm  1  can  be  seen  in  real  time  on  an 
oscilloscope  screen  with  large  signal  to  noise 

It  can  be  shown  that  to  a  good  approximation  the  ab¬ 
sorption  coefficient  o  for  a  linear  molecule  can  be 
written 

o  (cm"1)  (4.  94  X  10'  Ar7\  <1 ) 

where  r  (GHz)  is  the  transition  frequency,  At*  (MHz, 
Tori  )  is  the  linewidth,  /i  (D)  is  the  dipole  moment,  T 
(  K)  is  the  temperature,  and  /•„  is  the  fraction  of 
molecules  in  the  vibrational  state.7  This  relation  as¬ 
sumes  that  the  rotational  states  and  molecular  veloc¬ 
ities  are  in  equilibrium  and  can  be  described  by  a  tem¬ 
perature  T.  For  pure  rotational  spectra  p  and  i\,  are 
known  to  high  accuracy  and  Ar  can  be  accurately  mea¬ 
sured  from  the  high -resolution  microwave  data.  The 


I  IG.  j.  Total  IICN  concentration  derived  from  the  J  o—  t 
rotational  absorption  in  the  ground  vibrational  state. 


The  total  HCN  concentration  is  ar.  important,  but  un¬ 
known,  parameter  of  the  laser.  Chantry8  estimates  it  as 
1  pm  partial  pressure  for  typical  operating  conditions, 
and  Schotzan  and  Kneubuhl1  have  measured  with  a  mass 
spectrometer  the  relative  concentration  as  a  function  of 
discharge  current.  Figure  2  shows  results  calculated 
from  measurements  on  the  J  0-1  tx-ansition  of  the 
00  0  vibrational  state.  In  our  experiment,  the  HCN 
concentration  is  maximum  at  a  remarkably  low  current 
(typically  0.01  A)  and  decreases  substantially  at  higher 
currents. 

The  well-known  HCN  laser  mechanism  of  Lide  and 
Maki  depends  upon  the  near  degeneracy  of  J -  10  in  the 
11 ’0  and  04n0  vibrational  states  and  a  population  inver¬ 
sion  between  them  in  the  laser  plasma.  This  inversion 
has  often  been  attributed8’’  to  the  metastable  nature  of 
the  10"0  state  and  the  bending  stack  built  upon  it  Fig¬ 
ure  0  shows  the  absolute  populations  of  the  lO^O  and 
00 ’1  states  as  a  function  of  the  discharge  current.  It 
should  be  noted  that  these  states  are  essentially  equally 
populated  and  that  the  “metastable”  nature  of  10n0  is 
thus  shown  not  to  be  a  requirement  for  a  population 
excess. 


PRESSURE  (TORR) 


FIG.  I.  Pressure  dependence  of  the  excitation  of  the  lo°o  and 
00°1  vibrational  states  of  HCN. 
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!'.  u r i‘  4  shows  the  population  of  10  0  and  00  1  as  a 
tu  ict it >n  <il  pressure.  These  results  are  consistent  with 
a  model  in  which  resonant  deexcitation  between  the 
iO.m.Oi  and  .l.«  -3.0)  bending  stacks  is  an  important 
limit im;  factor  on  the  operating  pressure  of  the  HCN 
laser.  On  the  other  hand,  the  10.//.1)  stack  is  not  in 
close  resonance  with  another  family  of  states  and 
maintains  its  excitation  at  substantially  higher 
pressures. 

In  conclusion,  we  have  demonstrated  the  kinds  of 
information  that  can  be  obtained  in  a  very  direct  manner 
from  our  experiments.  It  is  clear  that  they  have  appli¬ 
cability  both  to  other  states  and  species  found  in  the 
HCN  laser  and  to  molecular  lasers  in  general.  These 
data  provide  both  constraints  and  tests  for  models  of 
the  complex  dynamics  of  these  systems  as  well  as 
suggestions  for  new  models. 


The  author  would  like  to  thank  the  Research  Corpora¬ 
tion  for  its  support  of  this  work. 
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NOTES 

Centrifugal  Distortion  in  Water  and  Hydrogen  Sulfide1 

It  is  the  intent  of  this  note  to  present  specific  results  pertaining  t«»  the  changes  in  molecular  geometry 
of  water  and  hyd rotten  sulfide  u  hich  arc  induced  by  centrifugal  effects  in  rotational  states  of  the  molecules. 
Hydrogen,  deuterium,  and  tritium  substitutions  are  considered.  A  general  ex  predion  tor  ccntritugally 
induced  geometry  changes  has  been  presented  by  Toyam  t  ct  <il.  1 1)  liut  specific  results  included  only 
the  s>  mmctric  forms  of  water,  115 )  and  IM).  A  sign  error  i<  pre>ent  in  that  work  but  appears  to  be  only 
typographical,  as  the  numerical  results  given  are  correet.  '1  he  term  n  sin  u  ir.  Ecj.  (.55  should  be  negative. 
Considerable  simplification  occurs  in  the  symmetric  cases  and  a  calculation  of  centrifugal  distortion  in 
111)0,  I)TS,  etc.,  must  proceed  from  the  most  general  approach. 


TABLE  I.  Coefficients  of  Eq.  (3)  for  As  ter  and  Hydrogen  Sulfide 


P 

Ar/r 

A? 

A  r/r 

T 

h20 

a 

1.016 

-10.742 

H2S 

0.300 

-2.878 

b 

0.293 

3.041 

0.212 

2.238 

c 

0.247 

0.000 

0.126 

0.012 

D?0 

a 

0.565 

-5.976 

d2s 

0.159 

-1.525 

b 

0.147 

1.521 

0.106 

1.120 

c 

0.126 

0.070 

0.064 

0.025 

T2° 

a 

0.415 

-4.390 

V 

0.112 

-1.075 

b 

0.093 

1.016 

0.071 

0.748 

c 

0.086 

0.090 

0.044 

0.029 

P 

Vi 

■Vr2 

At 

irl/r! 

V? 

is 

HDD 

a 

0.969 

0.673 

-9.192 

HDS 

0.372 

0.033 

-2.639 

b 

0.065 

0.293 

1.483 

0.003 

0.179 

0.705 

c 

0.111 

0.204 

0.030 

0.046 

0.088 

0.013 

HTO 

a 

0.955 

0.60G 

-9.160 

HTS 

0.374 

0.075 

-2.855 

b 

0.028 

0.232 

0.857 

0.002 

0.125 

0.370 

c 

0.066 

0.170 

0.034 

0.026 

0.0/4 

0.012 

DTO 

a 

0.541 

0.449 

-5.334 

DTS 

0.139 

0.071 

-1.359 

b 

0.072 

0.157 

1.117 

0.021 

0.116 

0.632 

c 

0.084 

0.117 

0.076 

0.037 

0.054 

0.022 

aAll  coefficients  have  been  r.ultiplieJ  by  10  ;  antes  are  in  radians;  r^  is 
associated  with  the  lighter  atori  in  the  * : yrretric  species. 
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The  general  expression  fur  the  tlistorliott  effects  are  given  in  Kef.  ( / )  in  matrix  form 

(R )=,-„<  -  F*  'G,r‘B-|.\,  (1) 

w  here 

<R)reti.  =  (  Ar,  )  (2) 

\(rir2)'Au/ 

and  \  is,  in  the  case  of  the  bent  triatomic  molecule,  a  nine  component  vector  describing  the  position  ot 
the  three  nuclei  in  a  body  fixed  Cartesian  coordinate  system  which  is  parallel  to  the  molecular  inertial 
axes.  Fa.  Ga,  B,  and  ■!»  are  defined  in  Kef.  (/).  As  was  shown  by  tile  authors  of  Kef.  (/),  it  is  convenient 
to  rewrite  lap  (1)  as 

=  p*U‘J,'  +  ps((V'  +  pt(/V.',  (o) 

where  a,  b.  and  c  designate  the  molecular  inertial  axes  and  (/V,  etc.,  are  the  expectation  values  of  the 
angular  momentum  component  opera  tors  along  the  molecular  equilibrium  axes.  1  hese  values  max  be 
easily  calculated  for  any  specific  rotational  level  of  a  given  molecule  (2).  The  vector  coefficients  p„. 
Pi,,  and  Pc  of  lap  pi),  which  are  defined  in  Kef.  t / 1 .  have  been  evaluated  and  are  presented  in  Table  1 
for  water  and  hydrogen  sulfide.  Coclticients  for  ICO  and  1M)  have  been  recalculated  for  completeness. 

The  force  constants  used  are  those  of  Cook,  I)e  l.ucia.  and  llelminger  (3.  7b  As  noted  by  the  authors 
of  Kefs.  (J)  and  (7),  model  errors  produce  significant  uncertainties  in  the  values  of  the  force  constants. 
The  constants  used  are  those  obtained  from  infrared  data  w  ith  equal  weighting.  These  uncertainties  are 
the  primary  limitation  to  the  accuracy  of  the  coclticients  of  this  work  and  those  of  Ref.  (/ ' .  Symmetric 
forms  of  water  and  hydrogen  sulfide  distort  symmetrically  and  only  one  coefficient  is  necessary  to  de¬ 
scribe  the  change  in  bond  length.  In  the  case  of  the  asymmetric  forms,  changes  in  both  bonds  must  be 
specified.  In  Table  I.  n  refers  to  the  bond  extending  from  the  apex  atom  to  the  atom  of  least  mass. 
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Millimeter  spectrum  and  molecular  constants  of  silicon  monoxide* 
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A  new  high  lempcralure  millimeter-wave  microwave  spectrometer  has  been  constructed  and  used  to  measure  a 
wide  range  of  vibrational  and  rotational  states  of  silicon  monoxide.  This  work  results  in  accurate  rest 
frequencies  for  all  of  the  intersteller  StO  maser  transitions  that  have  been  observed  as  well  as  accurate 
measurements  or  predictions  of  all  transitions  that  are  likely  to  be  of  astrophysical  interest  In  addition,  the 
Dunham  spectral  and  potential  constants  are  calculated  for  the  three  major  isotopic  species  For  ’Si  '<) 

V  21  7X7  453(11)  MM/.  V,,  151020(11)  Mil/.  )\  70.5(24)  kll/.  V„,  293KUJ)  kll/. 

u,  5  300(22)  \  10’  cm  a,  2.0X00(4 1  >.  a  5  7(7).  a.  0  (X5)  The  calculated  equilibrium  paiametcrs 
for  ‘Si  "O  are  H,  21  7X7  5(10)  Mil/,  u),  1252.(3)  cm  .  ui.x,  5  0b<7l)cni  '.  and  r,  I  500  73(4)  A 


INTRODUCTION 

III  1968.  Torring1  reported  the  laboratory  mea¬ 
surement  of  the  <7=0-1  rotational  transitions  of 
2aSi160  in  the  43-GHz  region  for  four  vibrational 
states.  r  =  0, 1,2.3.  Meanwhile,  a  maser  signal 
at  86  243  MHz.  observed  by  radioastronomers2 
in  the  nebular  source  Orion  A.  has  been  assigned 
to  the  <7=2-1.  /•  =  1  transition  of  28SilfiO  by  Lovas 
and  Johnson2  from  a  theoretical  prediction  of  the 
frequency  based  on  Torring’s  measurement  of  the 
<7=0-1  transition.  Since  that  time,  SiO  maser 
emission  lines  corresponding  to  the  transitions 
i  =l,  <7=1-0;  e  =  l,  <7=2-1;  r=l,  <7=3-2;  and 
r  =  2,  <7=1-0  have  been  observed  in  20  astronom¬ 
ical  sources.1 

With  the  present  work,  direct  laboratory  mea¬ 
surements  have  been  made  of  the  frequencies  of 
all  of  the  transitions  which  give  rise  to  the  im¬ 
portant  maser  signals  from  outer  space,  and  mea¬ 
surements  of  the  rotational  spectrum  of  28Si,60 
are  extended  to  303  926.96  GHz.  These  measure¬ 
ments.  shown  in  Table  I,  include  seven  rotational 
transitions  of  28Si160  in  five  vibrational  states. 
They  allow  precise  evaluation  of  the  centrifugal 
distortion  effects  which  influence  all  rotational 
and  vibrational  constants  of  the  molecule  and 
which  cannot  be  evaluated  from  the  <7=0-1 
transition  alone.  These  constants  make  possible 
the  accurate  prediction  of  all  transitions  that  are 
likely  to  be  of  astrophysical  interest.  Figure  1 
shows  the  laboratory,  real-time  oscilloscope 
display  of  the  most  commonly  observed  astro- 
physical  maser,  the  <7=2-1,  r  =  1  transition. 

LXPIRIWFNT 

A  sufficient  concentration  of  silicon  monoxide 
vapor  for  microwave  absorption  measurements 
was  produced  by  the  heating  of  a  homogeneous 


mixture  of  silicon  and  silicon  dioxide  to  1350  C 
in  a  high- temperature  reaction  chamber  which 
also  served  as  the  microwave  absorption  cell. 

This  combination  of  reaction  chamber  and  absorp¬ 
tion  cell  consisted  of  a  Mullite  MV30  ceramic 
tube  (9x45  cm)  concentrically  positioned  inside 
a  cold- rolled  steel  pipe  (25  x  105  cm)  which  was 
connected  to  a  pumping  line.  For  the  purpose  of 
mechanical  support  and  thermal  insulation,  the 
volume  between  the  ceramic  tube  and  the  outer 
steel  jacket  was  filled  with  A.  P.  Green  G-26 
insulating  firebricks.  The  heating  element  was 
constructed  from  ten  graphite  welding  rods  which 
were  electrically  connected  in  series  and  equally 
spaced  around  the  inside  of  the  Mullite  tube. 

A  tapered  horn  combined  with  a  Teflon  lens 
focused  the  microwave  radiation  into  the  reaction 
chamber  and  through  the  vapor  above  the  re¬ 
actants.  The  radiation  emerging  from  the  reac¬ 
tion  zone  was  then  focused  onto  a  sensitive  InSb 
radiation  detector  operated  at  1.6  K.  The  detected 
absorption- line  signals,  after  amplification,  were 
displayed  on  an  oscilloscope  for  the  frequency 
measurements.  Description  of  the  millimeter- 
wave  spectrometer  and  measuring  techniques  are 
given  in  an  earlier  publication.5 

DATA  AND  ANALYSIS 

The  most  accurate  formulas  for  analysis  of  the 
rotational  spectra  of  diatomic  molecules  are  pro¬ 
vided  by  Dunham  theory.6  Dunham  solved  the 
Schrodinger  wave  equation  for  the  diatomic  vi¬ 
brating  rotor  with  a  potential  function  of  the  form 

F(t) =  chn0t,2(\  +  fJt4  +  ti242+  •  •  • )  .  (1) 

where  £  =  (r -  re)/ra,  n0=w’2/4Bc.  and  nltn2, . . . 
are  the  potential  constants.  His  solution  for  the 
rotation- vibration  energies  is  expressed  by  the 

223 


:z  t 


M  A  N  S  O  N  ,  JR..  CI.AKK.  I)  K  i,  I  CIA.  A  IS  I)  G  (>  K  l>  V 


15 


TABLK  1.  Observed  transitions  of  silicon  monoxide. 


J—J  +  l 

V 

Measured 

frequencies 

(MHz) 

Difference  from 
calculated 
frequencies 
(MHz) 

2Jsi  ICq 

0—  1 

0 

43  423.76“ 

-0.04 

1 

43  122.03“ 

0.00 

2 

42  820.48“ 

-0.06 

3 

42  519.34“ 

0.01 

1  —2 

0 

86  846.96 

0.07 

1 

86  243.37 

0.02 

2 

85  640.46 

0.09 

2  —  3 

0 

130  268.61 

0.04 

1 

129  363.24 

-0.02 

3—4 

0 

173  688.31 

0.16 

1 

172  481.15 

0.09 

2 

171  275.28 

0.18 

3 

170  070.35 

0.07 

4  —  5 

0 

217  104.98 

0.09 

1 

215  595.95 

-0.09 

9 

214  088.54 

-0.06 

3 

212  582.60 

0.04 

4 

211  077.87 

-0.07 

5  —  U 

0 

260  518.02 

-0.09 

1 

258707.45 

-0.04 

C>  —  7 

0 

303  926.96 

-0.15 

-%i"’o 

0—  1 

0 

42  879.82  a 

0.00 

MSi,60 

0—  1 

0 

42  373.34“ 

0.00 

“Reference  t. 


series 

=  *  £  K  +  '•)V(^+1)I".  (2) 

I,  m 

in  which  the  coefficients  K,_ m  are  now  known  as 


FIG.  1.  Video  display  of  the  J  f  —  2,  v  =1  transition 
of  ;8Si,sO. 

Dunham’s  constants.  Each  F,.m  is  a  function  of 
the  molecular  parameters  re,a0,al,a2, . . .  which 
appear  in  Eq.  (1).  In  terms  of  the  familiar  spec¬ 
troscopic  constants, 

^  io ~  *  Yn~Bt.  1  ii  »  ^2i  ~ Y.  ' 

1  20  '  —  ^  02  ~  —  ^  12  ~  ~  »  ^03  ~  ^ e  ■ 

With  the  selection  rules  J  —  J+l  and  r  —  r  for  pure 
rotational  transitions,  Eq.  (2)  predicts  the  rota¬ 
tional  frequencies  to  be 

Vv,  2-m.i  =  2[Kot  +  YuO’  +  z>  +  l,2i('’+  i)2  +”*l(</+  1) 

+  4Y02(J  +  l)3*  •  •  •  ,  (4) 

where  only  terms  which  give  a  detectable  contri¬ 
bution  to  the  spectra  reported  in  this  paper  are 
retained.  For  the  less  abundant  isotopic  species 
the  Dunham  isotope  relation. 

V’,,m  =  (n/n')<U2l"),1Y,.„.  (5) 

is  accurate  to  within  experimental  error  for  all 
F,t  m  except  Vol.  In  Eq.  (5)  p'  denotes  the  reduced 
mass  of  the  less  abundant  species,  and  the  unprimed 
parameters,  V,  m  and  p ,  refer  to  the  parent  species 
(2"SileO). 

The  silicon  monoxide  microwave  spectrum  was 
analyzed  by  substitution  of  each  measured  ab¬ 
sorption  frequency  and  the  appropriate  J  and  r 
quantum  numbers  into  Eq.  (4).  In  this  manner  a 


TABLE  II.  Comparison  of  molecular  constants  measured  in  this  work  with  previous  values. 


Constants 

28Si 

This  work* 

,6o 

Others 

29Si  ,gO 

This  work  Others 

^Si 

This  work 

I6q 

Othe  rs 

Y oi  (MHz) 

21  787.453(11) 

21  787.462(25)” 

21  514.046(13) 

21  514.074(25)” 

21  259.494(13) 

21  259.425(25)” 

Y,,  (MHz) 

-151.026(11) 

-151.05(4)” 

-148.193(11) 

-148.22(4)” 

-145.570(11) 

-145.59(4)” 

r21  (kHz) 

70.5(24) 

76.(1)” 

68.7(24) 

74.(1)” 

67.1(24) 

72.(1)” 

F02  (kHz) 

-29.38(13) 

-29.9(1)”” 

-28.64(13) 

-29.2(1)”'” 

-27.97(13) 

-28.5(1)”'” 

He  (MHz) 

21  787. 5(1  d) 

21  514.1(10) 

21  259.5(10) 

we  (cm-1) 

1  252.(3) 

1  241.4(1) 

1  244.(3) 

1  237.(3) 

U),x,  (cm"1) 

5.96(71) 

5.9(1) 

5.89(71) 

5.82(71) 

“Numbers  in  parentheses  represent  one  standard  deviation. 
’’Reference  1. 
c  Reference  7. 
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TABLK  III.  Potential  constants  and  equilibrium  bond 
length  for  silicon  monoxide. 


This  work* 

Ref.  1 

«0 

5.390(22)  x  10s  cm'1 

5.298(1)  x  10s  cm"1 

"I 

-2.9899(41) 

-2.973(2) 

a  2 

5.7(7) 

5.58(10) 

«3 

-9.0(5) 

-3.32(75) 

1.509  733(40)  A 

1.509  732(40)  A 

‘Errors  are  in  standard  deviation. 


set  of  .V  simultaneous  equations  in  six  unknowns 
was  constructed,  where  is  the  number  of  ex¬ 
perimentally  measured  frequencies  and  where  the 
six  unknown  constants  are  Vol,  Vu .  >'21 ,  and  Y0., 
for  28SilliO,  plus  the  K01  values  for  2!>SilB0  and 
30Sil6O.  By  the  method  of  least  squares,  this  set 


of  equations  was  solved  for  these  constants.  The 
vibrational  constants  Yl0  and  Y20  were  calculated 
for  each  isotopic  species  by  use  of  the  following 
relations: 


v  ^  (4>'3  /_  >'  )l/a 

1  10  1  01'  *  02' 

(6) 

^  20  =  2  ^  0l(^2  “  T  *0  * 

(7) 

where  ft,  and  a.,  are  the  Dunham  potential  con¬ 
stants  which  appear  in  Eq.  (1). 

RESULTS  AND  COMM1  NTS 

Table  II  summarizes  the  results  of  this  analysis 
and  gives  a  comparison  of  the  rotational,  vibra¬ 
tional,  and  potential  constants  derived  from  our 
measurements  with  those  obtained  by  others. 


TABLE  rv 

Predicted  rotational  spectrum 

>f  silicon 

monoxide. 

28Sil60 

J9sil6o 

MSil60 

Transition 

Predicted 

Predicted 

Predicted 

V 

frequency  Uncertainty 

frequency  Uncertainty 

f  requency  l  ’nee  rtainty 

0 

0—  1 

43  423.798 

0.015 

42  879.820 

0.025 

42  373.340 

0.025 

1 

0  —  1 

43  122.027 

0.015 

42  583.710 

0.028 

42  082.409 

0.028 

•) 

0—  1 

42  820.539 

0.015 

42  287.874 

0.030 

41  791.800 

0.030 

3 

0—  1 

12  519.332 

0.015 

11  992.314 

0.030 

41  501.531 

0.030 

1 

0—  1 

42  218.407 

0.020 

11  097.028 

0.035 

41  211.405 

0.034 

5 

0  —  1 

11  917.701 

0.052 

II  402.018 

0.053 

40  921.007 

0.052 

0 

1  —2 

80  840.891 

0.028 

85  758.953 

0.050 

84  740.009 

0.050 

1 

1—2 

80  243.350 

0.020 

85  160.732 

0.050 

84  104.200 

0.050 

2 

1  —2 

85  040.373 

0.022 

84  575.001 

0.000 

8.3  583.000 

0.000 

3 

1  —2 

85037.959 

0.025 

83  983.940 

0.000 

83  002.391 

0.000 

4 

1  —2 

84  430.110 

0.052 

S3  393.309 

0.070 

82  422.259 

0.009 

5 

1  —2 

83  834.824 

0.103 

82  803.348 

0.100 

SI  842.003 

0.103 

0 

2  —  3 

130208.574 

0.037 

128  030.710 

0.070 

127  117.335 

0.070 

1 

2  —  3 

129  303.202 

0.027 

127  748.379 

0.085 

120244.721 

0.085 

2 

2  —  3 

128  458.790 

0.030 

120  800.873 

0.092 

125  372.912 

0.092 

3 

2  —  3 

127  555.170 

0.034 

125  974.192 

0.092 

124  501.908 

0.092 

4 

2  —  3 

120  052.402 

0.075 

1 25  088 .335 

0.105 

123  031.710 

0.104 

5 

2  —  3 

125  750.473 

0.153 

124  203.303 

0.158 

122  702.317 

0.154 

0 

3  —  4 

173  088.142 

0.046 

171  512.400 

0.104 

109  480.047 

0.104 

1 

3  —  4 

172  481.000 

0.035 

170  327.904 

0.118 

108  323.102 

0.117 

•) 

3  -  4 

171  275.105 

0.041 

109 144.023 

0.127 

107  100.750 

0.127 

3 

3—  4 

170070.278 

0.045 

107  902.381 

0.126 

105  999.411 

0.120 

1 

3  —  4 

108  800.579 

0.098 

1 60  781 .239 

0.142 

104  839.147 

0.140 

5 

3  —  4 

107  004.008 

0.202 

105  001 .190 

0.210 

103  079.957 

0.205 

0 

4  —  5 

217  104.891 

0.058 

214  385.351 

0.138 

21 1  853.275 

0.137 

1 

4  —  5 

2 1 5  590.037 

0.052 

212  904.800 

0.157 

210  398.918 

0.156 

2 

4  —  5 

214  088.594 

0.061 

21  1  425.023 

0.108 

208  945.902 

0.107 

3 

4  —  5 

212  582.500 

0.002 

209  947.820 

0.166 

207  494.230 

0.105 

4 

1  —  5 

211  077.930 

0.122 

208  471.392 

0.183 

200  043.899 

0.181 

5 

4  —  5 

209  574.722 

0.249 

200  990.339 

0.204 

204  594.912 

0.258 

0 

5—  0 

260  518.114 

0.084 

257  254.800 

0.181 

254  210.540 

0.180 

1 

5  —  0 

258  707.489 

0.086 

255  478.198 

0.200 

252  471.317 

0.204 

2 

5  —  0 

250898.557 

0.095 

253  703.185 

0.219 

250  727.099 

0.217 

3 

5  —  (i 

255  091.317 

0.094 

251  929.823 

0.210 

248  985.092 

0.214 

1 

5  —  0 

253  285.708 

0.152 

250  158.109 

0.232 

247  245.290 

0.228 

5 

5  —  0 

251  481.912 

0.298 

248  388.040 

0.322 

245  500.510 

0.315 
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Torring‘  could  not  determine  the  effects  of  cen¬ 
trifugal  distortion  from  his  microwave  measure¬ 
ments  alone;  but,  by  combining  his  accurate  mea¬ 
surements  of  the  J=  0—  1  transition  with  the 
earlier  measurements  of  the  vibrational  spectra 
of  SiO  by  Lagerqvist  and  Uhler,7  he  was  able  to 
make  corrections  for  these  effects  in  his  derived 
constants.  We  believe  that  the  vibrational  con¬ 
stant  <j)e,  which  we  have  derived  entirely  from 
microwave  measurements  (see  Table  II).  is  more 
accurate  than  the  values  obtained  more  directly 
from  vibrational  spectra.  The  inclusion  of 
Torring's  frequencies  for  the  <7=0—1  transitions 
with  our  results  did  not  alter  the  constants  de¬ 
rived  from  our  results  alone.  We  conclude  that 
his  measurements  are  at  least  as  accurate  as 
ours. 

Table  III  gives  the  Dunham  potential  constants 
aa,al.n2,a2  which  were  calculated  from  the 
J 'l  m 's  of  Table  II.  The  Yu  constant,  which  is 
necessary  for  the  computation  of  the  constants 
a2  and  aK  but  not  recorded  in  Table  II.  was  found 


to  be  zero  within  the  limits  of  experimental  er¬ 
ror.  A  contribution  to  the  uncertainty  in  the  at 
from  a  calculated  upper  limit  for  Yl2  is  included 
in  our  estimates  of  the  uncertainty  in  the  This 
contribution,  as  well  as  other  contributions,  was 
neglected  in  Torring's  uncertainties.  Although 
our  listed  uncertainties  are  somewhat  larger 
than  those  quoted  for  previous  work,  our  «4  val¬ 
ues  are,  in  fact,  substantially  more  accurate. 

The  values  for  Bc,  which  appear  in  Table  II  were 
obtained  from  Yol  by  application  of  the  Dunham 
correction.  The  wobble- stretch  correction  is 
zero  to  within  experimental  error,  but  a  calcu¬ 
lated  upper  limit  on  it  is  included  in  the  uncer¬ 
tainty.5 

Table  IV  shows  a  set  of  predictions  of  all  tran¬ 
sitions  that  appear  to  be  of  possible  astrophysical 
interest  as  well  as  the  uncertainties  in  these  pre¬ 
dictions.  For  the  higher  i\J  states,  these  pre¬ 
dictions  are  comparable  in  accuracy  with  the  ex¬ 
perimental  observations,  but  at  low  v,J  they  are 
substantially  better. 
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Millimeter  and  submillimeter  microwave  techniques  have  been  used  for  the  spectroscopic  study  of  an  HCN 
laser  plasma.  Forty-seven  rotational  transitions  in  12  excited  vibrational  states  have  been  observed 
Numerous  rotational,  vibrational,  and  perturbation  parameters  have  been  calculated  from  these  data  A 
discussion  of  experimental  techniques  is  included. 


I.  INTRODUCTION 

The  selective  excitation  of  vibrational  states  in  dis¬ 
charges  or  in  the  products  of  discharges  is  the  basis  for 
a  number  of  important  molecular  lasers.  Calculations 
based  on  oscillation  threshold  conditions  show  that  these 
highly  excited  states  are  sufficiently  populated  in  active 
laser  plasmas  for  study  by  the  techniques  of  quasioptical 
millimeter  and  submillimeter  microwave  spectroscopy. 
Most  molecular  lasers  are  active  on  small  molecular 
species  with  large  rotational  constants  and  favorable 
partition  functions.  As  a  result,  spectral  coverage 
throughout  the  millimeter  region  and  into  the  submilli¬ 
meter  region  is  required  if  for  no  other  reason  than  to 
reach  the  transitions  of  interest.  Microwave  studies  of 
these  highly  excited  vibrational  states  are  uncommon  due 
to  their  vanishingly  small  populations  at  room  tempera¬ 
ture. 

Many  interesting  and  complex  vibration-rotation  in¬ 
teractions  become  important  at  microwave  resolution  and 
can  be  studied  by  these  techniques.  In  addition,  they  are 
a  powerful  means  for  the  study  of  the  physics  of  the  com¬ 
plex  molecular  laser  systems.  We  have  given  a  pre¬ 
liminary  report  on  a  study  of  the  HCN  laser1  and  a  full 
discussion  will  be  published  elsewhere. 

HCN  has  been  extensively  studied  in  the  infrared  and 
has  been  the  prototype  molecule  for  investigations  of 
vibration-rotation  interactions  in  linear  molecules.2-6 
Rotational  transitions  in  the  ground  vibrational  state 
were  among  the  first  to  be  observed  by  millimeter  mi¬ 
crowave  spectroscopy1  and  transitions  in  low  lying, 
thermally  populated  vibrational  states  have  also  been  in¬ 
vestigated. 5,6,8  In  this  paper  are  reported  spectroscopic 
studies  in  the  millimeter  and  submillimeter  microwave 
region  of  12  excited  vibrational  states  of  HCN.  Included 
are  the  11*0  and  04°0  states  which  are  the  basis  of  the 
well  known  laser  emission  in  the  far  infrared. 

II.  EXPERIMENTAL 

The  millimeter  and  submillimeter  wave  region  of  the 
spectrum  is  in  many  ways  ideal  for  the  study  of  the 


11  This  work  was  supported  by  the  U.  S.  Army  Research  Office, 
Grant  No.  DAAG29-77-G-0007. 

“Travel  to  Duke  University  supported  by  the  Southern  Regional 
Education  Board  of  Atlanta. 


molecular  species  contained  in  laser  plasmas.  A  corol¬ 
lary  of  the  well-known  difficulty  of  producing  radiation 
in  this  spectral  region  is  that  laser  plasmas  emit  almost 
no  interfering  radiation.  It  can  also  be  shown  that  the 
noise  modulation  of  the  microwave  signal  due  to  fluctua¬ 
tions  in  the  index  of  refraction  which  result  from  plasma 
instabilities  decreases  systematically  with  increases  in 
frequency.  In  addition,  at  lower  microwave  frequencies, 
the  quasioptical  techniques  that  allow  the  remote  sam¬ 
pling  of  the  plasma  arediff  icult  and  rotational  line 
strengths  decrease  as  approximately  i/s.  At  infrared 
wavelengths  the  vibrational  transitions  are  orders  of 
magnitude  weaker  than  the  strong  pure -rotational  transi¬ 
tions  studied  in  this  work,  and  the  emission  of  noise  in¬ 
creases  sharply.  It  is  interesting  to  note  that  other 
spectroscopic  techniques  have  been  unable  to  establish 
the  existence  of  HCN  in  any  vibrational  state  in  the  laser 
plasma.  * 

Figure  1  shows  the  details  of  the  system  used  for  this 
study.10  Energy  in  the  millimeter  and  submillimeter  re¬ 
gion  is  produced  by  klystron  driven  crystal  harmonic 
generation.  The  output  of  this  generator  is  focused  by 
quasioptical  techniques  through  the  laser  plasma  and 
detected  by  an  InSb  photodetector  which  operates  at 
1.  4  °K.  Because  only  the  stretching  modes  are  pref¬ 
erentially  excited  in  a  laser  discharge, 1  the  lower  en¬ 
ergy  bending  modes  were  thermally  populated  in  a  cell 
at  600  °K.  The  sensitivity  of  this  technique  is  such  that 
all  studies  reported  in  this  paper  are  either  the  result 
of  real  time  measurements  on  an  oscilloscope  screen  or 
lock-in  techniques  with  r  «  1  sec. 

III.  THEORY 

In  principle,  the  energy  levels  and  rotational  transi¬ 
tion  frequencies  of  HCN  can  be  calculated  by  construc¬ 
tion  and  diagonalization  of  the  vibration-rotation  Hamil¬ 
tonian.  Although  this  procedure  leads  directly  to  the 
unperturbed  spectral  constants,  it  requires  iterative 
nonlinear  fitting  and  makes  prediction  of  unobserved 
spectral  lines  cumbersome.  As  a  result,  the  emphasis 
in  this  section  will  be  placed  upon  the  relation  between 
effective  constants  defined  by 

and  the  unperturbed  spectral  constants  of  the  several 
vibrational  states,  Bv,  Dv,  qv,  q'v,  etc. 
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FIG.  1.  Spectrometer  for 
the  study  of  laser  plasmas. 


A.  Effects  of  /-type  doubling 

The  basic  theory  of  /- type  doubling  has  been  developed 
by  Amat  and  Nielsen11’12  and  extended  and  applied  by 
Maki  and  Lide. 5  The  notation  used  here  is  consistent 
with  that  of  Ref.  5. 

For  linear  triatomic  molecules  vibrational  states  are 
specified  as  i\ ,  v[,  v%,  and  states  with  /  =  0,  1,  2,  3... 
are  referred  to  as  £,  it.  A,  4>, . . .  .  For  all  states  with 
vz  5  1,  there  exists  a  multiplet  of  vibrational  states  that 
is  degenerate  in  the  harmonic  approximation.  The  con¬ 
tributions  of  /-type  resonances  to  the  lifting  of  these 
degeneracies  can  be  calculated  by  the  solution  of  a  secu¬ 
lar  determinant  of  dimension  (r>2  +  l). 

For  =  1,  the  secular  determinant  is  of  the  form 
/•S-e  IlF,, 

5^1,  E°,-< 

with  IF,,  =  q„J(J +  1).  E°  is  the  unperturbed  energy  of  the 
01*0  state.  The  rotational  contribution  to  this  energy  is 

Er--Bv\j(J+\)-P\-Dv[j(J+\)-l*\*  .  (3) 

When  centrifugal  distortion  in  the  /-type  doubling  con¬ 
stant  is  included,  Eq.  (2)  gives  for  its  solutions 

E J  =  B„\J(J  +  !)-/*]*  1) 


-D,[W+1)-/T  •  (4) 

It  should  be  recognized  that  this  equation  contains  terms 

of  the  form  B„/2  and  Dv / 4  that  are  not  functions  of  J  and 
do  not  contribute  to  the  rotational  frequencies.  In  this 
paper  these  will  be  considered  to  be  part  of  the  vibra¬ 
tional  energy,  thereby  making  Eqs.  (1)  and  (3)  com¬ 
patible.  Equation  (4)  can  be  factored  into  coefficients  of 
</(</  +  1)  and  Jl(J  +  l)2  to  give 

B:„  =  Bv±{qv  +  2Dv  (5) 

/>‘.ff  =  //„*kJ  •  (6) 

For  z>2  =  2  the  secular  determinant  is  of  the  form5 


£°a-< 

0 

»*0 

£° -e 

iv2„ 

=  0 

(7) 

0 

W'zo 

with 

Wio  =  (qv/J2)[jl(J+  1 P  -  2J(J  +  1  )]l/t  . 

Solution  of  this  equation  gives 

=  6-i{62+4<?2u[j2(J+l)2-2jy4l)]}1/2  (8) 

=  (9) 
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+  1)1- W+l))}*'*  , 

(10) 

where  6  E°&-  E^.  The  rotational  contribution  to  /}£ 
and  £‘®  is  given  by  Eq.  (3). 

For  transitions  in  the  millimeter  and  submillimeter 
microwave  region,  these  expressions  can  be  expanded 
and  compared  with  Eq.  (1)  to  give 


A 

F°c-l</2  6)lJ2(J4l)2-2J(j4l)| 

(ID 

A 

An  -  2(</2  6) 

(12) 

A 

Aff-W2  6) 

(13) 

£°a+(</2  6)UV4  1)2-2J(J  +  1)| 

(14) 

A 

Aff  4  2</J  6  -  8  A 

(15) 

A 

Air  +  (<!  1  6) 

(16) 

(17) 

A 

An  -  8  A 

(18) 

A 

(19) 

For  rz  =  3,  the  secular  determinant  is  of  the  form 


w1t 

0 

0 

<Ai 

*!-< 

0 

0 

w’n 

A°-« 

«Ai 

0 

0 

w3l 

with  tV',,  (  .3</„  2|[«fV+l)2-8J(J+  l)+12|l/*.  Solution 
of  this  equation  gives 


Eu 

£.04i64^u-H[«-H'lJ244W*,}‘'2 

(21) 

E: 

A°.-‘64U'114‘{[6-W'n|2  4  4  W21}1'2 

(22) 

Eu 

A-M«-JW,.-£{[84  H'u|*  +  4l*'‘l}1/2 

(23) 

E»c 

Ei-46-|Wu  +  i{lfl+  Wll|2  4  4W|J1/2  , 

(24) 

where  6 

Ef  -  E? .  When  centrifugal  distortion  is 

in- 

eluded, 

for  states 

A  An-</„-6v2  '6  — 2 A 

(25) 

A  =  A« 

(26) 

A  An  +  qI  62 

(27) 

For  uc  states: 

A  An  +  <lv  -  6</2  -  2  A 

(28) 

A  Aff  +  Qv  ~  i(rll 

(29) 

A  =  Att  -<f  J  /62  . 

(30) 

Only  the  /  -  0  state  of  r2  4  was  studied  and  since  there 
is  no  I  degeneracy  for  this  state,  perturbation  theory 
may  be  used  to  calculate  the  effects  of  I  doubling. 

4  12y«  [J3(J^l)3-2JV^l)2l  ,  (31) 


and 

A  At  (-6ry‘  6 


(32) 


A  -  An 


3  </l  24y,,92 


A =  An —  12>,i (/l 1  &  . 


(33) 

(34) 


The  term  which  contains  yt,  results  from  the  expansion 
of  El  -  Ea  and  accounts  for  the  difference  between  the  ro¬ 
tational  constants  of  the  04°0  and  0420  states. 


B.  Effects  due  to  Coriolis  coupling  between  11 1  0  and  04°  0 

Because  a  near  degeneracy  exists  betw'een  the  rota¬ 
tional  levels  of  the  ll'O  and  04°0  vibrational  states  at 
J  10,  it  is  not  possible  to  fit  either  the  11*0  or  04°0 
data  on  the  basis  of  the  effective  constants  of  Eq.  (1). 
The  perturbation  between  the  states  can  be  given  by4 

A2  =  62  +  4J(J+  DIF2,  .  (35) 

where  6  is  the  separation  of  the  unperturbed  states,  A 
the  separation  of  the  perturbed  states,  and  Wai  the 
Coriolis  constants.  6  is  a  function  of  J  and  can  be  ex¬ 
pressed  as 

6  AE0- AR7(«/+ 1)4  AA/V+ l)2  ,  (36) 

where  A £0  represents  vibrational  energy  difference 
between  the  two  states. 


III.  RESULTS 

A.  Observed  frequencies  and  effective  rotational 
constants 

Table  1  displays  the  forty-seven  rotational  transitions 
in  12  different  excited  vibrational  states  that  have  been 
observed.  With  the  exception  of  the  rotational  transi¬ 
tions  in  the  04°0  and  ll'O  states,  all  may  be  represented 
by  the  effective  rotational  constants  of  Eq.  (1).  These 
effective  constants  are  shown  in  Table  II.  Also  included 
are  precise  ground-state  constants  which  were  calcu¬ 
lated  from  the  results  of  a  beam  maser  experiment  of 
De  Lucia  and  Gordy.13 

B.  Unperturbed  rotational  constants 

Since  the  rotational  levels  of  the  10°0,  20°0,  00°1, 

00°2,  and  10°1  are  essentially  unperturbed,  the  effective 
constants  of  Table  II  are  also  the  unperturbed  spectral 
constants  for  these  states. 

Table  111  shows  the  unperturbed  rotational  constants 
for  all  states  for  which  corrections  between  the  effective 
and  unperturbed  constants  are  necessary.  Equations 
(5)  and  (6)  show  that  the  constants  for  Ol'O  can  be  ob¬ 
tained  by  the  averaging  of  the  effective  constants  for  the 
two  components  and  subtracting  2 Dv.  A  similar  proce¬ 
dure  was  used  for  the  01*1  state.  Equations  (U)-(13) 
were  used  to  calculate  the  constants  of  the  02°0  and  the 
results  of  Eqs.  (14)  — (19)  were  used  for  the  0220  state. 
Because  the  J  =  0-  1  and  1  -  2  transitions  do  not  exist 
for  states  with  /  2,  there  is  no  experimental  redundancy 
within  the  0220  state  and  the  listed  standard  deviations 
reflect  estimated  experimental  uncertainty.  Equations 
(25)  — (30)  were  used  to  calculate  the  constants  for  the 
03‘0  state. 
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TABLE  I.  Observed  transition  frequencies  in  vibrationally  excited  HCN  (MHz).* 


State 

J  =  0-1 

j= i 

—  2 

J -2 

—  3 

J-3 

—  4 

Observed 

Obs.  -calc.  Observed 

Obs.  -calc 

.  Observed 

Obs.-calc.  Observed 

Obs.-calc. 

0l‘*0 

... 

177  238.73" 

0.04 

265  852.  68 

-0.02 

354460.33 

0.01 

0l“0 

... 

178  136.50" 

-0.02 

267199.37 

0.01 

356  255.71 

0.00 

02°0 

89  087.  69* 

0.04 

178 170.38 

0.01 

267  243. 15 

-0.07 

356301.33 

0.04 

023“0 

... 

... 

267 109.37 

356135.46 

022c0 

.  .  . 

... 

267120.02 

356 162.77 

03**0 

... 

177  698.78 

0.02 

266540.00 

-0.04 

355  371.69 

0.02 

03  “0 

.  .  . 

179  547.32 

-0. 10 

269312.89 

0. 16 

359  067.86 

-0.07 

04°  0 

179  127.26* 

268  663.29* 

10°0 

88  027. 28“ 

0.03 

176  052.38 

-0.01 

264  073.30 

-0.01 

352  087.89 

0.01 

20°0 

87  419.40“ 

0.06 

174  836.53 

-0.04 

262  249.  60 

0.01 

349  656.  27 

-0.00 

00°1 

88  006.69“ 

0.03 

176011.26 

0.03 

264011.53 

-0.08 

352  005.76 

0.04 

00°2 

174  725.29 

262  082.87 

10°1 

87  415.29 

0.07 

262  237.31 

-0.07 

349  640.22 

0.04 

01**1 

... 

264  005. 10 

351996.75 

0l“l 

... 

176  919.08 

265  373.12 

lllc0 

... 

264019.77* 

-0.09 

352  016.88* 

0.07 

ll“0 

.  .  . 

176  913.10 

265  364.36 

“Corrected  for  quadrupole  splittings.  '‘These  lines  have  also  been  observed  by  R.  C.  Woods 

“From  Ref.  6.  (private  communication). 

*See  footnote  b,  Table  II. 


C.  The  11 '0/04°  0  system 

The  11*0  and  04°0  states  must  be  considered  together 
because  of  the  Coriolis  coupling  between  them.  This  is 
an  exceedingly  important  and  interesting  system,  it  is 
the  basis  of  the  well  known  far-infrared  laser.1*  The 
absolute  frequency  measurement  of  the  laser  emissions15 


TABLE  II.  Effective  rotational  constants  for  excited  vibrational 
states  of  HCN  (MHz). 


State 

Constant 

a 

Constant  o 

00°0 

44  315.9757 

0.0004 

0.08724 

0.00006 

0l‘*0 

44  310.381 

0.013 

0.0888 

0.0005 

oi‘“o 

44  534.851 

0.004 

0.0902 

0.0002 

02°0 

44  544.233 

0.014 

0.2054 

0.0006 

022“0 

44  519.895 

0.0926 

022*0 

44  519.563 

-0.0245 

03**0 

44  425.768* 

0.013 

0.1350* 

0.0005 

03 ‘“0 

44  887.974* 

0.050 

0.1397* 

0.0021 

04  °0 

44  785.584" 

0.051 

0.45297" 

0. 00034 

10°0 

44  013.803 

0.006 

0.0881 

0.0003 

20°0 

43  709.846 

0.014 

0.0879 

0.0006 

00°1 

44  003.505 

0.  015 

0.08716 

0.0006 

00°2 

43  681.998 

0. 0844 

10°1 

43707.783 

0.036 

0.0863 

0.0014 

01**1 

44  002.465 

0.0898 

0l‘“l 

44  230. 503 

0.0917 

ll‘*0 

44  004.714* 

0.037 

0.08873* 

0.00016 

ll‘“o 

44  229.006 

0.0912 

“Interactions  with  03s0  produce  an  H,lt  that  makes  a  significant 
contribution  at  low  J.  To  reproduce  the  calculated  frequencies 
of  Table  lj  for  03lc0,  //,„=- 1.52 x It)-5  MHz,  and  for  O3'“o, 
//,„  =  + 1. 52  xlo"5  MHz. 

"Because  of  the  Coriolis  perturbation  between  the  rotational 
levels  of  the  04°0  and  11!*0  states,  these  constants  will  not 
reproduce  the  rotational  levels  without  the  use  of  Eqs.  (35)  and 
(36)  and  the  constants  in  Table  IV. 


coupled  with  the  pure  rotational  transition  measured  in 
this  study  make  possible  the  calculation  of  a  number  of 
parameters  to  high  accuracy,  including  the  centrifugal 
distortion  parameters  for  both  states,  the  Coriolis  in¬ 
teraction  constant  H'ol ,  and  the  energy  difference  be¬ 
tween  the  two  vibrational  states  A £0 .  These  constants 
are  shown  in  Table  IV  and  the  transition  frequencies 
calculated  from  them  are  compared  with  the  observed 
frequencies  shown  in  Table  V.  It  is  interesting  to  note 
that  vibrational  energy  difference  between  the  04°0  and 
11*0  states  has  been  calculated  to  ±  0.  0001  cm"1. 

In  addition  to  the  Coriolis  perturbation,  both  the  11‘0 
and  04°0  states  are  affected  by  I- type  doubling.  Equa¬ 
tions  (32)  and  (33)  provide  the  corrections  for  the  04°0 
state,  but  values  of  qv  and  6  are  required.  The  value  of 
qv  can  be  extrapolated  from  the  values  for  the  01‘0, 

02z0,  and  03*0  states  and  Maki  gives  £°  -  £ c  to  be 
15.  225  cm"1.  The  unperturbed  values  of  both  fi04o0  and 
£>04o0  are  shown  in  Table  III.  O04o0  is  a  very  sensitive 
test  of  qv  and  6.  Its  value  is  in  good  agreement  with  its 
expected  value  and  this  serves  as  a  confirmation  of  the 
values  of  qv  and  6  used  in  Eq.  (33). 


TABLE  III.  I’nperturbed  rotational  constants  of  the  perturbed 
excited  vibrational  states  of  HCN  (MHz). 


State 

Constants 

a 

Constant  a 

ol‘o 

44  422.438 

0.007 

0.0895 

0. 0003 

02°0 

44  544 . 004 

0.014 

0.0908 

0.0006 

0220 

44  519.112 

0.050 

0.0914 

0.0015 

03‘o 

44  6  56.332 

0.026 

0.0933 

0.0010 

04  °U 

44  784.866 

0.05 

0.0937 

0.0005 

Ol'l 

44116.303 

0.05 

0.0908 

0.0010 

li'o 

44 116.680 

0.05 

0.0900 

0.0010 
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TABLE  IV.  Spectral  constants  calculated 
from  the  04°0/lllc0  system  (MHz). 


Constant 

Value 

a 

112.27 

0.06 

AEq 

78 

663.1 

4.0 

UJggOg),,, 

44 

785.58 

0.05 

^o»uo).n 

0.452  97 

0.000  34 

(BlllCg),,, 

44 

004.71 

0.04 

lB||lCg),ff 

0.088  73 

0.00016 

D.  Vibrational  effects  of  moments  of  inertia 

The  relation  between  the  effective  and  equilibrium  con¬ 
stants  can  be  expressed  as 

B„  =  B.-£  “ifr’f +5<fj)  +  X/  + 
i  1/ 

+  +  2^*^  +  Til  f  2  *  (37) 

u» 

Table  VI  shows  a  comparison  between  the  constants  cal¬ 
culated  from  extensive  infrared  data  and  the  previously 
existing  microwave  data6  and  an  analysis  of  our  micro- 
wave  data  alone.  It  is  necessary  to  include  the  e222  term 
in  our  expansion  because  of  the  significantly  larger 
amount  of  microwave  data  now  available.  Previously, 
microwave  data  for  only  two  excited  bending  modes  were 
available  and  they  essentially  determined  fl2  and  r22 .  It 
is  interesting  to  note  that  makes  a  significant  con- 


TABLE  V.  Analysis  of  the  04°0/lllc0  system  (MHz). 


Transition 

Coriolis 

v '  —v" 

Observed 

Obs.  -calc. 

Perturbs- 

J'-J" 

tion 

04°0-04°0 

2-1 

179  127.26 

0. 12 

-0.69 

04°0— 04°0 

3-2 

268 663.29 

-0.08 

-1.15 

04°0-04°0 

9-8 

804  750.  91 

0.00 

-61.35 

04°0-04°0 

10-9 

894414.2* 

0.28 

524.42 

11*0— 04°0 
10-9 

890  760.7* 

-0.  27 

-327.87 

11*0— 04°0 
11-10 

964313.4* 

0.28 

-517.22 

2—1 

176016.82 

0.08 

0.69 

ll’o-li'o 

3-2 

264  019.77 

-0.09 

1.15 

n'o-n'o 

4-3 

352  016.88 

0.07 

1.79 

li'o-n'o 

11-10 

967  965.8* 

-0.28 

335. 07 

“Measured  frequencies  of  HCN  laser  emission  lines  from  Ref. 
14. 


TABLE  VI.  Equilibrium  rotational  constant  expansion  (MHz). 


Con¬ 

stant 

This  work 

Value  a 

Reference  6 

Value  n 

Be 

44  511.62 

0.03 

44  512.36 

0.23 

<*i 

300.  028 

0.03 

299.043 

0.57 

“2 

-108.369 

0.05 

- 106. 774 

0.075 

“3 

312.  983 

0.03 

313. 163 

0.270 

Til 

-0.893 

0.009 

-1.049 

0.300 

> 12 

-3.583 

0.033 

-4.  116 

0.  075 

>13 

6.448 

0.027 

6.  625 

0.210 

722 

0.  762 

0.027 

1.532 

0.015 

>23 

6.332 

0.033 

6.434 

0.075 

>33 

-4.521 

0.018 

-4.527 

0.  008 

6  222 

0.094 

0.003 

... 

.  .  . 

>11 

-  6. 224 

0.006 

-6.056 

0.010 

tribution  to  the  B„  of  the  lower  bending  modes  as  evi¬ 
denced  by  the  differences  between  the  Ba  ,  a2 ,  and  Y& 
of  the  two  analyses. 

Although  our  constants  result  from  a  fit  of  12  constants 
to  only  13  microwave  data  points,  the  constants  which 
result  accurately  predict,  with  the  several  exceptions, 
the  infrared  A B  which  were  used  in  the  previous  analysis 
and  the  weighted  addition  of  this  infrared  data  to  our 
analysis  makes  no  significant  change  in  our  calculated 
constants. 

The  AB’s  which  are  not  well  predicted  are  interesting. 
The  value  of  A B  for  0220  -  02°0  is  -  80.  72  ±  0.  30 x  1  O'5 
cm*1.5  The  microwave  value  for  this  A B  is  -  83.02 
±0. 17xl0*5  cm*1,  an  analysis  of  all  microwave  data  pre¬ 
dicts  -  83. 05±  0. 03x  10*5  cm*1,  and  an  analysis  of  all 
microwave  data  except  B0220  predicts  -83.09xl0*s  cm*1. 
Since  the  calculation  of  the  previous  value  requires  a 
complicated  and  indirect  calculation,  assumptions  about 
other  spectroscopic  constants,  and  since  a  more  recent 
version  of  this  calculation  now  gives  -  82.  23 ±  0.  24 x  10*5 
cm*1, 16  it  would  seem  reasonable  to  conclude  that  this 
does  not  represent  any  real  inconsistency  with  our  anal¬ 
ysis. 

Our  predictions  of  the  A B  for  the  12°0-01'0  band  dif¬ 
fers  by  about  three  times  the  stated  uncertainty  of 
±0.  7x  10"5  cm*1.  But  the  difference  between  this  A B 
and  the  AB  for  the  1220  -  01 ‘0  is  -  80.  80 ±  1. 1  x  10*5  cm*1. 
This  value  is  in  excellent  agreement  with  the  0220  -  02°0 
AB  of  Ref.  6,  but  differs  from  both  the  newer  ABo^ozOo 
of  Ref.  16  and  our  microwave  data. 

The  value  of  AB  for  05*0  -0l'0  of  1590±  15  cm*1  dif¬ 
fers  from  our  predicted  value  by  20  cm*1.  It  is  interest¬ 
ing  to  note  that  the  12°0,  1220  and  05'0  are  involved  in  a 
Coriolis  interaction  and  were  analyzed  together. 

Recently,  Maki16  has  made  some  highly  precise  mea¬ 
surements  on  hot  baude  of  HCN  and  has  results  which 
question  the  completeness  of  Eq.  (37).  For  example,  he 
finds  AB(12°0  -  02°0)  =  —  999. 26(31)x  10*5  cm  and 
AB(1220  -  02z0)  =  -  999.  95(29)x  10"5  cm.  In  the  context 
of  Eq.  (37)  these  two  values  should  be  the  same.  At 
face  value  these  results  require  terms  to  be  added  to 
Eq.  (37)  which  mix  v(  and  l.  However,  it  is  also  pos¬ 
sible  that  since  the  analyses  used  to  produce  these  re- 
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suits  require  the  fixing  of  energy  levels  and  spectroscopic 
constants  derived  from  microwave  and  other  infrared 
experiments  and  since  the  uncertainties  quoted  do  not 
include  systematic  errors  in  the  experiments,  errors 
in  the  fixed  energy  levels  and  constants,  or  model  er¬ 
rors,  the  actual  uncertainties  are  somewhat  larger  and 
no  inconsistency  exists.  It  should  be  pointed  out  that 
the  A B  predicted  from  the  microwave  numbers  and  Eq. 
(37)  generally  fall  within  the  range  of  values  for  the 
infrared  multiplet.  In  addition  to  providing  improved 
spectral  parameters,  this  calculation  is  necessary  for 
the  assignment  of  the  rotational  transitions  in  several  of 
the  vibrational  states  studied.  This  is  because  of  the 
near  degeneracy  caused  by  the  similar  values  of  Bell  for 
ll'^O  and  01“l  and  for  001,  llle0,  and  01lcl. 

IV.  CONCLUSIONS 

It  has  been  demonstrated  that  millimeter  and  sub¬ 
millimeter  microwave  spectroscopy  is  an  effective  and 
sensitive  probe  of  active  laser  plasmas.  Since  these 
plasmas  contain  large,  nonequilibrium  populations  in 
excited  vibrational  states,  these  studies  also  yield  sub¬ 
stantial  amounts  of  new,  high  precision  microwave  data. 
These  results  make  possible  the  calculation  of  the  ro¬ 
tational  spectrum  of  HCN  over  a  wide  range  of  vibra¬ 
tional  and  rotational  states.  For  most  states  these  cal¬ 
culations  are  especially  convenient  because  of  the  re¬ 
lations  provided  between  the  effective  and  unperturbed 
constants. 

It  should  also  be  noted  that  this  is  a  technique  that  can 
be  used  for  the  study  of  the  otherwise  inaccessible  ex¬ 
cited  vibrational  states  of  many  molecules.  Because  the 
kinetic  and  rotational  temperatures  remain  approximate¬ 
ly  ambient,  this  technique  is  substantially  more  sensi¬ 
tive  than  the  hot  cell  method,  especially  for  highly  ex¬ 
cited  states. 
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"Forbidden"  Millimeter-Wave  Transitions  in  Arsine  1 

Jiavid  A.  Helms  and  Walter  Gordy 

Department  oj  Physics,  Duke  University,  Durham,  North  Carolina  27706 

A  millimeter-wave  spectrometer  having  a  sensitivity  of  4  X  10-10  cnr1  in  the  2-mm  region 
has  been  used  for  observation  of  the  “forbidden”  transitions  J  — >  J ,  K  =  ±4— >  ±1  and 
J  —r  J ,  K  =  ±5  — *  ±2  in  AsH3.  A  comprehensive  computer  analysis  was  made  of  the  fre¬ 
quencies  measured  in  this  work  together  with  available  microwave  frequencies  of  other  transi¬ 
tions.  This  analysis  provides  accurate  values  of  the  rotational  constants,  nuclear  quadrupole 
couplings,  and  effective  structural  parameters  of  the  molecule.  The  spectral  constants  B0  and 
Co  (in  MHz)  are  112  470.597  and  104  884.665,  respectively. 

I.  INTRODUCTION 

In  a  previous  paper  (/)  we  reported  observations  of  forbidden  AK  =  ±3  microwave 
transitions  in  PH3  and  PDj.  Our  measurements  of  A'  =  3  — *  0  transitions  in  PH3  and 
K  =  3— >0  and  K  =  ±4— >  ±1  transitions  in  PI)3  confirmed  the  assignments  of  the 
lower-frequency  A  =  ±2  — >  Tl  transitions  of  PH3  and  of  PD3  reported  by  Chu  and 
Oka  ( 2 ),  and  improved  significantly  the  molecular  constants,  especially  Co.  A  similar 
investigation  is  here  described  for  arsine. 

Transitions  in  AsH3  obeying  the  AJ  =  0,  AK  =  ±3  selection  rules  were  first  observed 
by  Chu  and  Oka  (J).  They  reported  eight  transitions  obeying  the  selection  rules  AJ  =  0, 
K  =  ±2— >  T 1,  AF  —  0  in  the  frequency  range  of  21.16  to  21.84  GHz.  The  large  As 
quadrupole  coupling  splits  the  K  =  ±2  — >  =Fl  line  into  two  doublets,  each  component 
of  which  is  an  unresolved  doublet.  Chu  and  Oka  combined  their  value  of  Co  with  the 
known  value  of  Ii0  (3)  and  with  other  information  to  obtain  the  zero-point  effective, 
the  zero-point  average,  and  the  equilibrium  structures.  Recently,  Olson  ct  al.  {4) 
reported  their  observation  of  ‘‘perturbation-allowed”  transitions  of  the  type  A\K  —  l\ 
=  ±3  in  the  iq  and  iq  infrared  bands  of  AsH3.  They  combined  their  A|A'  —  / 1  =  ±3 
infrared  measurements  with  the  K  =  =fc2  — >  =F  1  microwave  measurements  of  Chu 
and  Oka  and  with  the  normally  allowed  AK  =  0  microwave  measurements  of  Helminger 
et  al.  (J)  to  obtain  both  a  set  of  ground  state  molecular  constants  and  an  associated 
zero-point-effective  structure.  Olson  et  al.  then  used  the  constants  for  the  vt  and  iq 
bands  given  by  Sarka,  ct  al.  (5)  with  their  constants  for  the  >q  and  >q  bands  to  determine 
the  equilibrium  structure.  Our  detection  of  the  A  =  ±4  — ►  ±  1  and  the  K  =  ±5  — »  ±2 
forbidden  transitions  was  assisted  by  a  preliminary  prediction  of  their  frequencies  from 
the  results  of  Helminger  et  al.,  of  Olson  et  al.,  and  of  Chu  and  Oka. 

1  This  work  was  supported  by  the  U.  S.  Army  Research  Office,  Grant  No.  DAAG29-77-G-0007. 
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The  need  for  observation  of  the  higher  K  =  ±4—*  ±1  and  K  =  ±5  — >  ±2  milli¬ 
meter-wave  transitions  described  in  this  paper  was  indicated  by  Chu  and  Oka  (2),  who 
reported  that  their  spectrometer  was  not  sufficiently  sensitive  at  the  higher  frequencies 
to  detect  them.  Measurements  of  these  transitions  confirm  the  assignments  of  the 
frequencies  measured  by  Chu  and  Oka  and  provide  increased  accuracy  in  the  values  of 
the  spectral  constants,  particularly  those  influenced  by  the  distortion  caused  by  rotation 
about  the  symmetry  axis.  For  observation  of  very  weak  transitions  such  as  these  in  the 
shorter  millimeter-wave  region,  we  designed  an  exceptionally  sensitive  millimeter -wave 
spectrometer,  which  is  described  in  the  earlier  paper  (/). 

The  0-branch  A  A'  =  ±3  transitions  represent  only  one  class  of  forbidden  transitions 
in  symmetric  top  molecules  included  in  the  general  theory  of  Watson  (6,  7)  and  others 
(A,  9).  Oka  (10)  has  recently  written  an  excellent  review  of  the  experimental  and 
theoretical  developments  of  this  relatively  new  type  of  spectra.  In  the  ground  vibrational 
state,  the  off-axis  dipole  moment  which  couples  the  radiation  field  to  the  rotation  is 
induced  bv  centrifugal  distortion.  The  selection  rules  are  determined  by  the  molecular 
symmetry.  In  As  Hi,  the  rotation  about  the  A-axis  of  inertia  induces  a  small  dipole 
moment  (~ 3.938  X  10~5  D)  perpendicular  to  the  symmetry  axis  which  leads  to  a  slight 
admixture  of  the  A  levels  and  a  breakdown  of  the  familiar  A  A'  =  0  selection  rule.  The 
threefold  symmetry  about  the  c-axis  limits  this  admixture  to  K  levels  differing  by  three 
units.  It  is  of  interest  that  centrifugal  distortion  due  to  rotation  about  the  symmetry 
axis  can  likewise  generate  a  small  dipole  component  along  the  symmetry  axis.  This 
component  is  of  no  particular  consequence  in  molecules  such  as  PH3,  PD3,  and  AsH3 
which  already  have  permanent  dipole  moments  along  the  symmetry  axis,  but  it  can  give 
rise  to  forbidden  A/  =  ±1  transitions  in  molecules  such  as  CH,  which  has  no  permanent 
dipole  moment.  In  fact,  A7  =  1  rotational  transitions  in  CH,,  SiH,,  and  GeH,  have  been 
observed  by  Rosenberg  el  al.  (11  14)  in  the  far-infrared  region.  Kreiner  et  al.  (15)  have 
used  Laser-Stark  spectroscopy  to  measure  the  distortion-induced  dipole  moment  of 
GeH,.  Forbidden  transitions  of  CH,  in  the  ground  vibrational  state  have  also  been 
detected  in  the  microwave  region  by  Holt  el  al.  (16).  At  a  recent  symposium,  Kagann 
et  al.  (17)  reported  their  observations  of  over  140  0-type  A  A'  =  ±3  microwave  transi¬ 
tions  in  OPF3. 


II.  EXPERIMENTAL  DETAILS 

The  millimeter-wave  spectrometer  and  the  experimental  procedure  used  in  these 
measurements  are  described  in  the  earlier  paper  (/).  The  high  sensitivity  of  this  spec¬ 
trometer  (  =  4  X  10-1"  cm  1  at  144  GHz  for  a  time  constant  of  1  sec)  made  possible  the 
measurement  of  these  weak  AA’  =  ±3  transitions  with  an  estimated  accuracy  of  three 
parts  in  107. 


III.  THEORETICAL  BACKGROUND 

As  a  result  of  centrifugal  distortion,  the  rotational  wavefunction  <Pj.k  of  a  molecule 
having  C3s  symmetry  mixes  slightly  with  the  wavefunction  4'j.k±i  in  the  ground  vibra¬ 
tional  state.  The  term  in  the  Hamiltonian  responsible  for  the  interaction  may  be  ex- 
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pressed  (IS)  as 

if  =  — —  [(7+3  +  +  j.u+*  +  y_*)]t  (i) 

where  r„„  is  a  centrifugal  distortion  constant  and  ./.,  7  +  ,  and  7_  are  rotational 
operators. 

In  the  absence  of  the  perturbing  II'  term,  the  rotational  levels  are  all  doubly  degenerate 
for  each  value  of  A,  (A'  +  0).  The  energy  levels  corresponding  to  K  ^  3  are  raised  or 
lowered,  but  they  remain  degenerate  with  symmetry  E.  The  A'  =  3  level,  however,  is 
split  into  .1 1  and  ,1>  components.  The  designations  .1 1.  .4»,  and  E  indicate  the  symmetry 
of  the  rovibrational  wavefunction  (excluding  spin).  In  the  case  of  As  Hi,  each  of  the 
energy  levels  is  split  further  into  (21  +  l)  components  by  the  quadrupole  interaction  of 
the  As(I  =  |)  nucleus.  The  symmetry  selection  rules  for  (7-type  A  A  =  ±5  transitions 
(19)  are  +  «-»—,  along  with  .  1 1 «-» .1  ••,  or  E  <-»  E. 

The  formulas  giving  the  energy  levels  may  be  expressed  (20,  21)  for  A  X  3  as: 

E(F,  J,  A) 

— - —G(J,  K)  -f- iQC'.v +  (Ck — Cx)  A'2  J  (7-1-1)]] 

h 


A'2 

i - 1 

7(7+1) 


1(7,  ],/•)  + 


Ah 


X  [A(A+ 1 )-/(/+ 1  )-7  (7+ 1 )]- .-(A, 

|  (2A'-3)![7(7+l)-A'(A--l)][7(7+l)-(A'-l)(A--2)][7(7+l)-(A'-2)(A'-3)] 

xj - - - 

f  C(7,  A)-G(7.  A -3) 

(2A'+3)21]7(7+1)-A'(A-+1)J7(7+1)-(A'  +  1)(A'+2)X7(7+1)-(A-+2)(A'+3)]| 


;  (2) 


C'J.  A")  — G(7,  A' +3) 


and  for  A  =  3  as 

/•-(/•',  J,  3) 


=  (-■(./,  3)  +  J[C'.v  +  (CK  -  C\)K-  7(7  +  1)] 


X  [F(A  +  1)  -  /(/  +  1)  -  ./(./  +  1)]  -  ,\)q  ^3 

/ IOt x,j  f 

+(ir) 

+  | /to  -  is  / PV.  -  <-V,  3)]J 


A- 


7(7  + 1) 

81[7(7  +  1)  -  12J7(7  +  1)  -  20j/(7  +  1)  -  50] 
G(J,  3)  -  G(J,  b) 


l)  IV.  I  F) 


X  7(7  +  1)[7(7  +  1)  -  2X7(7  +  1)  -  b],  (3) 
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TABLE  I — Continued 


Transition 

Observe  1 
Frequency 
v  (MHz) 

v  t  c- 
*  Calc 

a 

max 
■  • 

00" 10  cm’1) 

J  —  J  K  * 

4  -  t  1  K  -  F  (C) 

(continued) 

11 

23/2 

107  564. 50 

0.02 

22. 1 

11 

21/2 

107  606.  26 

-0. 05 

22.  1 

11 

26/2 

■ 

-0. 05 

22.  1 

11 

13/2 

107  5.  80.  14 

-0. 02 

22.  1 

12 

25/2 

10C  354. 89 

-0.02 

21. & 

12 

23/2 

10C  356. 38 

0.04 

21.5 

12 

27/2 

IOC  306. 72 

0.0l> 

21.6 

12 

21/2 

ioc  see,  1? 

0.  11 

21. & 

13 

27/2 

105  072.41 

-0.06 

19.4 

13 

25/2 

105  073.  65 

0.06 

19.4 

13 

29/2 

10:  082.50 

-0.02 

19.4 

13 

23.2 

105  083. 61 

-0. 03 

•  19.4 

14 

29/2 

103  725. 04 

-o.oc 

10.3 

14 

27/2 

103  725.  85 

-0.  10 

10.  3 

14 

31/2 

103  723.  tiO 

16.  3 

14 

25/2 

103  734.  <V 

-0.06 

10.  3 

li. 

31/2 

102  322. 06 

12.9 

lb 

29/2 

102  322.95 

12.9 

15 

33/2 

102  329. 74 

0.05. 

12. 9 

J  —  J  K  -  * 

F  -  F  (c) 

J  .»  8 

F  «  17/2 

155  785.04 

20.  3 

8 

1d/2 

15!’  791.36 

0.02 

20.  3 

8 

19/2 

155  820.  52 

-0. 02 

20.  3 

8 

13/2 

156  626.92 

20.  3 

11 

19/2 

154  443.  53 

0.  10 

28.  4 

9 

17/2 

28.  4 

l* 

21/2 

164  4  .  ■ 

-0. 03 

20.  4 

15/2 

470.  *i0 

-0. 04 

28.  4 

10 

21/2 

152  972.  !.$» 

-0. 0” 

34.7 

10 

19/2 

162  975.92 

-0.  10 

34.7 

10 

23/2 

1d2  995.93 

0.01 

34. 7 

iO 

17/2 

152  999.27 

-0. 03 

34.7 

( continued ) 


The  V(J,  1,  F)  which  appears  in  the  contribution  from  the  quadrupole  inteiaction 
may  be  expressed  as: 

fC(C  +  1)  -  /(/  +  1 )/(/  +  1) 

YU,  I,  F)  - - ,  (6) 

2/(2/  -  1)(2/  -  1)(27  +  .1) 

where 

C  =  /•’(/•'  +  1)  -  /(/  +  1)  -  J(J  +  1),  and  where  F  =  J  +  I.  (7) 

(/  =  §  in  the  case  of  AsHj  considered  here.)  The  K  =  3  levels  contain  the  term 
zkhoJ  (J  +  1 )[/(./  +  1)  —  2J[i(7  +  1)  —  6],  first  derived  by  Nielsen  and  Dennison 
(22)  in  their  explanation  of  the  anomalous  effects  in  K  =  3  lines  observed  by  early 
microwave  spectroscopists.  Unfortunately,  we  were  unable  to  evaluate  the  splitting 
constant  h„  for  As  If  because  the  frequencies  of  the  (7-type  K  =  3— *0  series  occur 
below  the  operating  range  of  our  spectrometer. 

The  line  strengths  of  the  forbidden  AK  =  ±3  transitions  (6,  23)  are  given  by 

S  =  i  («.")„„  2(J  =F  K)(J  TA  -  1)(./  TA  -  2)(J  ±K  +  1)0/  ±K  +  2) 

X  (J  ±  K  +  3)(2J  +  \)/J(J  +  l),  (8) 
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TABLK  I  Continual 


where 

(0/0., ff  =  0/‘  +  h*rxxxitx./ [2 l'(H  -  C)]  (») 

is  the  effective  dipole  moment  inducing  the  transition.  The  term, 


is  an  induced  dipole  component  which  results  from  admixture  of  the  ground  vibrational 
state  with  excited  vibrational  modes.  The  last  term  of  Kq.  (9)  results  from  centrifugally 
induced  mixing  of  rotational  levels  of  the  ground  vibrational  state.  In  this  expression, 
r,,„  is  the  centrifugal  distortion  constant  in  //'  of  Eq.  (1),  n,  is  the  permanent  dipole 
moment,  and  C  and  />  are  the  s|)ectral  constants  with  respect  to  the  symmetry  axis 
and  the  perpendicular  axis,  respectively. 

The  contribution  of  the  "vibrational  intensity  borrowing”  [first  term  in  Eq.  (9)J 
is  small  in  comparison  with  that  of  the  “rotational  intensity  borrowing"  [second  term 
in  Eq.  (9)3  (2).  For  l’H;i,  the  relative  contribution  of  these  two  terms  to  the  forbidden 
line  strengths,  Kq.  f8),  is  (1.6  8.3)2  =  0.0.17,  and  it  is  estimated  that  the  relative 
contribution  is  even  smaller  in  AsHa.  Thus,  for  Asll-t  in  the  ground  vibrational  state 
the  contribution  of  the  vibrational  mixing  to  the  intensities  of  the  A  A  =  ±3  transitions 
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I'  li;.  I  Single  scan  of  ()-ty])c  K  ±5  — *  ±7  multiplrt  for  ./  14. 


is  negligible  in  comparison  with  that  from  rotational  mixing.  YVc  used  the  value  of 
r,,„  from  the  analysis  of  the  transition  data  to  calculate  the  value  of  (Ox'T).-tt  used  in 
the  calculation  of  the  intensities.  The  peak  absorption  coefficients,  listed  in  Table 
I  were  calculated  from  the  formula, 

1.147  X  10' 1,1  /  hv 0\ 

a„,„(cm  ')  -  eu'  (  1 - ) 

AgT*  \  2kT/ 

X  exp[  —  /'.(/•', ./,  K)  kT](B'-C)^KK„S[(l\  ./,  K  ±  41  —  (/•',  J,  A.)],  (11) 
where  .S'  is  given  by  Kq.  (8),  T  =  4(H)  K,  g„  is  the  reduced  nucieat  statistical  weight 


TA1IL.E  11 
Molecular  Constants 


Milo'  it'll* 

V  iltii' 

t’taihl'iril 
l  wlat  i'ti 

R. 

11"  la'.  !t4»" 

0. 030 

mu.- 

r 

104  884.  *•<’!' 

0. 043 

MM 

t>V 

0. 0030 

MM.' 

y>K 

-n.  7  h*4 

0.  0004 

MM.' 

a.  n:v 

0. 002 1 

MU'. 

t. ; .  / 

hi  0.  t'VO 

0. 003 
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factor,  is  the  statistical  weight  factor  for  each  A'  level,  and  Sr  =  20  MHz  is  the  lint- 
breadth  at  a  pressure  of  1  Torr.  This  value  for  the  Sr  is  the  same  as  that  used  by  Chu 
and  Oka  (J)  for  AsH;i. 


IV.  RESULTS 

Twenty-seven  A  =  ±4— >  ±1  and  thirty-six  K  =  ±5— »  ±2  transitions  were  mea¬ 
sured  in  this  work.  The  observed  frequencies  are  given  in  Table  I.  Figure  1  shows  a 
chart  recording  of  a  single  scan  of  the  Q-type  K  =  ±5  — *  ±2  multiplet  for  J  =  14.  The 
quadrupolt-  doublet  is  well  resolved. 

A  combined  computer  analysis  was  made  of  the  available  frequencies  for  different 
types  of  transitions  measured  by  us  and  by  others  as  listed  in  Table  I.  The  analysis 
consisted  of  a  least-squares  fitting  of  thirteen  molecular  constants  to  the  observed 
frequencies  of  AsH.,.  The  program  calculates  the  difference  in  the  energy  levels  by  use 
of  Eqs.  (2)  (7)  with  assumed  values  of  the  molecular  constants.  Transition  frequencies 
are  calculated  and  compared  to  the  measured  transition  frequencies  to  generate  new 
estimates  of  the  molecular  constants.  This  iterative  procedure  is  repeated  by  the 
computer  until  convergence  is  attained.  The  degree  of  fitting  of  the  complete  frequency 
set  is  shown  by  column  3  of  Table  1.  The  molecular  constants  for  AsH3  obtained  from 
the  combined  analysis  are  listed  in  Table  II. 

The  two  independent  zero-point-effective  structural  parameters  for  AsH3,  r„  and  ao, 
were  calculated  from  the  H0  and  the  resulting  values  are  listed  at  the  bottom  of 
Table  II.  The  standard  deviations  recorded  for  these  values  reflect  only  those  un¬ 
certainties  in  the  measurement  of  H„  and  Co  and  do  not  indicate  the  possible  errors  in  the 
absolute  values  which  depend  on  Planck’s  constant. 

Received:  August  19,  1977 
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Centrifugal  Distortion  Analysis  of  the  Ground  Vibrational  States 
of  H,,:0  and  fi.'O1 


Paul  Hklmingek2 

Department  of  Physics,  L'niversily  of  South  Alabama,  Mobile,  Alabama  3MSS 

AND 

Frank  ('.  I)k  Lucia 

Department  of  Physics,  Duke  l'niversily,  Durham.  Sorth  Carolina  277<>6 

Centrifugal  distortion  analyses  of  combined  infrared-microwave  data  sets  for  1 F/  70  and 
H  ,'H)  have  been  performed.  These  analyses  make  possible  critical  evaluation  of  data  points 
and  results  in  substantially  improved  energy  levels.  The  I  raley  Kao  interpolation  rule  for 
the  prediction  of  Hd'O  energy  levels  is  examined  and  extended. 


INTRODUCTION 


Several  years  ago  we  published  centrifugal  distortion  analyses  for  a  number  of 
isotopic  species  of  water.  F.xcept  for  the  H2";0  (I),  Hd70  (2),  and  FI/'O  (J)  species 
these  were  based  entirely  upon  microwave  data.  Because  for  these  species  it  was 
not  possible  to  measure  enough  microwave  lines,  it  was  necessary  to  include  infrared 
data  in  our  calculations.  Although  at  the  time  of  our  analyses  extensive  Hdf>0  in¬ 
frared  data  was  available,  only  a  modest  amount  of  Ild'O  data  (4.5)  and  virtually 
no  H/'O  data  had  been  published.  In  this  paper  we  report  new  analyses  of  H>170 
and  Hd'O  that  include  substantial  amounts  of  new  infrared  data  in  addition  to  our 
earlier  microwave  data. 

Because  a  centrifugal  distortion  analysis  has  substantially  fewer  degrees  of  freedom 
than  a  purely  algebraic  combination  difference  analysis  and  because  the  former  are 
more  sensitive  to  "bad”  data  points  which  can  be  eliminated,®  more  accurate  and 
reliable  energy  levels  can  be  obtained.  In  addition,  the  independent  analyses  of  the  three 
isotopic  species  make  possible  a  critical  evaluation  and  refinement  of  the  Fraley-Rao 
interpolation  rule  for  the  prediction  of  H2170  energy  levels. 

1  This  work  is  supporter!  by  the  U  S  Army  Research  Office,  Grant  PAAG29-77-G-0007. 

2  Travel  to  Duke  University  supported  by  Southern  Regional  Kdueational  Hoard. 

;  For  an  enlightening  discussion  of  a  mathematician’s  view  on  the  rejection  of  data,  see  F.  J. 
An-i i im be  and  J.  \V  Tukey.  Ter  hnometrics  5,  141  (1963). 
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ll2‘!0  ANALYSIS 

We  have  combined  4<>  of  the  far-infrared  lines  reported  by  Winther  (0),  31  com¬ 
bination  dillcrences  calculated  from  the  data  of  loth  t'l  ill.  (S),  20  combination  dif¬ 
ferences  calculated  front  the  data  of  Caniy-l’evret  <’/  ill.  (')),  and  our  seven  microwave 
lines  (_’)  in  a  weighted  Watson  analysis.4  The  assigned  weights  were  inversely  pro¬ 
portional  to  the  square  of  the  expected  uncertainty  and  ranged  front  It)7  for  the 
microwave  lines  to  one  for  W inther  s  weakest  lines.  In  order  to  carefully  screen  the 
data,  each  of  the  104  data  points  was  removed  (in  groups)  front  the  analysis  and 
predicted  on  the  basis  of  the  remaining  data.  With  the  exception  of  a  few  lines  at 
high  J ,,  the  uncertainty  in  each  of  the  predictions  was  <  ±0.01  cm  1  and  most  were 
much  better,  lor  all  the  lines  contained  in  the  data  set.  the  agreement  between  the 
calculated  uncertainly  in  the  prediction  and  the  expected  experimental  uncertainty 
was  satisfactory  except  for  the  data  points  listed  in  Table  1.  While  only  one  of  these 
data  points  is  dramatically  bad  (the  rest  are  oil  by  about  0.05  cm'1),  thev  are  all 
substantially  worse  than  the  data  retained  in  the  analysis.  The  spectral  constants 
which  result  front  our  analysis  are  shown  in  Table  11  and  the  energy  levels  calculated 
from  them  in  Table  111.  Tour  places  are  retained  in  the  energy  levels  because  it  has 
been  found  that  the  analyses  successfully  predict  unanalyzed  microwave  lines  to 
this  accuracy. 


Hs‘*0  AMI  ll;‘»o  ANALYSTS 

The  lfj"0  energy  levels  shown  in  Table  111  result  from  an  analysis  very  similar 
to  the  1I,170  analysis  discussed  above  except  that  each  of  the  data  subsets  is  larger 
and  a  total  of  240  data  points  were  analyzed.  At  the  highest  it  was  not  possible 
to  establish  with  certainty  the  reliability  of  the  infrared  data  points  because  of  the 
scattered  nature  of  the  data.  However,  at  lower  ./,  the  deviations  of  several  points 
from  the  lit  were  inconsistent  with  their  assigned  weights  and  they  were  eliminated 
from  the  lit.  These  are  shown  in  Table  IV. 

I  he  IT  "O  energy  levels  listed  in  Table  111  are  those  which  result  from  the  original 
analysis  (/>  of  15  microwave  lines  and  the  infrared  data  available  to  us  at  that  time. 


TAB1.K  I 


Data  Points  Removed  from  lbl7()  Analysis 


- 

f.n-t) 

Ref. 

•  - 

0.049 

, 

ft41  *  ”43 

-0  *'6? 

t  Winmer  iRef.  «•,)  >f  on>* 


of  me  transition*;  f  Hof.  0  which  leads  to  the  <’j(.  -  8^ 

*  in  i’i.  a  liff^rerv-  t;oenw;  *o  be  . . iivt  is  ,n- 

•  1  .  !•’  1  in  .  'ir  unlyu.i;;. 

‘  l  «>r  a  ilctaili'il  ili>»  ussinn  of  our  analysis  lit  liniiiuo,  sec  Ri*f.  (/I 
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TABLE  II 


Rotation  and  Distortion  Constants  of  the  Ground  Vibrational  State  of  Water  (MHz) 


Constant 

h216° 

h2170 

«2  0 

Value  3 

a 

Value3 

0 

Value  a 

a 

a 

835  840.  288 

0.  50 

830  282.791 

825  366.  844 

1.0 

y 

435  351.  717 

0.45 

435  357.  289 

2.( 

435  356.  68: 

0.8 

u 

278  138.  700 

0.  45 

277  505. 735 

2.6 

276  948. 996 

0.8 

aj 

37.59422 

0.02 

37.59414 

0.09 

37.  54579 

0.017 

ajk 

-172.9128 

0.  17 

-171.8532 

0.3 

-171.  1683 

0.  13 

\ 

973. 29052 

0.  10 

960. 46634 

0.5 

949.  88229 

0.2 

6J 

15.21040 

0.010 

15. 34663 

0. 04 

15.23181 

0.011 

6k 

41.0502 

0. 15 

41.4538 

1.2 

3b. 8850 

0.2 

HfdO2) 

1.  56556 

0.020 

1.82110 

0.  12 

1.  53906 

0.019 

-4.2081 

0.6 

-3.  1738 

0.4 

-5,09 

0.  10 

-6.  42468 

2.6 

-5.21896 

0.  12 

Man 

3. 733028 

0. 008 

3.745123 

0. 06 

3.  624767 

0.017 

hj-UO3) 

7. 79579 

0.09 

10. 18803 

0.6 

8. 04641 

0.  14 

hJK 

-2.  5165 

0. 11 

3. 7859 

2.4 

-2.4187 

0.25 

V<10  > 

1.0971 

0.03 

1.  1495 

0.  13 

1.0486 

0. 04 

Ljk.(103) 

-3. 0647 

0. 46 

-4.  6449 

0.6 

W(1°2> 

1.02952 

0.  12 

1.46700 

0.  16 

(102) 

-2.  340138 

0.08 

-1.517327 

0.  17 

-2.617359 

0.  11 

•s-ooe) 

-2.59 

0.8 

'■K' 

-4.257 

2.6 

»K-  (to3) 

- 1.  3546 

0. 10 

-0.  326 

0.2 

-1.32006 

0.  13 

Py-  (106) 

5.  19841 

0.07 

4.71660 

1.5 

5.72244 

0.4 

pk-  ocn 

3. 7603 

1.0 

7 The  large  •  -nicer  of  iigits  retained  In  the  higher  order  constants  is  required  in 
order  to  reproduce  the  energy  levels  to  experimental  uncertainty. 


More  recent  energy  levels  (10),  based  upon  the  high  resolution  experimental  data 
of  Guelachvili  (and  apparently  upon  our  microwave  data,  too),  result  in  essentially 
the  same  energy  levels  over  the  J ,  range  of  the  original  analysis. 


I  RALEY  RAO  INTERPOLATION  RULE 
The  Fraley- Kao  (4)  interpolation  rule  can  he  stated  as 

v\i  —  i'le 

- -  k  (1) 

>*18  —  •'IS 

I  able  V'  shows  the  result  of  the  application  of  this  rule  to  the  observed  microwave 
transitions  of  Hjl70.  It  is  obvious  that  large  deviations  exist.  A  similar  effect  has  been 
observed  for  microwave  transitions  of  I  TO  (11).  However,  critical  inspection  of  the 
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TABLE  III 

Energy  Levels  of  the  Ground  Vibrational  States  of  Water  (cm  l> 


h,16o 

h„17o 

1  8 

«2  0 

23. 794 o 

23. 7735 

23.  :  . 

37. 1372 

36. 9311 

3€.  48" 

l\l 

42. 371 0 

42. 1670 

42. 0235 

70.0911 

70. 0047 

69.9276 

7y.  4965 

7y.2273 

’3.9  367 

95.  1762 

91.9706 

0 

0 

.-,11 

13  *.  .*019 

134. 1453 

133.  •,  759 

2‘ 

*20 

13S.  16-52 

138. 4312 

134.7833 

3rt„ 

130. 7622 

136. 8378 

. 

142.2788 

141. 902 1 

141..  682 

3!^ 

173.366:' 

173.1102 

1 6832 

206.3019 

208. 4820 

204.  •;  '*1 

212. 1568 

*  11.  ,3:  9 

210.799:- 

. 

25.*'.  2199 

2b0.  9615 

282.0950 

330 

289.  ilyl 

283.  7or 

•. 

4o. 

222. 0536 

221.6212 

221.2344 

224.  8386 

224.1  ■: 

223.  6288 

jl3 

276.  1981 

275. 1305 

•  s  : 

300. 3028 

299.  1369 

7  .< 3. 6206 

4oo 

313.7801 

31:  . 0786 

314  •  . . 

-it" 

382.  5177 

3oG.  8068 

3  79. 292  5 

383.  6432 

382. 1772 

380.  03: 

•,38.  1083 

185. 208*.'. 

,82.6445 

,68.  1348 

,3: . '  381 

482.6736 

to 

^05 

32:’.34d9 

32-;  *613 

3?..;  0472 

326.  '3261 

325. 6803 

31  .2161 

396.  67  67 

39a.  3612 

41*'.  20y4 

41:  .  1269 

51  ;.  1 689 

446. 5114 

,4;'.  7931 

■ 

•  02.966.' 

:  '02 .  !  807. 

500.  59 73 

532 

508.  5127 

907.  1/64 

509.  ’293 

8 

010. 1 loc 

60  7.  1629 

• 

610. 342s 

60".  101  t 

»>04.  9  la 

74". 0744 

737.6190 

33.  <  81 

c,r  1 

742. 0777 

737. 6225 

733. 6854 

nu 

eG1 

446. 6976 

443.  71  A3 

444.  346 " 

447. 2531 

446. 2446 

■  15.  U 

542. 9076 

541. 9946 

:  41.1810 

552.9121 

1.6058 

.4515 

6r>\ 

602.  7  M3 

601.9593 

601.2382 

tf6** 

648. 9786 

647. 0694 

645  3837 

061.:-  489 

*>5o.  9894 

658.6111 

6 

786. 7262 

753.  7094 

751.0384 

757.  7616 

754. 8182 

752.  1 699 

888.  6029 

864.  Odor 

680.  03O'; 

888.  6368 

664.  1228 

880. 1189 

1045.  0623 

1038.  7644 

1033. 19*’ 

660 

1048. 0627 

1033.  7648 

1033. 1991 

7r_ 

$66. 2-Mtf 

r>84.  94^3 

m83.  7782 

586.4800 

885.  1610 

583. 9869 

7IT 

704.2166 

702.  6611 

701 . 6952 

709. 6092 

706.  0092 

706.  5992 

72r> 

782.  1110 

731 . 3750 

780. 4530 

7=-’ 

816.  'Ml 

61 4.6007 

312. 7623 

^34 

842.  3560 

340.  8691 

839.  5C 

927.7436 

924. 6402 

921 . 3983 

7lj 

931.2371 

923. 3024 

925.  2: 

1089.6810 

1088. 0641 

1050.9953 

1C:  >.».  8398 

1055.266" 

1051.2086 

7“ 

1210. 1996 

1209. 32 

1204. 17.  v 

7«; 

1216.2043 

1709.  8? 

1204.  1621 

1394. 8190 

1386.41 

1378. 9899 

70 

1394.  81 9 1 

1386.41-. 

1378.9899 
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TABLE  III — Continued 


energy  levels  of  Table  III  show  that  they  obev  the  rule 


where  k  =  0.5307  — *  0.5275.  This  results  directh  from 


0.5307 


with  the  change  in  rotational  energy  due  to  isotopic  substitution  given  bv 


\n  + 


and  A/I  ~  0  for  an  on  axis  isotopic  substitution.  As  a  result  at  high  K 
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TABLE  V 

Values  of  ihe  Fraley  Kao  Constant  for  Observed  Microwave  Transitions 


Transition 

k 

o 

I 

cT 

0.  530807 

9  —  9 

HI  “02 

0. 532543 

“3  —  9 

M3  “20 

0.532021 

''’M-321 

0. 536887 

423  “  330 

0. 531228 

515~  422 

0. 493528 

et6~  ‘23 

0.  523757 

dE,  dA  =  (Pa)  is  large  and  <)/■■  dC  ~  (Pc)  is  small,  k  ~  0.5307.  Conversely  at 
low  A'_i,  k  ~  0.5275. 

Equation  (2)  can  be  rewritten  for  transitions  as 

—  Els  EW‘  —  Ei6u 

- - =  k‘  +  (*“  -  kl) -  (4) 

I'is  —  ■'le  ^is  —  E16 

where  the  superscripts  refer  to  the  upper  and  lower  energy  levels  of  the  transition. 
The  correction  term  is  large  in  regions  where  the  dominant  terms  of  Eq.  (3)  are 
changing,  especially  for  transitions  whose  frequency  changes  little  upon  isotopic- 
substitution.  The  lines  of  Table  Y  which  show  the  greatest  deviations  from  Eq.  (1) 
are  the  same  lines  that  Eq.  (4)  predicts  to  have  these  variations. 

COMPARISON  WITH  OTHER  ENERGY  LEVELS  AND  CONCLUSIONS 

Our  original  1 1  TO  analysis  was  based  upon  our  measured  microwave  data  and 
distortion  constants  calculated  from  the  ITlrO  and  I1YO  analyses.  A  comparison 
between  the  energy  levels  which  resulted  from  this  and  those  of  Table  111  show 
exceedingly  close  agreement  (typically  0.002  cm  ')  up  to  the  5u  level.  At  higher  J, 

the  agreement  is  still  good  (typically  0.002  cm  1  at  low  K  ,  up  to  0.02  cm  1  at 

high  K  i)  except  for  the  6fl,  and  66(i  levels.  These  are  lower  by  0.17  cm  1  in  our  latest 
analysis.  This  is  because  our  earlier  analysis  of  HYO  was  perturbed  by  infrared 
energy  levels  for  6<u  and  66i  which  were  high  by  0.30  cm  1  compared  to  the  I I  YO 

TABLE  VI 

Comparison  of  Observed  Energy  Levels  with  those  Calculated  from  Eq.  (2)  (cm-1) 

Difference 
0.0000 
-0.0001 
-0.0001 
-0.  0004 
-0.0014 
-0. 0019 
-0.0007 
-0. 0099 


State 

Observed 

•  ; 

l10 

42.  1870 

42.  1870 

220 

135.4312 

135.4313 

330 

283. 7677 

7^3.  7  6"  8 

4  40 

485.2361 

485. 2365 

^50 

37.623! 

737.  6239 

660 

1038. 7648 

1038.  7667 

770 

1386.4179 

1336.  •  - 

3  80 

1778.  2960 

1778.  3059 
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energy  levels  shown  in  Table  111.  For  medium  values  of  J,  the  energy  levels  of  Refs, 
(hi  and  (<V)  differ  from  the  levels  of  Table  III  by  substantially  more  (typically 
0.01  cm-1)  than  do  our  earlier  energy  levels.  The  energy  levels  of  Ref.  ( 6 ),  which 
extend  to  higher  ./,  than  those  of  Ref.  (<?),  deteriorate  somewhat  relative  to  the  values 
of  Table  Ill,  but  all  agree  to  within  0.05  cm  1  or  better. 

It  is  to  be  expected  that  our  latest  analysis  is  better  than  previous  analyses  because 
it  contains  the  earlier  data  sets  as  subsets  and  because  our  analysis  technique  makes 
possible  rather  sensitive  tests  for  bad  data  points  and  substantially  reduces  the  degrees 
of  freedom.  It  is  perhaps  somewhat  surprising  that  our  earlier  analysis  was  so  good. 

Since  all  of  the  energy  levels  of  the  three  isotopes  were  calculated  independently, 
one  check  of  their  accuracy  and  also  of  the  accuracy  of  Eq.  (2)  can  be  accomplished 
by  the  use  of  the  IE1  "O  and  HslsO  energy  levels  in  this  equation  to  calculate  the 
energy  levels  of  iF170.  Since  it  is  most  difficult  to  get  good  energy  lex’els  at  high  ALt 
and  since  the  isotopic  splittings  of  the  energy  levels  are  also  greatest  there,  the  most 
stringent  test  would  be  the  Jj,<  levels.  The  close  argument  shown  in  Table  VI  between 
the  energy  levels  calculated  directly  from  the  H>l70  data  and  the  levels  calculated 
via  Eq.  (21  confirms  both  the  accuracy  of  this  relation  and  the  quality  of  the  energy 
levels  for  all  three  species. 

Received:  November  18,  lb 7 7 

Mole  added  in  proof.  The  problems  in  our  earlier  analysis  of  I  D'H )  have  recently  been  attributed 
(J.-M.  l'laud,  C.  Came  1‘evret,  and  R.  A.  Toth, Mol.  Speclrose.  65,  216  (,1977))  to  the  heavy  weights 
assigned  the  microwave  data  in  that  analysis.  In  fact,  as  shown  above,  the  errors  in  the  literature  of 
0J0  cm  1  in  the  infrared  energies  of  and  h6i  perturfted  several  of  the  higher  order  distortion  constants 
and  some  of  the  energy  levels  closely  associated  with  06o  and  6m.  All  analyses  reported  in  this  paper  con¬ 
tain  heavily  weighted  microwave  data,  and.  as  demonstrated  above,  excellent  agreement  among  the 
isotopes  result.  All  checks  of  isotopic  agreement  were  performed  after  each  isotopic  analysis  was  final, 
and  no  adjustments  at  subjective  points  in  the  data  analysis  were  carried  out  with  isotopic  agreement  as 
an  objective. 
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Millimeter-Wave  Spectrum,  Centrifugal  Distortion 
Analysis,  and  Energy  Levels  of  HNCV 
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One  hundred  and  eleven  new  rotational  transitions  of  HNO;,  have  been  measured  in 
the  millimeter-microwave  region.  This  data  set  is  the  basis  for  a  centrifugal  distortion 
analysis  and  for  the  calculation  of  a  complete  set  of  energy  levels  through  J  =  50.  The 
constants  which  result  from  this  analysis  are  (in  megahertz):  A  =  13  011.0287  ±  0.0057. 

H  =  12  099.8611  ±  0.0057.  C  =  6260.6391  ±  0.0006.  A,  =  (14.038  ±  0.026)  x  10  Aja 
=  (-20.1780  ±  0.0037)  x  1 0‘,  A*  =  (7.4153  ±  0.0108)  x  10  h,  =  (1.1828  ±  0.00039) 
x  IO-\  K  =  (-20.5648  ±  0.0046)  x  10  \  H,  =  (  -9.84  -  3.82)  x  10  \  =  (-9.933 

±  0.367)  x  10  ",  f/hJ  =  (1.03  ±  0.12)  x  10  :.  h,  =  (-9.241  ±  0.127)  x  10  “.  hJK 
=  (-1.398  ±  0.032)  X  10\  and  /i„  =  (1.135  ±  0.027)  x  10 

INTRODUCTION 

Nitric  acid  is  a  near-oblate  asymmetric  rotor  with  rotational  constants  of  the 
order  of  10  GHz  and  moderate  centrifugal  distortion.  The  presence  of  this  species 
in  the  atmosphere  and  its  basic  role  in  chemistry  have  stimulated  numerous 
spectroscopic  studies.  The  rotational  spectrum  of  HNO.,  has  been  studied  in  the 
centimeter-wave  region  (/-.?).  and  Kaushik  and  Venkateswarlu  ( 4 )  have  per¬ 
formed  a  centrifugal  distortion  analysis  of  these  data.  They  noted  that  the  rela¬ 
tively  large  standard  deviations  of  the  resulting  coefficients  and  the  indeter¬ 
minacy  of  must  be  ascribed  to  the  data  set  available.  Conventional  infrared 
spectroscopy  reveals  only  broad  vibrational  bands  because  of  the  dense  rota¬ 
tional  fine  structure  of  HNO:t  (5-8).  Recently,  however.  Brockman  ct  at.  (V)  have 
published  an  unassigned  diode  laser  spectrum  of  HNO;i  in  the  I  l-pim  region  that 
clearly  resolves  the  rotational  fine  structure.  Several  other  diode  laser  studies 
of  HNO,  are  also  in  progress  (10). 

In  this  paper  we  report  the  measurement  of  1 1 1  rotational  transitions  in  the 
region  between  80  and  300  GHz.  This  spectral  coverage  results  in  a  well- 
conditioned  data  set  which  makes  possible  an  accurate  centrifugal  distortion 
analysis  and  calculation  of  energy  levels  through  ./  =  50.  These  ground-vibra- 
tional-state  energy  levels  can  serve  as  a  basis  for  the  analysis  of  the  complex 
rotation- vibration  spectra  which  are  now  being  recorded  as  well  as  for  the  study  of 

I  1  This  work  was  supported  by  Army  Research  Office  Grant  DAAG29-77-G-0007. 

-  Present  address:  Laboratorio  di  Spettroscopia  Molecolare.  Consiglio  Nazionale  delle  Ricerche, 
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»  1  Nato  Senior  Postdoctoral  Fellow. 
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HNO:,  by  direct  microwave  techniques.  A  preliminary  report  of  this  work  has  been 
given  recently  (//). 

II  EXPERIMENTAL  DETAILS 

We  have  previously  detailed  our  general  experimental  technique  {12.  13).  A 
brief  description  of  the  specific  configuration  for  this  work  follows.  Millimeter- 
wave  energy  was  produced  by  King  and  Gordy  {14)  crystal  harmonic  generators 
driven  by  OKI  klystrons  in  the  35-GHz  range.  This  energy  was  focused  by  quasi- 
optical  techniques  through  a  4-m-long  absorption  cell  made  of  10-cm-diameter 
KIMAX  glass  pipe  and  detected  by  a  1.6-K  InSb  photoconduction  detector. 
The  HNO.,  vapor  was  taken  directly  from  a  standard  laboratory  mixture  of  con¬ 
centrated  nitric  acid. 

Although  most  of  the  lines  reported  here  are  strong  video  oscilloscope  lines, 
source  modulation  and  lock-in  techniques  were  used  to  record  the  weaker  high-./ 
lines.  For  this  work  the  35-GHz  klystron  was  phase  locked,  via  a  phased- 
locked  X-13  Varian  klystron  and  transfer  oscillator,  to  a  crystal  oscillator  whose 
phase  was  continuously  monitored  against  WWVB. 

Ill  RESULTS  AND  DISCUSSION 

We  have  used  Watson's  reduced  centrifugal  distortion  Hamiltonian  {15)  and 
the  computational  and  statistical  techniques  that  we  have  previously  discussed 
{16.  17)  for  the  analysis  of  the  rotational  spectrum  of  HNO;,.  A  "bootstrap  ' 
assignment-analysis  procedure  was  used.  At  each  step  lines  that  would  provide  a 
balance  between  a  maximum  of  new.  independent  information  for  the  analysis 
and  a  minimum  risk  of  assignment  error  were  selected  for  measurement.  In 
practice  this  amounted  to  the  selection  of  lines  whose  prediction  uncertainties 
were  several  megahertz.  This  procedure  was  iterated  until  all  of  the  approximately 
3000  lines  of  significant  strength  below./  =  50  were  predicted  with  an  uncertainty 
of  I  MHz  or  less. 

Table  I  shows  the  131  rotational  transitions  included  in  our  analysis  and  Table 
II  shows  the  spectral  constants  which  result.  In  our  analysis  Hk  was  found 
to  be  only  marginally  determined  and  was  subsequently  eliminated  from  the 
constant  set.  The  rms  deviation  of  the  analysis  was  0.124  MHz  with  or  without 

TABLE  II 

Spectroscopic  Constants  of  HNO;,  (MHz) 

a  .o.hi':i,  .  <-■:  c  c  2eo.  cast  «  o.ocor 

K  -20.  !790t  0.0037)- to*5  /y  =  f7.  4!53t  0.010S).  to*3 
0. 00039) .  to*3  6k  ■  (-20.  0.0040).  10*3 

Hjk  -  I-9.S33  t  0.307).  to*8  !!,„  -11.03*  0.  12)-  10‘7 
hJK  -1-1. 398  *  0.032).  to*  ty  =  (1.  1354  0.  027).  10*8 
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//A  as  a  variable.  The  spectral  constants  have  been  used  to  generate  the  energy 
levels  shown  in  Table  III.  Because  of  the  large  number  of  rotational  lines,  the 
calculated  spectrum  is  not  reproduced  here.  It  is  available  either  from  the  Editorial 
Office  of  this  Journal  or  from  the  second  author. 

An  important  point  should  be  made  about  the  use  of  these  data  in  conjunc¬ 
tion  with  other  data  sets,  especially  infrared  rotation- vibration  data  sets. 
Although  accurate  ground-state  energy  levels  are  a  convenient  starting  point 
for  spectra  assignments,  from  a  statistical  point  of  view  it  is  preferable  to  use  the 
observed  microwave  data  directly  in  the  final  spectral  analysis.  The  appropriate 
weightings  for  data  are  then  reasonably  straightforward.  On  the  other  hand,  if 
energy  levels  or  spectral  constants  are  used,  the  appropriate  statistical  treat¬ 
ment  requires  the  consideration  of  both  the  uncertainties  in  the  spectral  constants 
and  in  the  correlation  matrix  (/#).  This  is  especially  important  in  large  centrifugal 
distortion  analyses  because  of  the  significant  correlations  among  the  constants. 

Received:  August  14.  1978 
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Millimeter  wave  measurements  of  the  rotational  spectra  of 
CIF,  BrF,  BrCI,  ICI,  and  IBr 

Robert  E.  Willis,  Jr.a)  and  William  W.  Clark,  lllb> 

Department  of  Physics,  Duke  University,  Durham,  North  Carolina  27706 
(Received  7  December  1979;  accepted  18  January  1980) 

The  rotational  spectra  of  all  twelve  stable  isotopic  species  of  CIF,  BrF,  BrCI,  ICI,  and  IBr  were  observed 
and  measured  in  the  millimeter  wave  region  by  means  of  a  sensitive  microwave  spectrometer.  Transitions 
were  detected  over  a  wide  range  of  frequencies  for  molecules  in  both  the  ground  vibrational  state  and 
several  excited  states.  The  rotational  spectrum  of  each  molecule  was  split  by  the  nuclear  quadrupole 
interaction.  Altogether,  250  new  lines  were  measured.  These  correspond  to  136  pure  rotational 
transitions.  Values  of  the  Dunham  coefficients  y„„  Y,„  Y2I,  Y}„  Ym,  Y:1,  and  T0J  were  obtained  from  a 
computer  analysis  of  the  measured  frequencies.  From  these  coefficients  a  number  of  equilibrium  constants 
were  derived  to  significantly  greater  accuracy  than  in  previous  work.  In  particular,  the  equilibrium 
distance,  rr,  was  found  to  two  or  three  more  significant  figures. 


INTRODUCTION 

The  mixed  halide  diatomic  molecules  described  in  this 
report  have  all  been  studied  previously  through  tech¬ 
niques  of  microwave  spectroscopy  by  various  workers1-8 
whose  results  are  presented  in  the  N.B.S.  tables  of 
diatomic  molecules.9  However,  most  of  their  work  was 
confined  to  the  centimeter  wave  region  and  was  con¬ 
cerned  with  nuclear  quadrupole  splittings  of  a  few,  low, 
rotational  levels.  Consequently,  only  a  small  number 
of  rotational  constants  have  been  reported.  The  present 
work  in  the  millimeter  wave  region  has  generated  mea¬ 
surements  of  136  new  rotational  transitions  for  mole¬ 
cules  in  a  number  of  excited  vibrational  states.  This 
has  increased  both  the  quantity  and  the  accuracy  of  the 
rotational/vibrational  constants  which,  in  turn,  has  per¬ 
mitted  the  derivation  of  the  equilibrium  constants  Be, 

Ci>e,  u„xe,  and  re  as  well  as  the  Dunham  potential  con¬ 
stants  a0,  alt  n2,  and  n3.  In  addition,  this  work  has 
yielded  rotational  constants  which  are  significantly  dif¬ 
ferent  from  previously  reported  values  for  six  out  of  the 
12  isotopic  species  studied. 

THE  SPECTROMETER 

The  spectrometer  used  for  this  study  was  similar  in 
many  respects  to  the  ones  used  in  other  investigations 
of  this  laboratory. 10,11  Briefly,  it  consisted  of  a  kly¬ 
stron-driven  crystal  harmonic  generator,  a  quasi-free- 
space  absorption  cell,  and  a  helium -cooled  InSb  photo¬ 
detector.  Most  of  the  transitions  were  measured  in  real 
time  on  an  oscilloscope  screen,  although  a  few  were  ob¬ 
served  on  a  chart  recorder  with  a  lock-in  amplifier  and 
a  phase -locked,  slow -sweep  klystron. 12 

For  these  experiments  the  absorption  cell  was  a  con¬ 
tinuously  evacuated  2  in.  x  4  ft.  nickel  pipe  with  Teflon 
end  plates  to  permit  entry  of  the  microwave  radiation. 
The  cell  was  encased  by  two  1400  W  Linberg  2724 -KSP 
semicircular  heating  elements.  This  arrangement  al¬ 
lowed  the  temperature  of  the  cell  to  be  raised  so  that 
excited  vibrational  states  could  be  populated.  The 
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flanges  and  end  plates  were  sealed  with  neoprene  O  rings. 
Because  of  reactions  with  the  halogens,  these  rings 
had  to  be  regreased  or  replaced  frequently.  A  nitrogen 
cold  trap  was  used  to  protect  the  forepump  from  the 
gases;  nevertheless,  the  oil  had  to  be  changed  occasion¬ 
ally,  especially  during  the  use  of  iodine  compounds. 

CHEMICAL  PREPARATIONS 
CIF  and  BrF 

Chlorine  trifluoride  C1F3  was  used  in  the  preparation 
of  CIF  and  BrF.  This  is  an  extremely  reactive  gas13 
even  at  the  pressures  used  in  the  experiment  (~  5*  10"2 
Torr);  therefore,  a  few  special  precautions  were  taken 
for  the  fluorine  compounds:  (1)  The  neoprene  O  rings 
were  wrapped  with  Teflon  tape  and  (2)  copper  tubing  and 
brass  fittings  were  used  for  connectors  in  place  of 
rubber  or  glass  tubing. 

Chlorine  monofluoride  CIF  was  observed  as  a  decay 
product  of  C1F3  as  the  latter  was  was  simply  allowed  to 
flow  through  the  cell.  This  method  was  chosen  over  the 
use  of  CIF  gas  directly  because  C1F3  is  less  expensive, 
easier  to  handle,  and  somewhat  safer  than  CIF  gas. 
However,  the  reaction  C1F3-C1F  +  F2  proceeds  very 
slowly  at  the  low  pressures  of  the  absorption  cell.  As  a 
result,  the  CIF  lines  were  relatively  weak.  Adding 
chlorine  gas  to  the  stream  caused  the  line  strength  to 
go  down  rather  than  up.  Warming  the  cell  to  about 
100  °C  increased  the  reaction  rate,  but  above  100  °C, 
the  line  strength  decayed,  dropping  to  1/4  of  the  maxi¬ 
mum  strength  at  only  200  °C.  The  rapid  signal  decay  at 
elevated  temperatures  made  it  difficult  to  observe 
higher  vibrational  states.  This  problem  was  common 
to  all  the  molecules  in  this  study  but  was  most  severe 
for  CIF  because  of  its  small  molecular  weight. 

Bromine  monofluoride  BrF  has  been  found  only  in  the 
gaseous  state  in  equilibrium  with  Br2  and  BrF3.  For 
this  reason  BrF  was  made  by  simultaneously  pumping 
room -temperature  quantities  of  Br2  and  C1F3  through 
a  cell  that  varied  in  temperature  from  25  to  450  °C. 
Optimum  line  strength  occurred  with  the  cell  at  200  °C, 
and  the  lines  disappeared  at  500  °C.  Temperatures  as 
high  as  450  "C  were  used  to  populate  excited  vibrational 
states.  Often  strong  lines  were  observed  when  the  bro- 
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mine  flow  was  off.  This  agrees  with  the  fact  that  in  an 
equilibrium  mixture  of  Br2,  BrFj,  and  BrF,  the  rela¬ 
tive  fraction  of  BrF  drops  off  sharply  at  Br2  is  added. 

BrCI,  ICI,  and  Brl 

The  compounds  not  involving  fluorine  were  all  easier 
to  make  and  handle  than  the  above  two.  Bromine  mono¬ 
chloride,  BrCI,  was  made  when  one-third  mole  each  of 
bromine  and  chlorine  were  combined  in  a  22  1  Pyrex 
flask.  At  room  temperature  such  a  mixture  will  be 
57%  BrCI.  The  BrCI,  which  boils  at  20  °C,  was  then 
allowed  to  flow  through  the  cell.  The  ground-state  lines 
were  strongest  at  room  temperature  and  200  °C  was  the 
highest  temperature  used  to  populate  excited  vibrational 
states. 

Since  the  vapor  pressure  of  iodine  is  too  high  for  it  to 
be  placed  directly  in  the  cell,  the  iodine  compounds 
were  also  made  in  the  22  1  flask.  Addition  of  one -third 
mole  of  chlorine  to  excess  iodine  resulted  in  a  rapid 
exothermic  reaction  producing  ICI,  a  red  liquid  or  gas 
(melting  point  27  °C),  as  well  as  IC13,  a  dark  yellow 
solid.  The  ICI  was  then  allowed  to  flow  through  the  cell. 
All  measurements  were  made  at  room  temperature  be¬ 
cause  ICI  decomposes  at  97  °C. 

Iodine  monobromide  IBr  was  made  in  a  similar  man¬ 
ner  from  a  mixture  of  iodine  and  bromine  in  the  22  1 
flask.  This  produced  IBr,  dark  grey  crystals  that  melt 
at  41  °C  and  sublime  at  50  °C,  as  well  as  IBr3,  a  brown 
liquid.  The  flask  was  heated  slightly  with  an  infrared 
lamp  to  produce  flow  rate  sufficient  for  measurement. 
The  IBR  lines  did  not  appear  until  the  IBRj  was  pumped 
out  of  the  flask.  Since  IBR  decomposes  at  116  °C,  the 
cell  was  not  heated.  This  did  not  prevent  observation 
of  higher  vibrational  states,  since  in  molecules  as  heavy 
as  IBr,  excited  vibrational  states  are  well  populated  at 
room  temperature. 

ANALYSIS  OF  ROTATIONAL  SPECTRA 

The  most  accurate  formulas  for  analysis  of  the  rota¬ 
tional  spectra  of  diatomic  molecules  are  provided  by 
Dunham  theory.14  Dunham  solved  the  Schrodinger  wave 
equation  for  the  diatomic  vibrating  rotor  with  a  potential 
function  of  the  form 

F(£)  =  cha0t,z(l  +a,£  +  a2£2  +•••),  (1) 


vv,J~  2fT'oi  ♦  >'ii(r  ♦  J)  *  T21(»  <  l)1  ♦  Fj,(r  ♦  i?  +■  ••  I 

x  U  +  1)  +  4( Fo«  +  l',j(r  ♦  J) - |(^  ♦  1 13 

♦  >'„(./ ♦l)1f(./*2)>-.7'|*  •••  ,  (4) 

where  only  terms  which  give  a  detectable  contribution  to 
the  spectra  of  tKs  report  are  retained. 

Tht  one;  irum  of  each  species  was  analyzed  by  sub¬ 
stitution  of  each  reduced  experimental  rotational  fre¬ 
quency  (see  quadrupole  analysis)  and  the  appropriate  J 
and  v  quantum  numbers  into  Eq.  (4).  Hus  produced  a 
set  of  simultaneous  equations  which  were  solved  by  the 
method  of  least  squares  for  the  K,  „’s.  Although  it  is 
possible  to  analyze  each  individual  isotopic  species  in¬ 
dependently,  use  of  the  Dunham  isotope  relationship 

(5) 

allows  all  the  isotopic  data  for  each  molecular  species 
to  be  analyzed  together.  This  results  in  significantly 
better  values  for  the  Y,  „,  especially  in  the  case  of  less 
abundant  species  for  which  few  transitions  are  observed. 
Equation  (5)  is  accurate  to  within  experimental  error  for 
all  Ylm  except  Tol 


Tbe  K03  was  calculated  from  the  Dunham  relation 


and  its  contribution  to  the  measured  frequencies  re¬ 
moved  before  the  least-squares  analysis.  Although  its 
value  can  be  calculated  directly  from  the  spectra,  the 
procedure  used  here  results  in  a  more  accurate  y03  as 
well  as  improved  values  for  the  other  constants. 

After  determination  of  the  Dunham  coefficients  the 
first  four  potential  constants  of  Eq.  (1)  can  be  found  by 
the  relations 


_  _ 
aoSW.  =  ~De 


o  1  : 


-1  = 


_ 1  n  i 

3(->o2T01)i/z 


(7) 

(8) 
(9) 


°2  =  (_“f3‘2)  +  *  +  77,)  +  T  . 

'73  =  -l^i  +  Ti(3  +  13t7,)  -^[4  +377,(1  +77,)1  -  1.  (10) 
D  1 02  5  £ 


where  £  =  (r  -  r,)/r„  and  770,  t i,,  t72,  ...  are  the  potential 
constants.  His  solution  for  the  rotation-vibration  ener¬ 
gies  is  expressed  by  the  series 

Ev,,  =  hY,  TI>m(i-+ |)'kU+l)l"  ,  (2) 

l  ,m 

in  which  the  coefficients  Ylm  are  known  as  the  Dunham 
coefficients.  Each  Y,  „  is  a  function  of  the  molecular 
parameters  r„  770,  77,,  t72,  ...  which  appear  in  Eq.  (1). 
In  terms  of  the  familiar  spectroscopic  constants, 

Tio  -  uv  l'o,  =  Bj  F,,  a  -  a,-  F2,  =  y,-  , 

T2o  =  —  <+'*  v*'  y<a  —  -  fV  T,2  =  -  0,-  F03  =  H,  .  (3) 

For  pure  rotational  transitions  the  transition  frequen¬ 
cies  are  given  by 


The  relations  given  in  Eq.  (3)  and  subsequently  used 
in  Eqs.  (7)-(10)  are  not  exact.  The  exact  relations  are 
power  series  in  (R,/u',)2.  Since  for  the  mixed  halides 
this  ratio  is  of  the  order  10‘6,  these  relations  may  be 
considered  accurate  to  the  experimental  uncertainty  ex¬ 
cept  for  F01,  which  can  be  calculated  to  a  higher  accu¬ 
racy  from  precision  spectral  data.  The  precise  rela¬ 
tionship  between  1'01  and  Be  is 

Y0i=  Bjl  +C3(R€/u),)2l  =  Rf  +  fi3  .  (11) 

This  difference  between  Be  and  1’01  is  known  as  the  Dun¬ 
ham  correction  (?3). 

Two  other  small  corrections  to  I'0,  result  from  elec¬ 
tronic  effects  which  are  not  included  in  the  Dunham 
theory.  These  arise  from  the  nonspherical  distribution 
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of  the  atomic  electrons  around  their  respective  nuclei 
(5,)  and  the  wobble -stretch  effect  (6,).ls,le  Hie  first  of 
these  may  be  calculated  from  knowledge  of  the  molecu¬ 
lar  g  factor  by 

6 1  =  —  (m ,/tnp)gj  B,  .  (12) 

Although  fl.  cannot  be  calculated  directly  from  known 
molecular  parameters,  its  reduced  mass  dependence  of 
(1  q)2  allows  its  calculation  if  two  or  more  isotopic 
species  are  observable.  Since  J>,  and  P3  also  have  a 
U  q)2  dependence,  the  net  effect  of  all  three  correc¬ 
tions  may  be  calculated  from  the  system  of  equations. 

(B.-  >'0I),  -C/pf  =  0  ,  (13) 

where  C  -  P,  6.  v  P3  and  i  is  an  isotopic  species  label. 

The  values  of  Be  thus  calculated  provide  a  very  ac¬ 
curate  determination  of  the  equilibrium  distance  by 
means  of  the  definition  of  By. 

B,=  It  8i72p>'  .  (14) 

The  zero-point  vibration  frequency  iv,  can  next  be  cal¬ 
culated  from  the  relation 

u>,  (4 Bl  f;,)t/!  '  f4B8  (-  r<e)|,/!  (15) 

and  u>,v,  is  obtained  from 

u'.v,  3 1  o\(fi3  —  i (l\ )  2  . 

NUCLEAR  QUADRUPOLE  INTERACTIONS 

Hyperfine  structure  which  occurs  in  the  diatomic 
halides  causes  tlie  observed  frequencies  to  deviate 
slightly  from  those  given  in  Eq.  (4).  These  deviations 
are  caused  by  nuclear  electric  quadrupole  interactions 
with  the  electric  field  gradient  at  the  nucleus.  This 
splits  the  rotational  energy  levels  into  2/  *  1  sublevels, 
where  /  is  the  nuclear  spin.  The  actual  energy  levels 
are  found  from  perturbation  theory  and  to  first  order 
are  given  by1' 

Ej.,., f-  Ej.v-rQnY(J,l,F)  .  (16) 

where  E j  is  the  rotational  energy  given  by  Eq.  (2), 
ri/Q  is  the  quadrupole  coupling  constant,  and  Y(J,l,  F ) 
is  Casimer’s18  function.  The  F  in  Eq.  (16)  is  the  total 


i  l 

Kli;.  1.  I’ompMrison  of  Che  theoretical  ami  experimental 
quail  rupoL*  structure  of  'PrK.  «/  9  *10,  r  2,  at  -09979.  S-l 
Mil/. 


TABLE  1.  Dunham  coefficients  and  equilibrium  con¬ 
stants  of  chlorine  monofluoride. 


**Cl'*F 

s,Cl'*F 

t01  (MHz.) 

15483.2702(36)* 

15  188.8247(30) 

T„  (MHz) 

-129.6497(50) 

-  125.9087(49) 

T,,  (kHz) 

-520.3(1.7) 

-500.7(1.0) 

T0;  (kHz) 

-26.91(2) 

-25.  90(2) 

T,j  (ktl/.) 

-0.  u<2) 

-0.13(2) 

Tj3  (Hz) 

-  5. 47(3)  v  10-’ 

-5. 16(3)  x  It)"5 

8,  (MHz) 

15  484.4304(70) 

15  189.  9412(70) 

r,  (A) 

1.62S341(4)b 

1.628341(4)*’ 

(cm"1) 

783.553(30) 

776.007(30) 

(cnT1) 

5.045(87) 

5.020(84) 

al  neerta  inties  represent  one  standard  deviation. 
*rhe  uncertainty  accrues  from  Planck’s  constant. 


angular  momentum  quantum  number  which  can  assume 
tlie  values  (J  +  1).  (Jt/-I) . IJ  -  /I.  Each  rota¬ 

tional  transition  is  thus  split  into  a  number  of  hyperfine 
lines,  some  of  which  are  doubly  degenerate  or  nearly 
degenerate.  Figure  1  shows  a  typical  rotational  transi¬ 
tion  which  is  split  by  the  quadrupole  interaction. 

For  J  —  </  +  l  rotational  transitions,  the  selection 
rules  for  quadrupole  structure  are  AF'-O,  ±  1.  For  the 
high  J  values  observed  in  this  report,  the  SF  -  +  1 
transitions  dominate  in  strength.  Therefore,  the  devia¬ 
tion  in  line  frequency  from  Eq.  (4)  is  given  by 

Ai - iV/VflVt  1,/,  D-1V./.F)]  .  (17) 

Since  the  frequency  shift  given  in  Eq.  (17)  decreases 
as  1  J 2,  tlie  quadrupole  structure  was  unresolved  at 
high  values  of  J.  Nevertheless,  quadrupole  structure 
was  observed  in  lower  J  transitions  for  all  of  the  mole¬ 
cules  in  this  report  and  had  to  be  considered  before  the 
rotational  analysis  was  attempted. 

The  frequencies  of  the  unsplit  rotational  lines  were 
determined  in  two  ways.  First,  through  measurement 
of  the  hyperfine  splitting  in  this  work,  a  reduced  fre¬ 
quency  could  be  determined  from  Eq.  (17).  Tlie  second 
method  was  to  use  the  n/Q's  obtained  by  the  earlier 
lower  frequency  work,  which  yielded  values  one  to  two 
orders  of  magnitude  more  accurate  than  the  ci/Q's  cal¬ 
culated  from  this  work.  This  latter  method  gave  a 
slightly  different  center  frequency  for  each  hyperfine 
component;  the  average  of  these  frequencies  was  then 
averaged  with  the  frequency  resulting  from  the  first 
method  to  yield  a  final  value  for  the  reduced  rotational 
frequency. 

For  the  species  in  which  both  nuclei  have  a  quadrupole 
moment,  the  weaker  quadrupole  moment  produced  no 
observable  splittings  or  shifts  but  did  cause  consider¬ 
able  broadening  in  some  cases.  Tlie  effect  was  strong¬ 
est  for  low-frequency  transitions  in  IBr  and  BrCI. 
Severely  broadened  transitions  which  could  not  be  re¬ 
solved  were  avoided  in  most  instances  because  of  diffi¬ 
culty  in  determination  of  the  center  frequency. 

SPECTRAL  CONSTANTS 

The  actual  transition  frequencies  are  not  listed  here 
but  are  recorded  in  the  Ph.D.  thesis  of  R.  E.  Willis. 19 
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TADl.K  11.  Dunham  coefficients  and  equilibrium 
constants  of  bromine  monofluoride. 


,,6r,,F 

“Br**!’ 

i’a,  (MHz) 

10  687.4502(35)* 

10  010.  8470(80) 

*•„  (MHz) 

—  77.  8029(79) 

-77.3040(78) 

Y:x  (kHz) 

-  201 . 0(45) 

-  199.7(45) 

i’j,  (kHz) 

-0.  2(8) 

-0.  1(8) 

I'd;  tkll/i 

-  12.  049(10) 

-  1  1.084(10) 

T,;  (Hz) 

-66(1) 

-00(11) 

I'm  (Hz) 

-3.  937(7)  x  10° 

—  3.  881(7)  v  10's 

Bt  (MHz) 

10  668.  0970(69) 

10  010.  0882(00) 

r,  (A) 

1.758987(4)* 

1. 758987(4)* 

s-,  tom'1) 

000.  070(28) 

068.07(28) 

s',  V,  tcm'1) 

8.  808(70) 

8.  844(78) 

*1  ncertainties  listed  are  one  standard  deviation. 
bl'he  uncertainty  accrues  from  Planck’s  constant. 


The  molecular  constants  are  given;  they  can  be  used  for 
calculation  of  the  frequencies  to  high  accuracy  C  0.1  MHz\ 

CIF 

Previous  microwave  measurements  of  CIF  have  been 
made  by  White*  and  Gilbert.2  They  measured  a  number 
of  quadrupole  lines  of  the  J  -  0—  1  rotational  transition 
in  the  0. 1  vibrational  states.  The  current  work  has 
produced  measurements  of  42  new  transitions:  32  of 
,5C1F  and  10  of  the  lesser  abundant  S7C1F.  These  mea¬ 
surements  yielded  22  reduced  rotational  frequencies 
covering  the  rotational  states  from  J  -  2  up  to  J  1 1  in 
the  0. 1,  2  vibrational  states.  Hie  frequencies  ranged 
from  92  GHz  at  J  -  2  -  3  to  339  GHz  at  J  10  -  1 1 . 

The  Dunham  coefficients  obtained  from  these  measure¬ 
ments  are  shown  in  Table  1  along  with  the  equilibrium 
constants.  Comparison  of  calculated  frequencies  with 
the  measured  frequencies  yields  a  root  mean  square 
deviation  of  35  kHz  for  both  isotopes  combined. 

From  the  previous  measurements  and  an  infrared 
measurement  of  ay, 20  three  Dunham  coefficients  (F01, 

I  n.  were  obtained  (calculated  by  Ref.  9'.  These 
values  are  in  rough  agreement  with  our  values,  except 
for  of  57CIF  which  differs  by  5  MHz. 


TABLE  IV.  Dunham  coefficients  and  equilibrium 
constants  of  iodine  monochloride. 


,7,l’Vi 

,7TIs:C'l 

T„,  (Mill) 

3422. 3638(20)* 

8277.4173(10) 

rit  (MHz) 

- 15.9751(43) 

-  14.9710(40) 

y,,  (MHZ) 

-  35.  2(26) 

-32.3(24) 

r„  (kHz) 

-0.0(4) 

-0.5(3) 

To:  (kHz) 

-  1.2081(0) 

-  1.  108(0) 

(11?) 

-  5.  0(0) 

-  5.  5(6) 

1  os  ( 11/ ) 

-  2.639(3)  x  10'4 

—  2. 31 8(3)  x  UP4 

B ,  (MHz.) 

8422. 4275(38) 

3277 . 47 58(30) 

>,  lA) 

2.320962(6)* 

2.320962(6)* 

(cm"1) 

884 . 2  00  v 00) 

376.064(94) 

a1,  v,  tcm'1) 

1.501(20) 

l . 437(25) 

U  ncertainties  listed  are  one  standard  deviation, 
^ihe  uncertainty  accrues  from  Planck’s  constant. 


BrF 

Previous  microwave  measurements  of  BrF  have  been 
made  by  Smith  <•/  <tl.3  As  in  the  case  of  CIF,  only  the 
lowest  rotational  transition  in  the  0, 1  vibrational  states 
was  measured.  In  the  current  work  56  new  transitions 
have  been  measured:  20  of  79BrF  and  26  of  slBrF. 

These  measurements  produced  28  reduced  rotational 
frequencies  covering  the  rotational  states  from  J  3  to 
J  13  and  the  vibrational  states  0, 1,2,3.  Hie  measured 
frequencies  ranged  from  85  to  276  GHz  and  the  mis 
deviation  of  the  fitting  was  37  kHz. 

Table  11  shows  the  Dunham  coefficients  and  equilibri¬ 
um  constants  obtained  from  this  work.  The  earlier 
measurements  gave  values  for  1'01  and  Fu  (calculated 
by  Ref.  9l.  While  our  values  of  l-,,  are  in  agreement, 
our  values  of  If,,  differ  by  about  50  MHz  for  both  isotopes. 

BrCI 

For  BrCI,  63  new  transitions  were  measured:  28  of 
7aBr55Cl.  17  of  *'Br5SCl,  10  of  7*Br57Cl,  and  eight  of 
**Br57Cl.  The  number  of  reduced  rotational  frequencies 
obtained  was  33,  covering  rotational  states  from  ./-  13 
to  ./  -  34  and  vibrational  states  v  -  0  to  r  3.  The  fre¬ 
quencies  of  the  measurements  ranged  from  122  to  309 


VA1U  K  111.  Dunham  coefficients  and  equilibrium  constants  of  bromine  monochloride. 


;,HrsVl 

«Hrs»Cl 

’Dr 'VI 

’’Br’Vl 

(MHz) 

4570.  9087(13)* 

4530.3140(13) 

4399. 8510(13) 

4305. 1979(12) 

»  ,,  (MHz) 

-  23. 1420(29) 

-  22. 8800(29) 

-  21.8553(27) 

-  21 . 5970(27) 

1'.,  (kHz) 

—  01.9(18) 

-  01 . 0(18) 

-57.4(17) 

—  50. 5(10) 

(kHz) 

-0.9(3) 

-0.  9(3) 

—  0.  8(3) 

-0.  8(3) 

Ka:  (kHz) 

-2.1527(7) 

-2.  1202(7) 

-  1.9945(7) 

-  1 . 9632(7) 

r,:  diz) 

-  10.7(8) 

-  10.  5(8) 

-  0.  7(8) 

-  9.  5(8) 

I'm  illzi 

-4. 05013)  v  U>‘4 

-  4. 551(3)  M0*4 

-  4.  152(3)  v  10*4 

-  4  .  Oaatfi)  x  10'4 

fl,  tMllz) 

4571 . 0205(26) 

4530. 3793(25) 

4399. 9082(25) 

4365.  2527(25) 

>,  (A) 

2.  136124(5)* 

2.  136124(5)* 

2.  136124(5)* 

2.  136124(5)* 

s',  (cm*1) 

444.322(02) 

442.034(02) 

435. 920(01) 

•134.206(61) 

v,  tcm'1) 

1 . 850(20) 

t .  849(20) 

1.78'  19) 

1.771(21 

*1  pcertaintles  tiste.l  are  one  standard  deviation. 

*  i  he  uncertainty  accrues  from  Planck’s  constant. 
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1  API  V  V.  I'unhnm  coefficients  and  equilibrium 
constants  ot  iodine  monobromide. 


,:T  *Br 

'•  I4’  Hr 

1  MHz) 

1 703.  7 9S20(30'* 

1077.  S3 9o 7(30' 

V„  (MHz) 

-  5.  90413(41) 

-  76982(46) 

>  :i  (kHz' 

—  13. 00(13' 

-  13.23(13' 

T„  (kli/1 

-0.  14(2' 

-0.  13(2) 

T,v  (Hr) 

-304.91(2) 

-  293.70(2) 

111.  ) 

-  1 . 30(2' 

-  1.44(2) 

r„  (Hr) 

-  3.  9868111)  x  10' 

-  3.8075(11'  M0* 

(MHz' 

1703. SI  72(0' 

10  77 . 87 80(0' 

>.  v\) 

2.  179(178(61* 

2. 40907S(0'b 

w,  (enT1) 

20S .  09(9' 

200. 04(9) 

...  v,  lem'l 

0.  826(0' 

0.814(0' 

*l'noerlamiies  listed  are  one  standard  deviation. 
N'he  uncertainty  aeerues  from  Planck’s  constant. 


GHz.  and  the  rms  deviation  of  the  fitting  was  36  kHz. 

The  Dunham  coefficients  and  equilibrium  constants  of 
BrCI  are  shown  in  Table  111. 

Smith  <•/  ill.'  did  earlier  work  on  BrCI  as  well  as 
BrF.  As  for  the  previous  molecule,  they  measured 
a  number  of  quadrupole  lines  of  the  J  -  0—  1  rotational 
transition  in  two  vibrational  states  10.  If.  11ns  yielded 
two  Dunham  coefficients  for  each  species  )01  and  Vu 
(calculated  by  Ref.  91.  Except  for  one  instance,  their 
numbers  are  in  good  agreement  with  our  more  accurate 
rallies.  Their  value  of  F01  for  :,Br;,CI  differs  by  100 
MHz  from  ours.  This  may  be  a  misprint  because  it 
clearly  violates  the  isotopic  relationship  of  Eq.  15'  al¬ 
though  it  does  appear  in  both  the  original  article  and  in 
the  review  work  of  Ref.  9. 

ICI 

Weidner,5  Townes  and  Herbst  and  Steinmetz' 

have  all  studied  ICI.  These  earlier  studies  produced 
measurements  of  over  150  quadrupole  lines  m  three 
rotational  and  two  vibrational  states,  These  measure¬ 
ments.  in  turn,  have  given  very  accurate  values  of  1\„ 
and  l'„.7'a  The  current  work  on  1C1  has  yielded  five 
additional  Dunham  coefficients,  as  shown  m  Table  IV. 
lliese  were  obtained  from  50  new  quadrupole  lines  143 
of  lS5Cl  and  eight  of  l5V0  of  23  new  rotational  frequen¬ 
cies.  The  current  work  covered  the  rotational  states 
from  J  12  to  J  44  and  vibrational  states  r  0  to  r  -  3. 
The  frequencies  involved  ranged  from  85  to  300  GHz. 

The  rms  deviation  of  the  fitting  was  higher  for  1C1  (50 
kllzl  than  for  any  other  molecule  in  this  study  because  of 
difficulty  in  resolution  of  the  quadrupole  structure. 

IBr 

Like  ICI,  IBr  is  a  heavy  diatomic  with  a  complicated 
quadrupole  structure  at  low  .1.  .I.iseta"  has  measured 
over  150  quadrupole  lines  of  two  rotational  transitions 
(./  4-5  and  ,/  5-61  in  three  vibrational  states  (r  0, 
1,21.  in  this  study  we  have  measured  40  new  quadrupole 
lines  (19  of  l”Hr  and  21  of  I*‘Br'  in  30  rotational  states 
and  five  vibrational  states  (0,  1,2,3,41.  rhe  rotational 
suites  covered  the  range  from  ./  43  to  ,/  96  and  the 
frequencies,  from  150  to  325  GHz.  Die  transitions 
above  J  53  -  54  revealed  no  quadrupole  structure. 


i  API  I  Y|.  Dunham  potential  constants  ol  the  diatomic  mixed 
halides. 


(cm*1' 

rt: 

•a 

ill-’ 

2. 9712(261*  v  10' 

-3.  1171(9) 

3.03(11) 

-  8.  8(8) 

Hr  1-’ 

3.  1303(20'  v  10 

-  3.  2893(1  O' 

0. 29(13) 

-  9.5(H) 

Hn'l 

3.  2309(9)  v  10 

-3.4590(3) 

6. 86(8) 

-11.2(7) 

U'l 

3.2338(10'  M0S 

-  3.0188(3' 

7.00(13) 

-  18.2(1.:!) 

IHr 

...  1 7338(20'  v  10 

-3.7305(2) 

7.70(3' 

-12.6(3) 

*1  nccrtainties  represent  one  stands  I'd  deviation. 


In  Table  V  are  listed  the  Dunham  coefficients  obtained 
from  this  work.  Jaseta8  was  able  to  obtain  values  for 
i;„  and  l’n,  but  our  numbers  are  not  m  agreement  with  his. 

POTENTIAL  CONSTANTS 

Shown  m  Table  VI  are  the  Dunham  potential  constants 
as  calculated  from  Eqs.  (71  through  (10l.  To  within  ex¬ 
perimental  uncertainty,  these  constants  are  isotopically 
invariant,  a  result  which  confirms  an  initial  assumption 
of  Dunham  theory.  It  is  of  interest  to  note  that  the  val¬ 
ues  are  similar  for  all  the  molecules,  generally  increas¬ 
ing  slightly  with  increasing  mass. 
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ABSTRACT  Hyperfine  structure  due  to  :US  in  its  natural 
abundanee  of  0.76%  has  been  measured  in  the  electron  spin 
resonance  of  free  radicals  produced  by  x-irrad iat ion  of  single 
crystals  of  N-acely  l-l  -cysteine  at  77  K.  These  measurements 
proved  that  the  radicals  produced  at  77  K  n  ith  principal  g  r  allies 
of  1.990. 2.006.  and  2.214  arc  monosulfide  radicals  with  the  3/r 
unpaired  electron  density  of  0.70  on  the  S.  They  are  believed  to 
he  negatively  charged  molecules  ROILS  II  or  neutral  RCHjSIL 
radicals  in  which  90%  of  the  spin  density  of  the  captured  elec¬ 
tron  is  concentrated  in  a  </-p  hybrid  orbital  on  the  S.  As  the 
temperature  is  raised  to  300  K.  these,  as  well  as  the  carbon- 
centered  radicals  produced  at  the  lower  temperature,  are  mostly 
converted  to  neutral  disulfide  radicals  RCILSS  like  those  ob¬ 
served  in  irradiated  cystine. 

Free  radicals  w  ith  electron  spin  density  concentrated  on  sulfur 
(sulfur-centered  radicals)  have  now  l>een  detected  with  electron 
spin  resonance  (ESR)  in  many  irradiated  organic  and  bio¬ 
chemical  compounds  containing  sulfur.  Their  principal  g  values 
are  usually  in  the  range  of  2.003-2.065.  However,  a  sulfur- 
centered  radical  with  an  exceptional  g  tensor  having  principal 
values  ranging  from  1  99  to  2  29  was  observed  by  Akasaka  ( 1 ) 
in  single  crystals  of  l  .-cysteine- HC'.L  irradiated  with  1.5  MeV 
at  77  K  (one  MeV'  equals  1.6  X  10-1  ’  joule)  Theg  tensor  was 
found  to  have  approximately  axial  symmetry  about  the  larger 
value.  In  1965,  the  structure  of  the  crystal  w  as  unknow  n,  and 
it  w  as  therefore  impossible  to  relate  the  observed  g  tensor  to  the 
molecular  structure.  In  the  following  year.  Box  et  al.  (2)  ob¬ 
served  radicals  with  principal  g  values  of  g„  =  l .985,  g,  = 
2.004,  and  g„  =  2.251  in  single  crystals  of  cysteine-HCl  irra¬ 
diated  with  UV  at  77  k  The  similarity  in  g  tensors  supported 
the  conclusion  that  these  radicals  were  probably  of  the  same 
form  as  those  observed  by  Akasaka  In  the  meantime,  the 
structure  of  the  cysteine-f  K  )l  crystals  was  obtained  from  x-ray 
diffraction  by  Ayyar  and  Srinivasan  (3)  This  made  it  possible 
for  Box  cl  al.  (2)  to  deduce  that  the  largest  g  value  had  the  di¬ 
rection  of  the  CIS  bonds  of  the  undamaged  molecule.  Later,  this 
radical  was  also  observed  by  Kayushin  t’l  al.  (4).  Recently 
Saxebytl  and  Herskedal  (5)  have  observed  sulfur  radicals  w  ith 
principal  g  values  of  1.990,  2.006,  and  2.214  (denoted  as  “R|" 
radicals)  in  single  crystals  of  N-acetyl-l.-cysteine  irradiated  at 
77  K  with  4.0  MeV  electrons. 

To  distinguish  between  the  different  types  of  sulfur  radicals 
and  to  gain  more  quantitative  information  about  them,  we  in¬ 
itiated  a  program  for  measurement  of  the  i!S  hyperfine  struc¬ 
ture  in  irradiated  single  crystals  of  sulfur  compounds  (6-8). 
I’nlike  the  g  tensor,  the  VIS  coupling  gives  a  reliable  distinction 
la-tween  disulfide  and  monosulfide  radicals  and  provides 
quantitative  values  for  the  spin  densities  on  particular  S  atoms. 
Because  of  the  difficulty  of  obtaining  sufficient  quantities  of 
i  orti|Kiituds  with  enriched  V,S  to  grow  the  required  single 

•  i  i.  t.ilitm  ENlt.  electron  spin  resonance 

-  1  ■  m  reprint  requests  should  Ire  addressed 


crystals,  we  decided  to  concentrate  our  efforts  on  improvement 
of  our  spectrometer  sensitivity  so  that  reliable  measurement 
could  Ire  made  on  crystals  with  liS  in  its  natural  abundance  of 
0.76%.  To  our  gratification,  we  found  that  we  were  able  to 
detect  the  ;V1S  lines  for  the  hyperfine  components,  with  ample 
sensitivity  (signal  response  alxmt  25  times  greater  than  back¬ 
ground  noise  level)  to  make  reliable  measurements  of  the 
coupling  for  certain  irradiated  cry  stals  In  some  instances  where 
the  i;S  splitting  is  small  or  the  spin-lattice  relaxation  time  is 
unfavorable,  the  !iS  components  are  masked  by  w  ings  of  the 
much  stronger  absorption  lines  of  the  normal  isotopic  species 
Other  factors  of  importance  to  detectability  of  the  WS  lines  are 
the  degree  of  g  anisotropy  of  the  various  species  of  radicals  in 
the  irradiated  sample  and  the  number  of  magnetically  distin¬ 
guishable  orientations  for  each  species.  We  have  now  succeeded 
in  measuring  the  i!S  hyperfine  structure  of  the  radicals  with 
the  unusual  g  tensors  in  crystals  of  With  cysteine-HCl  (H-O)  and 
N-acetyl-l.-cysteine.  The  measurements  prove  that  both  are 
monosulfide  radicals  w  ith  similar  electronic  structures.  A  de¬ 
scription  of  the  results  for  N-acetyl-l.-cysteine  is  given 
Irelow . 

EXPERIMENTAL  DETAILS 

The  N-acetyl-l.-cysteine  cry  stals  were  grown  by  slow  evapo¬ 
ration  from  saturated  ILO  solutions  as  was  done  by  Saxebysl  and 
Herskedal  (5).  They  were  given  200  Mrad  dosages  of  50  keV 
x-ravs  at  77  K  and  were  observed  w  ith  an  X-batrd  (9.3  GHz) 
spectrometer  like  that  employed  in  our  earlier  work  (6).  One 
Mrad  equals  10‘ J/Kg. 

Samples  were  annealed  at  300  K  to  destroy  a  carbon-centered 
radical  (5).  Because  the  linewidth  was  too  great  to  permit  ob¬ 
servation  of  the  Ri  sulfur-centered  radical  at  room  temperature, 
the  ESR  sjx-ctra  were  obtained  at  77  K. 

Although  the  detailed  structure  of  acetyl  cy  steine  is  unknown, 
our  preliminary  x-ray  diffraction  shows  that  the  cry  stal  is  tri¬ 
clinic  with  one  molecule  in  unit  cell  An  orthogonal  uvw  coor¬ 
dinate  system  coinciding  with  the  principal  directions  of  the 
g  tensor  was  chosen  for  tlu-  hy  perfine  measurements.  Our  u  axis 
corresponds  to  the  Z  axis  employed  by  Saxebol  and  Herskedal 
and  our  u  and  u  axes  are  about  28  from  their  Y  and  ,Y  axes, 
respectively. 

THE  MONOSULFIDE  RADICAL 

This  is  the  sulfur  radical  (R|)  with  the  exceptional  g  tensor 
described  aliove  It  is  produced  by  the  irradiation  at  77  K  and 
is  relatively  stable  to  temperatures  well  above  77  K.  having  a 
room-temperature  half-life  of  the  order  of  I  day  However,  the 
radical  decay  s  rapidly  al  temperatures  of  400  k  and  above  Our 
measurements  confirm  the  unusual  g  tensor  observed  by  Sax- 
ebytl  and  Herskedal  We  did  not  improve  their  rather  accurate 
values:  g„  =  I  990,  g,  =  2.006,  and  g„  =  2  214 
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The  ESR  specirum  of  this  radical  has  an  essentially  isotropic 
proton  hyperfine  structure  consisting  of  four  lines  arising  from 
two  protons  with  couplings  of  22  G  and  29  G.  (One  gauss  equals 
1  X  1 0~*  tesla  )  This  structure  is  superimposed  on  each  i!S 
hy  perfine  component  and  provides  a  ready  identification  for 
them.  Fig.  1  shows  a  typical  pattern  observed.  In  the  central 
region  is  the  characteristic  |>attern  for  the  abundant  “-S  species, 
with  the  receiver  gain  reduced  by  a  factor  of  1/500.  Beyond 
the  breaks  in  the  curve,  on  either  sitle.  where  the  gain  is  500 
times  greater,  are  components  of  the  V’S  species  as  indicated. 

The  »S  hyperfine  splitting  was  found  to  be  maximum  along 
the  axis  of  g„.  the  minimum  principal  g  value,  and  to  be  very 
nearly  axially  symmetric  about  this  axis.  Fig.  2  shows  plots  of 
the  observed  displacements  of  the  V’S  components  as  a  function 
of  angular  deviation  from  the  axis  of  g„,  as  far  as  the  compo¬ 
nents  could  lx*  followed. 

There  is  a  small  (3  MHz)  low-field  shift  of  the  center  of  the 
’  ’S  patterns  from  the  ’-S  spectra  due  to  a  second-order  Breit- 
Rabi  type  of  interaction.  However,  this  interaction  did  not 
measurably  influence  the  component  spacings  which  were 
fitted  to  the  first-order  formula  expressed  in  the  principal  axes 
system  as 

A  =  +  g,-A,-l  -  +  g  a  ji  -t1  ■  111 

where 


< -  GAIN  *  500  — !|< -  GAIN  »  1  - 1(1 —  GAIN  »  500  — > 

Ku:  I.  Second-derivative  HSR  spectrum  for  the  monosulfide 
radical  Ki  observed  with  magnetic  field  in  the  tut  plane  '2(1°  from  the 
it  axis  (direction  of  minimum  tt  I-  The  four-line  proton  pattern  for  the 
unpaired  spin  localized  in  a  ;  S  atom  is  shown  in  the  center  of  the 
spectrum  at  reduced  gain.  The  -  '  and  +  /  ;  S  components  can  he 
seen  to  the  right  and  left  at  high  gain.  The  ± 1  components  are  com 
pletelv  obscured  by  the  stronger  S  pattern  and  radicals  stable  at 
room  temperature. 


g  =  +  g  -V  :  +  g. [2j 

Here  the  I  s  are  the  cosines  of  the  angles  the  external  field  makes 
with  the  orthogonal  uvw  coordinate  system.  A  least-squares 
fitting  of  Eq.  1  to  the  observed  hyperfine  splittings  w  itlill  in 
the  tit  and  titc  planes  yielded  (see  T able  1 )  principal  values  of 
199  ±  1,  38  ±  3.  and  32  ±  6  MHz  for  the  y,S  coupling.  Note  that 
the  coupling  is  of  approximately  axial  symmetry  about  the 
largest  element.  We  therefore  set 

A,  =  A ,  =  199  MHz.  [3] 

.1.  2-1,1  +  .4,1/2  =  35 MHz.  1 4 j 

Substitution  of  these  values  in  the  axially  symmetric  coupling 
formulas, 

.1,,  =  A  +  2 B  and  .4,  =  A .  -  H.  1 5 1 

yields 

A  =  90 MHz  and  K  =  54.7  MHz 

where  A/  is  the  isotropic  component  and  B  is  the  anisotropic 
component  of  the  ’  ’S  coupling.  These  results  are  summarized 
in  Table  1. 

We  can  safely  assume  that  the  anisotropic  component  B  arises 
primarily  from  3 p  orbital  spin  density  on  the  S.  Although  the 
g  tensor  gives  evidence  for  some  contribution  from  the  3 d  or¬ 
bital,  the  3 d  spin  density  is  probably  not  greater  than  0.20,  and 
the  coupling  by  a  3 d  electron  is  an  order  of  magnitude  less  than 
that  of  a  3 p  electron.  With  the  known  atomic  3 p  orbital  cou¬ 
pling  of  75  MHz,  the  3p  orbital  spin  density  is  thus  found  from 
the  anisotropic  coupling  to  be  p  j,,  =  54.7  78  =  0.70  Similarly, 
the  isotropic  component  with  the  know  n  3s  orbital  coupling 
indicates  that  p =  90  2700  =  0  03  on  the  S  This  3s  spin 
density  is  of  the  order  of  that  expected  from  a  a  ir  configuration 
interaction. 

Although  the  ’*S  coupling  proves  that  the  spin  density  is 
concentrated  mostly  on  a  single  sulfur,  the  g  tensor  of  the  mo- 
nosulfide  radical  observed  in  acety  I  cysteine  indicates  that  this 
radical  is  not  of  the  same  type  as  the  neutral  radical  RCHgS 
observed  at  sub-room-temperatures  in  irradiated  cystine  (6). 
Akasaka  (I)  has  proposed  that  the  approximately  axial  sym¬ 


metry  in  g  and  the  anomalously  large  g„  =  2  29  \\  hich  he  ob¬ 
served  for  irradiated  cysteine-HGl(HgO)  at  77  K  might  Ire  at¬ 
tributed  to  a  radical  of  the  type  R(  4  IjS  in  w  hich  there  is  a  near 
degeneracy  of  the  p,  and  p„  orbitals  However,  this  model  does 
not  account  for  the  fact  that  the  minimum  g  is  less  than  the  free 
spin  value  nor  for  the  large  deviations  of  the  other  principal  g 
values  from  those  of  the  RC4GS  radical  in  cystine-HCl  We 
propose  that  the  two  observed  monosulfide  radicals  are  basically 
distinct,  that  the  one  observed  in  acetyl  cysteine  utilizes  five 
valence  orbitals  of  the  S  w  hereas  the  one  found  in  irradiated 
cystine  utilizes  only  four.  The  five-orbital  S  radical  may  !«•  a 
negatively  charged  molecule  (RC.1GSX)-  or  a  neutral  radical 
RCfGSXj  in  which  the  X  is  li  in  the  radicals  discussed  here. 

Probably  the  radical  observed  in  l.-cysteine-HCl(HvO)  at  77 
K  is  the  negatively  charged  molecule  (RCIGSH)-  with  the  odd 
electron  concentrated  on  the  S  In  N-acetyl-t .-cysteine  the  sulfur 
radical  is  somew  hat  more  stable  than  it  is  in  cysteine  and  has 
a  lower  g,,  (2.214  as  compared  w  ith  2.29).  Hvdrogen  bonding 
involving  the  S.  as  indicated  in  structure  (1)  may  help  to  stabilize 
the  negati  .c  charge  on  the  S,  or  perhaps  an  internal  shift  of  an 
H  produces  a  more  stable  form  (II), 
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in  which  the  negative  charge  is  on  an  oxy  gen.  The  dotted  lines 
indicate  hydrogen  bonds.  A  possible  alternative  is  transfer  of 
a  proton  from  the  GOOIl  of  an  adjacent  molecule  in  the  cry  stal 
to  the  negatively  charged  SH  group  of  (I)  to  form  the  neutral 
radical  (III) 
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Kill.  2.  Plots  of  the  :  S  hyperfine  splitting  tor  the  monosulfide  radical'with  H  imposed  in  the  uc  (left  >  and  uu  (right  I  planes.  The  splittings, 
plotted  in  MHz  because  of  the  large g  variation  in  the  utc  plane,  are  measured  from  the  center  of  the  S  resonance.  The  circles  are  data  points: 
the  solid  lines  are  computer  fits  to  the  splitting  between  the  — %  and  +  i (_>  components.  The  t  !S  pattern  is  shifted  about  3  MHz  (downfieldl  from 
the  center  of  the  ;'S  pattern  by  second-order  effects. 


Without  specific  know  ledge  of  the  locations  of  the  atoms  in  the 
crystal  one  cannot  predict  which  form  is  the  more  probable. 

Although  a  radical  of  type  (1),  (II),  or  (111)  would  have  an¬ 
isotropic  SH  proton  hyperfine  interactions,  approximate  cal¬ 
culations  indicate  that  these  would  not  cause  resolvable  splitting 
and  that  the  resulting  spread  would  fall  within  theoliserved  line 
widths.  Furthermore.  Saxelxd  and  Herskedal  proved  with 
deuterated  crystals  that  the  observed  four-line  spectrum  was 
associated  with  the  OF  protons  rather  than  exchangeable 
hydrogens. 

The  important  characteristic  of  the  forms  (1),  (II),  or  (III)  is 
the  utilization  of  five  orbitals  of  the  S.  It  appears  that  an  electron 
captured  by  RSH  g<x*s  into  a  p-d  hybridized  orbital  of  the  S. 
whereas  one  captured  by  RSSR  goes  into  an  antibonding  a* 
orbital  formed  by  a  linear  combination  of  p  orbitals  of  the  two 
sulfurs  (7).  Evidently  a  similar  configuration  utilizing  an  anti¬ 
bonding  orbital  of  the  CS  or  SH  IhuhI  of  cysteine  has  higher 
energy  than  one  utilizing  a  fifth  atomic  orbital  on  the  S  formed 
by  a  hybridization  of  lid  and  3 p  orbitals.  Thus,  the  "three- 
electron  Ixind"  d<x*s  not  appear  to  lx*  involved  in  the  formation 
of  three  bonds  to  the  sulfur,  but  rather  the  3 d  shell  is  employed 
to  expand  the  numlier  of  Ixmding  orbitals  of  the  S  to  accom¬ 
modate  the  additional  Ixmd. 

The  many  unknown  parameters  prevent  determination  of 
the  specific  amounts  or  kinds  of  (/-orbital  hybridization  that 


Table  1 .  "S  couplings  and  S  spin  densities 
of  tile  monosulfide  radical  K1 


Radical  R  1 

*‘S  coupling  (MHz) 

S  spin  density* 

KOI. S  II 

A„  =  199  •  3 

or 

A,,  =  38  •  3 

A,  =  90 

Pa,  =  0.03 

RC’H.SII 

.4,,.  -  32  '  (i 

II  =  53 

Ps„  =  0.70 

occur  in  the  (RCIFSI I  )  or  RCIFSIF  radical,  but  normaliza¬ 
tion  of  the  unpaired  electron  density  with  the  3p  spin  density 
of  70%  and  the  10%  positive  spin  density  on  OIF  revealed  by 
the  proton  coupling  indicates  that  the  3 d  character  of  the  un¬ 
paired  electron  orbital  is  approximately  20%.  This  amount  is 
adequate  to  account  for  the  observed  reduction  of  g„  below  the 
free-spin  value  and  for  the  other  unexplained  features  of  the 
g  tensor.  The  3s  spin  density  of  3%  on  the  S  is  omitted  from  the 
normalization  liecause  it  is  likely  due  to  spin  polarization  of  the 
a  I xinds  and  is  hence  counterbalanced  by  an  equal  and  op|x»ite 
spin  density  of  these  Ixinds. 

The  principal  g  values  may  lx*  predicted  from  the  second- 
order  formula 


where  </>„  includes  all  orbitals  which  are  not  orthogonal  to  </></ 
=  1. 1 1  (/in),  i.e.,  for  which  (</>„  |  </>(/)  0.  Thus,  in  this  problem 

</>„  must  include  the  3 d  orbitals  on  the  S. 

The  ground  molecular  orbital  of  the  unpaired  electron  max 
lx*  expressed  as 

'/>,.  =  ps'  (^s"  +  Tai  ).  1 7 1 

where  ps  is  the  unpaired  electron  density  on  the  S  and  2)  a,  -  = 
1  —  ps  is  that  on  the  other  atoms.  Because  the  interactions  on 
the  other  atoms  i  do  not  contribute  significantly  to  the  aniso¬ 
tropy  in  g.  they  can  lx*  neglected  in  the  evaluation  of  the  last 
term  of  Eq.  (i  As  an  approximation  based  on  the  observed  spin 
density  and  S  model  deseriix'd  above,  we  express 

i/s"  =  (4/3)'  (3p(l  +  (1/3)'  <3 </,,).  [8 

( lontributions  of  the  3 p  and  3 d  components  of  the  hybrid  can 
be  calculated  separately.  Because  T ,  1 3p ,)  =  0  and  |3r/v„) 
=  i(3(/,;).  w  ith  I/,  neglected,  we  can  set 

/.,(</>„)  =  ps'  -'(/., j^s")  =  i/tj  (1  "»)'  (({</,  )  |9| 


*  A  '.id  spin  density  of  0.20  is  also  predicted  (see  text  I. 
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and  from  Eq.  6  can  obtain 


The  last  form  follows  from  the  normalization  (</,;  | </, ; )  =  I  and 
from  the  orthogonality  oi  all  other  combinations.  Because  the 
(/  shell  is  less  than  half-filled,  the  d  spin-orbit  coupling  A,/  is  here 
positive  in  sign  Thus,  the  contribution  of  the  d  character  re¬ 
duces  g„  lielow  the  free-spin  value  of  2.0023,  in  agreement  with 
observation  We  do  not  know  the  specific  value  of  the  spin-orbit 
coupling  constant  for  the  ,'W  orbital  ol  S  nor  /■.',/  —  TV  I  lowcvcr. 
with  Eq.  10  wc  can  evaluate  experimentally  the  ratio  A,/  (£,/ 
—  To)  I  rom 

g  -  g  =  -2A  ,,y  |  ,‘x  T,  -  T'„l  |  =  -0.0123  |lll 

Although  the  I  Minding  orbitals  may  have  some  s  hybridization 
and  [M'rhaps  also  some  d  character,  we  suggest  as  a  simplifying 
approximation  that  the  other  occupied  orbitals  ol  the  S  are: 

0,  =  (12)'  •'[(!  2)'  (.?/>,  >  +  (I  10)'  (3/i,) 


+  (2  5)'  (3d,  i 

'  +  4n  |.  .121 

0.  =  (1  2)1  1  — 1 1  2)'  (5/1.) 

-  (1  10)'  "(3/1,) 

+  (2/5)'  (3m/,  , 

'  +  •*/„  1*  1131 

0  =  ( 1  2)'  |(3/i  i  +  4,  1. 

|14| 

(/>,  =  ( ;>>•). 

1 15  | 

All  these  orbitals  are  assumed  to  have  electron  pairs  Orbitals 
0,  and  </>j  would  bond  the  two  1 1  atoms  of  structure  (II)  or  (III) 
while  </>i  would  form  the  a  bond  to  the  ( )  and  0 1  would  have  an 
unshared  pair 

Eq.  6  with  the  orbitals  ol  Eqs.  7,  and  12-15  gives: 


w  hen  it  is  assumed  that,  on  the  average,  one  electron  ol  each 
bonding  pair  is  on  the  S  With  the  term  in  A,/  given  by  Eq.  1 1, 
with  /is  =  />  i/i  +  />.w  =  d  70  +  0.20.  g /  =  2.0023.  T  i  =  Ts,  anil 
A,,  =  — 3S2  cm  ’,  these  expressions  reduce  to: 

g  =  unto  +  550  (/■:  -  /■;  >.  l  is  | 

g.  =  l.tttH)  +  550  t  T  -  T',,1.  1 19 1 

With  the  measured  values  g,  -  2  (KMi.  g„  =  2  211.  these 
equations  give  TV  —  TV  =  31.000  cm  1  and  P.\  —  To  =  2100 
cm  1  =  /V-  TV  In  addition  to  the  other  approximations,  it  is 
assumed  that  the  spin  density  of  0. 10 on  the  (  II  •  group  con¬ 
tributes  nothing  to  the  g  anisotropy  Despite  the  small  energy 
difference  between  T|  or  TV.  anil  TV.  tin*  Boltzmann  factor 
exp(  — AT  kT)  is  000031  at  300  k.  and  the  calculated  spin 
density  in  orbital  0/  or  <t>2  is  only  0  0003.  not  enough  to  give  an 
observable  SI  I  proton  splitting,  nor  a  measurable  deviation  in 
the  axial  svmmetry  of  the  V,S  coupling. 

In  structure  (I),  the  negative  charge  is  on  the  S.  and  there  are 
two  electrons  in  the  S  orbital  forming  the  S  II  O  bridge 
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Eli;.  :t.  KSH  spectrum  of  stable  disulfide  radicals  oliserx'ed  at 
room  temperature,  lu  tin*  left  -hand  part  of  the  spectrum  are  shown 
the  S  patterns  observed  with  low  amplification.  The  crystal  is  ori 
cnletl  with  the  magnetic  field  in  the  direction  of  the  minimum  g  value 
IhrlheSt  radical.  A  portion  oft  he  "S  spectrum  for  Si  is  shown  at  high 
gain  in  the  right  hand  portion  of  the  spectrum.  The  spectrum  reveals 
tiuetpial  coupling  to  two  sulfur  atoms.  S' 1 '  and  S'  '.  wit  bin  t  he  radi 


If  the  same  orbitals  arc  assumed  with  0i  forming  the  SI  I  bond 
and  0g  having  an  unshared  pair,  and  il  the  approximation  T  i 
*  fV  is  made,  then  a  similar  calculation  gives  g,  the  same  as 
before  but  g„  =  1.990  +  820  (T  i  —  T0).  The  observed  g„  then 
indicates  that  K\  —  TV  =  3(i(K)  cm  1 

Toeliminate  unknow  n  parameters  in  these  solutions,  xve  w  ere 
forced  to  make  rather  severe  simplifying  assumptions  W  hat¬ 
ever  the  nature  ol  hybridization  of  the  orbitals.  TV  —  /Valid  T.j 
—  To  must  lie  small  and  T  ;  —  To  must  be  large  to  account  for 
the  observed  g,  andg„  values  Similar  conditions  are  rciptired 
for  the  ttiixlel  proposed  by  Akasaka  ( I ).  xvlio  assmncil  T,„  s  T,„, 
«  T,, 

In  ih c  earlier  work  (f>),  we  found  that  the  isotropic  (  lb 
proton  coupling  A  in  the  neutral  BOIES  and  BObSS  radicals 
could  lit*  calculated  with  the  relation 

.1  =  .1,  /y  cos  o  20 

with  .Ao  =  31.5(l.  If  one  assumes  the  same  1 1C  'll  angle.  125  . 
for  the  present  radical  as  for  IU  IDS  and  a  /is  =  0.90.  the  com¬ 
bined  3/1  and  Ad  density .  Eq  20.  requires  a  somew  hat  higher 
value. . A  o  =  3(i.(i(),  to  account  for  the  oliscrx  ed  proton  couplings 
of  22.1  (I  and  29.4  C  Perhaps  this  signifies  a  more  effective 
hvpercoujugation  by  tin*  /i -d  hybrid  orbital,  but  this  conclusion 


■ii*''  ?o°  cf  2Cf  atf 

ANliLE  OF  Field  from  gu  DIRECTION 


Kul  I  I’lotuftbi*  S  byperfim*  splitting  versus  angle  of  the  ex 
lornal  field  from  tin*  direct  inn  of  the  minimum  principal  g  value  for 
the  disulfide  radical  SI.  Circles  represent  measured  splittings  from 
tile  S  resonance;  the  curves  are  computer  fits  to  the  -  '  components 
lor  each  of  tile  two  sulfur  sites.  The  splitting  is  axially  symmetric 
about  the  g  ,  direction. 
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Table  2.  "S  couplings  and  S  spin  densities 
of  the  disulfide  radical  St 


,JS  coupling  ( MHz) 

S  spin  density 

Radical  SI 

•V 

B 

3s 

3 p 

RCH.StDSO) 

SO) 

29.7 

23.0 

0.011 

0.29 

S<2> 

59.  6 

14.8 

0.022 

0.57 

is  uncertain.  Whatever  the  cause,  the  CH*  couplings  are  sig¬ 
nificantly  larger  for  the  current  radical  than  for  the  K(  .1  US  and 
RGIUSS  radicals 

THE  DISULFIDE  RADICALS 

Fig.  3  shows  the  ESR  spectrum  for  one  orientation  of  an  acetyl 
cysteine  crystal  which  has  been  aged  for  1  month  after  irra¬ 
diation.  In  the  low-gain  region  there  are  lour  distinct  signals 
of  sulfur-centered  radicals  of  the  normal  isotopic  species:  a 
triplet  indicated  by  T1  and  three  singlets  indicated  by  SI,  S2. 
and  S3.  The  principal  g  values  for  T1  are  2.003,  2.024,  2.053; 
for  SI,  2.002.  2.025,  2.064;  for  S2,  2.002.  2.025.  2.054;  for  S3, 
2.003,  2.026,  2.054.  They  fall  into  the  tw  o  sets  of  g  values  earlier 
observed  for  disulfide  radicals  (6-9).  differing  significantly  only 
in  the  maximum  value  which  is  either  close  to  2.054  or  to  2.064. 
For  this  reason,  all  four  signals  are  presumed  to  arise  from  di¬ 
sulfide  radicals  although  this  structure  could  lx’  verified  by  viS 
hyperfine  structure  only  lor  the  SI  signal. 


The  crystal  orientation  for  Fig  3  is  such  that  the  applied 
magnetic  field  is  along  the  direction  of  g„.  the  minimum  g 
value  for  the  SI  signal  The  signals  to  the  right,  observed  under 
high  gain,  are  1!S  hyperfine  com|x>nents  corresponding  to  the 
SI  s[X'cies.  We  were  able  to  observe  these  liS  components  over 
a  sufficient  range  of  crystal  orientations  (see  Fig  4)  to  obtain 
the  principal  values  of  the  coupling  as  given  in  Table  2.  The  two 
different  coupling  values  for  ,:iS  prove  the  disulfide  structure 
of  the  radicals  as  explained  in  the  first  paper  of  this  series  (6). 
The  spin  densities  on  the  two  sulfurs  obtained  from  the  coupling 
values  are  also  listed  in  Table  2. 
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ABSTRACT  Hyperfine  structure  due  to  170  has  been  mea¬ 
sured  in  irradiated  zein  and  edestin  after  exposure  to  gaseous 
oxygen  with  l70  concentrated  to  24%.  The  observations  prove 
that  the  free  radicals  produced  by  ionizing  irradiation  under 
vacuum  at  300  K  are  converted  by  the  oxygen  to  protein-per¬ 
oxide  free  radicals  X — O*1* — 02i,  with  the  unpaired  electron 
density  in  a  rr-type  orbital  predominantly  on  the  peroxide  group. 
From  the  observed  couplings,  the  2 p,  spin  densities  on  zein- 
peroxide  are  found  to  be  0.29  and  0.45  for  O"1  and  O*21,  respec¬ 
tively;  those  on  edestin-peroxide  are  0.26  and  0.4S  for  0*11  and 
O*21,  respectively. 

In  the  first  observations  of  electron  spin  resonance  (ESR)  signals 
in  irradiated  proteins  (1),  a  proton  hyperfine  doublet  produced 
from  bone  collagen  was  found  to  be  converted  to  a  singlet  when 
the  specimen  irradiated  in  vacuum  was  exposed  to  oxygen.  It 
was  postulated  that  that  O2  combined  with  the  free  radicals  in 
the  irradiated  sample  to  produce  a  protein-peroxide  free  rad¬ 
ical.  Since  this  initial  finding,  many  examples  of  this  oxygen 
effect  have  been  observed  for  other  proteins  (2)  as  well  as  for 
biochemicals  of  other  types  (3).  In  certain  proteins  such  as  zein 
and  edestin,  chosen  for  the  present  l70  study,  the  conversion 
of  the  proton  multiplet  to  the  oxygen  singlet  is  rapid,  with 
complete  conversion  to  the  singlet  occurring  within  minutes 
after  exposure  of  the  sample  to  atmospheric  oxygen.  In  others, 
most  notably  silk,  the  oxygen  effect  is  slow,  requiring  days  of 
exposure  for  partial  conversion.  This  slower  conversion  is  be¬ 
lieved  to  result  from  steric  factors  that  hinder  or  prevent  the 
molecules  of  oxygen  from  reaching  the  free-radical  centers  in 
fibrous  proteins  such  as  silk.  The  oxygen  effect  is  not  observed 
for  sulfur-centered  radicals  produced  by  irradiation  of  proteins, 
probably  because  molecular  oxygen  does  not  combine  with  the 
sulfide  free  radicals. 

The  earlier  studies  of  the  oxygen  effect  in  irradiated  proteins 
were  made  with  1bO,  which  has  no  hyperfine  structure.  Al¬ 
though  there  seems  to  be  little  doubt  that  the  radicals  formed 
upon  admission  of  oxygen  are  peroxide  radicals,  the  present 
study  of  1  'O-labeled  radicals  unquestionably  confirms  this  in¬ 
terpretation  and,  in  addition,  provides  new  information  alxmt 
the  electronic  structure  of  the  protein-peroxide  radicals. 

In  many  irradiated  proteins,  including  edestin  and  zein,  the 
free- radical  site  giving  the  F.SR  pattern  in  the  absence  of  oxygen 
is  believed  to  be  the  C„— R  group  in  structure  I 


H  R  o 

_  M  _l_ 


H  K  o 


formed  by  loss  of  a  If  from  a  carbon  of  the  polypeptide  chain. 
For  edestin  under  vacuum,  a  proton  doublet  like  that  for  the 


C„I1  fragment  in  irradiated  polyglycine  is  formed  (2),  which 
indicates  that  R  =  If  and  that  the  radical  is  formed  in  the  gly¬ 
cine  unit.  A  similar  doublet  is  observed  for  evacuated  zein,  but 
w  ith  a  superimposed  quartet  like  that  for  polyalanine  (2).  This 
indicates  formation  of  radical  sites  in  lx>th  glycine  and  alanine 
(R  =  CH.j)  units  of  zein.  Relative  intensities  of  the  superimposed 
patterns  indicate  that  formation  in  the  alanine  unit  is  somewhat 
less  than  that  in  the  glycine  unit.  Both  patterns  are  converted 
to  singlets  when  normal  oxygen  ( "’Os)  is  admitted  to  the  irra¬ 
diated  sample. 

Upon  exposure  of  the  irradiated  samples  to  oxygen,  peroxide 
radical  II  is  evidently  formed.  When  this  occurs,  as  the  present 
study  confirms,  the  spin  density  of  I  is  effectively  transferred 
to  the  captured  O2,  where  its  interaction  with  the  protons  of  the 
R  group  is  insufficient  to  produce  a  resolvable  proton  splitting. 
Although  lfiO  has  no  hyperfine  structure,  the  substituted  itO 
has  nuclear  spin  /  =  5/2,  which  gives  rise  to  a  sextet  ESR 
splitting  bv  each  of  the  inequivalent  oxvgens  in  peroxide  radical 
II 

Previously,  nO  hyperfine  structure  has  been  observed  for 
randomly  oriented  peroxide  radicals  of  irradiated  polytetra- 
fluoroethylene  (4)  and  for  oriented  peroxide  radicals  C10H 1  iOO 
in  single  crystals  (5).  Also,  the  isotropic  component  of  the  170 
hyperfine  structure  has  been  observed  for  some  simpler  per¬ 
oxide  free  radicals  in  liquid  solutions  (6,  7).  These  earlier  ob¬ 
servations  have  been  helpful  to  us  in  our  analysis  of  the  l70 
hyperfine  structure  for  the  protein-peroxide  radicals. 

EXPERIMENTAL  PROCEDURES 

A  balanced-bridge,  X-band  ESR  S|>ectrometer  operated  at  a 
frequency  of  9000  MHz  was  used  for  these  measurements. 
Evacuated  tubes  made  of  Spectrosil  quartz  held  the  protein 
powders,  which  were  placed  in  a  B0Co  7-ray  source  for  irra¬ 
diation  at  300  K.  Doses  were  approximately  3  X  105  rads  (300 
grass).  Afterward,  the  tubes  were  annealed  in  a  gas/oxvgen 
flame  to  eliminate  any  free  radicals  formed  in  the  quartz. 

Exposure  of  the  irradiated  samples  to  normal  gaseous  oxygen 
or  to  1 '()  concentrated  gas  was  accomplished  in  this  way.  The 
tube  was  connected  to  a  flask  of  the  oxygen  gas,  the  glassware 
between  them  was  evacuated,  and  the  flask  was  immersed  in 
liquid  nitrogen.  When  the  valve  on  the  flask  was  opened,  the 
oxygen  at  its  va|x>r  pressure  at  77  K  (about  130  torr)  flowed  into 
the  sample  tube.  After  the  sample  was  disconnected  from  the 
flask,  it  was  ready  for  observation  in  the  spectrometer  Con¬ 
version  of  the  radicals  formed  under  vacuum  to  peroxides  re¬ 
quired  2-4  hr.  The  enriched  oxygen  gas  contained  24%  ,70 
atoms. 


Abbreviation:  ESR,  electron  spin  resonance. 

*  Present  address:  Schixil  of  Engineering,  Duke  University,  Durham, 
NC  27706. 
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OBSERVED  SPECTRA 

Fig.  1  shows  a  comparison  ot  the  second-derivative  spectra  of 
irradiated  zein  after  exposure  to  normal  oxygen  gas  and  after 
similar  exposure  to  gas  vv  ith  24rr  substitution  of  1  ()  atoms  for 
lhO  atoms.  The  spectrometer  amplification  for  the  two  curves 
was  approximately  the  same  The  stick  diagrams  at  the  base 
correspond  to  the  theoretically  exjxvted  1 '()  hy  [X'rfine  patterns 
of  the  two  forms.  \  "X)  l7Oand  X— 17() — "X)  Hyperfine 

eom|X)nents  for  the  less-probable  species.  X — 1  •()  1  (),  pros od 

to  lx-  t<x>  weak  for  detection  The  pattern  vv  ith  the  \x  ider  span 
is  that  lor  X  — "’()  ’’O.  Fig.  2  shows  the  curve  obtained  for 

irradiated  edestin  after  exposure  to  the  '"O-substitnted  gas 
The  nature  ol  the  F.SR  signals  and  ol  the  I  roe  radicals  formed 
upon  irradiation  in  the  absence  of  oxy  gen  are  described  in  the 
Introduction.  For  both  proteins,  the  primary  radicals  are  be¬ 
lieved  to  have  the  general  structure  1  and  the  (XToxide  radicals 
(secondary)  to  have  structure  II.  The  hyperfine  structure  and 
g  values  for  these  protein  peroxide  radicals  are  explained  in 
the  follow  ing  paragraphs. 

HYPERFINE  STRUCTURE 

The  analysis  of  FSH  h\  (X'rfine  structure  for  randomly  oriented 
tree  radicals  is  now  well  know  n  (ref.  S.  pp.  353-44 1 V  The  sha|X‘s 
of  the  observed  absorption  curves  depend  on  several  factors, 
the  most  important  of  which  are  the  anisotropy  in  the  g  tensor 
and  in  the  nuclear  coupling  tensors  In  the  protein  peroxide 


Km.  I.  Second  derivative  KSH  spectrum  of  y  irradiated  zein 
after  exposure  to  gaseous  hiO-  1 1  pprr)  and  alter  exposure  to  oxygen 
gas  enriched  with  24r<  1  ()  I /.nicer). 
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Kir.  2.  Second  derivative  K KK  spectrum  of  y  irradiated  edestin 
after  exposure  to  oxygen  gas  enriched  with  24rl  '  4). 


radicals  considered  here,  the  isotropic  com|x>ncnt  ;\f  of  the  1  () 
coupling  is  ot  comparable  magnitude  to  the  anisotropic  cou¬ 
pling  /I  and  isol  like  sign  Conseiiuently.  the  resultant  coupling 
(A/  —  II)  for  the  magnetic  field  ol  orientations  perpendicular 
to  the  p „  orbital  of  the  unpaired  electron  is  very  small.  As  a 
result,  the  critical  absorptions  corresponding  to  the  perpen¬ 
dicular-field  orientation  are  not  resolvable  and  are  not  separable 
from  the  strong  "X)  single!  absorption.  To  obtain  a  value  for  II 
one  must  observe  the  much  w  eaker  critical  absorptions  corre¬ 
sponding  to  orientation  ol  the  field  approximately  parallel  to 
the  coupling  />„  orbital  of  the  1 '()  atoms. 

The  nature  of  the  1  ()  by  pet  fine  spectra  is  illustrated  by  the 
ideali/eil  curves  ol  Fig  3.  in  w  Inch  line-broadening  factors  are 
neglected  except  those  ol  the  anisotropy  itt  the  nuclear  coupling 
ami  in  g.  Also,  lor  simplification  of  llirdiufirMii  gas  well  as  A 
is  assumed  to  be  axially  sy  mmetrie  about  the  directions  of  the 
coupling/),  orbital.  The  weak  (X'aks  in  the  observed  derivative 
curves  (Figs  1  and  2)  correspond  to  the  outer  terminations  of 
the  ditferent  1  ()  component  absorptions  If  the  coupling  con¬ 
stant  is  expressed  in  magnetic-field  units,  these  outer  peaks 
occur  at  the  critical  field  values. 

II  =  lif'u  (g  ti)  A  A  Mi.  |1| 

w  here  i’u  is  the  o|X’rating  frequency  and  g  =  g„  is  the  effective 
g  \ alue  for  II  along  the  coupling  /)„  orbital,  and 

A  =  A,  +  2/1.  [2 1 

w  here  A  is  the  effective  coupling  constant  for  II  along  the 
coupling  orbital  The  much  stronger  unresolvable  inside- 
component  terminations  occur  at 

II  ,  =  hi'o  (£  i  ft)  +  A  j  Mp  |3| 

w  here  £  ,  is  the  effective  £  \  alue  for  II  imposed  in  the  plane 
perpendicular  to  the  p,  orbital,  and 

3  ,  =  A,  -  II  14 1 

\llhough  £  has  a  fixed  value.  £  j  ranges  in  value  from  g,.  to 
g„  However,  lor  simplicity  in  Fig  3.  g,  was  chosen  as  the 
particular  v  alue  lor  g  The  quantum  numbers  arc  Af/  =  ±5  2. 
A3  2.  A I  2.  corresponding  to  the  spin  I  —  5  2  for  1 '()  The 
coupling  constants  have  the  same  sign  and  are  both  negative 
because  the  nuclear  magnetic  moment  ol  *’()  is  negative 

The  fact  that  the  II  i  components  cannot  be  resolved  or 
separated  from  the  unsubstituted  "’()  absorption  indicates  that 
\f  -  II  With  A/  =  II.  the  observed  component  separation  of 
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Fie.  3.  Idealized  theoretical  pattern  of  F.SK  hyperfine  structure  for  a  coupling  l7()  of  randomly  oriented  peroxide  radicals  having  axial 
symmetry  in  coupling  with  .4  i  =  A/  -  H  =  2.5  ('.  and  with  .4  =  —VI  (!.  as  for  l7Ol21  of  edestin  peroxide.  For  simplicity,  the  g  tensor  is  also  assumed 
to  have  axial  symmetry  about  the  coupling  2p  orbital  with  g  =  g„  and  g  L  =  g,.  for  edestin.  lane-broadening  factors  other  than  anisotropies 
in  g  and  A  are  neglected. 


71  G  for  ''O'21  in  edestin-peroxide  gives  with  Eq.  2  the  value 
Af  =  II  =  23.7  G  for  the  terminal  oxygen;  and  the  observed 
component  separation  of  42  G  for  l7()'n  gives  Af  =  I)  =  14  G 
for  the  central  oxygen  However,  we  believe  that  more  accurate 
values  for  R  are  obtained  with  the  assumption  that  A/(l7()(l’) 
=  1(5  G  ami  Ay(’  ‘O'2’)  =  22  G.  as  measured  for  the  radicals 
(Gf f;,)2(XX)  and  GaiisfCI I  ;)A XX)  in  liquid  solution  (7  V  These 
values  are  close  to  the  isotropic  coupling  observ  ed  in  diverse 
peroxide  radicals  (5-7)  in  which  the  (V>  is  Itonded  to  a  carlxm 
With  these  values  assumed  lor  Af  in  With  zoin  and  edestin.  the 
anisotropic  coupling  constants  calculated  vv  ith  Eq  2  are  as  listed 
in  Table  1.  The  couplings  A  in  MHz  arc  related  to  the  com¬ 
ponent  separations  in  gauss  by 

A  (MHz)  =  1.40#  A//(G).  |5| 

PRINCIPAL  g  VALUES 

The  principal  g  values  of  these  protein  peroxide  radicals  cannot 
lx'  accurately  derived  from  the  randomly  oriented  radicals,  but 
approximated  values  have  been  deducetl  by  the  following 
procedure.  The  weak  [X’ak  to  the  left  of  the  strong  central  ab¬ 
sorption  in  the  unsubstituted  "’()  peroxides,  show  n  for  zoin  in 
Fig.  1  ( ’ pper .  evidently  corresponds  to  the  maximum  g  value, 
g,,..  The  value  determined  by  this  peak  is  gu  —  2044  for  the 
zein  peroxide.  The  value  similarly  obtained  for  the  edestin  is 
2.040.  These  values  are  in  the  range  of  those  derived  for  the 
peroxide  radical  in  single-crystal  measurements  I  ref.  5.  p.  002). 
For  example,  in  the  CiollnOO  radical.  g„  =  2.045  and  2.039 


Table  1 .  1 X)  coupling  constants  for  protein  peroxides. 
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Compo- 

Kffeet  ive 

Con- 

nent 

couplings.* 

Spin 

Parent 

pling 

separations 

MHz. 

densities^ 

protein 

atom 

AH  ((’.) 

•V 

a 

P2» 

pi/. 

Zein 

■7<)"| 

46  *  2 

t-)4.r> 

(-142 

0.010 

0.29 

■"O'21 

08  ±  2 

( -  )62 

(  -  )0f> 

0.012 

0.4.r> 

Kdestin 

IT0m 

42  *  2 

( -  )4f> 

(-137 

0.010 

0.20 

71*2 

(— )62 

( -  )69 

0.012 

0.48 

*  Signs  art*  not  measured  but  art*  theoretically  negative. 
f  Calculated  with  atomic  coupling  values  =  — 1028  MHz  and  It: 

=  -144  MHz  (ref.  S.  pp.  :t:l7.  338). 

1  Values  assumed  from  measurements  of  liquids  (see  text  t. 


for  the  two  crystal  sites  (5).  The  gu  values  are  expected  from 
theory  to  lx-  directed  approximately  along  the  O — O  bond  axis. 
Peaks  corresponding  to  the  intermediate  and  minimum  g  val¬ 
ues.  g,.  and  g„.  respectively,  cannot  lx-  separated  for  the  un¬ 
substituted  lhO  radicals;  both  of  them  are  ev  idently  within  the 
strong  central  absorption  with  the  center  corres|x>nding  to  2.007 
for  zein  and  to  2  (KKi  for  edestin.  However,  the  minimum  value 
g„  is  expected  from  theory  to  have  the  direction  of  the  p„  or¬ 
bital.  that  lor  vv  hich  the  1  ()  hy  perfine  components  are  mea¬ 
sured.  Thus,  it  slight  second-order  shifts  are  neglected,  the  gu 
value  corresponds  to  the  center  of  the  A  Mi  hyperfine  pattern 
like  that  shown  in  Fig  3  We  have  assumed  that  the  interme¬ 
diate  values  g,  correspond  to  those  for  the  strong  central  ab¬ 
sorption  for  the  unsubstituted  "’()  radicals  The  unresolved 
A  !  M  |  pattern  is  depicted  as  a  closely  spaced  multiple!  centered 
near  g,  in  Fig.  3.  This  multiplet  ev  idently  falls  under  the  strong 
central  "A)*  absorption.  For  orientations  of  the  magnetic  field 
along  or  near  g„  .  the  1  •()  components  are  too  weak  for  obser¬ 
vation  The  approximate  g  values  derived  by  these  procedures 
and  assumptions  are  listed  in  Table  2.  They  are  close  to  the 
values  2.(X)2b.  2.<KXi5,  and  2  03 S  derived  by  Ghe  and  Tench  (4) 
for  the  chain  radicals  — GIG  GF(OO) — GFg —  at  77  K 

CONCLUSIONS 

The  combined  spin  densities  on  O'"  anil  O'-’  deriv  ed  for  edestin 
and  for  zein.  0.74.  are  not  unity .  as  would  lxl  expected  for  an 
unpaired  electron  entirely  in  a  7r  orbital  on  the  two  oxygens  for 
vv  hich  orbital-overlap  distortions  are  neglected  This  discrep¬ 
ancy  might  be  interpreted  to  indicate  some  spreading  of  the  tt 
orbital  of  the  unpaired  electron  to  the  carlxm  or  to  another  atom 
ol  the  protein  Also,  some  reduction  in  the  effective  R  values 
and  derived  />,  spin  densities  is  expected  to  result  from  rather 
large  torsional  oscillations  of  the  |X'ro\ide  group  alxiut  theC  — O 
bond  Such  oscillations  would  reduce  the  observed  maximum 
coupling  below  that  expected  for  a  static,  parallel  alignment 
of  the  ;>  orbital  as  w  as  assumed  in  the  derivations.  If  the  entire 
reduction  is  attributed  to  these  oscillations  and  the  2p,  spin 

Table  2.  Approximate  g  values  of  protein  peroxide  radicals 


Parent 

protein 

Hu 

A»Y 

Hu 

Zein 

2.001 

2.007 

2.044 

Kdestin 

2.002 

2.000 

2.040 

! 
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density  on  the  two  oxygens  is  norinalizeil  to  unity,  the  iiulicateil 
2f)„  spin  densities  are  0.39  and  0.61  for  ()<n  and  0<-t  of  zein 


The  small  2s  spin  densities  are  attributed  entirely  to  spin-[xi- 
larization  effects  and  do  not  influence  the  calculations  of 
Pip- 

These  measurements  of  170  hyperfine  structure  verify  that 
the  ESR  signals  produced  by  exposure  of  these  irradiated  pro¬ 
teins  to  molecular  oxygen  are  due  to  peroxide  radicals,  as  earlier 
postulated. 
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Army  Research  Office. 


f 'roc.  Nall.  Acad.  Sci.  USA  77  (1980) 

1.  Gorily,  W.,  Aril,  W.  B.  &  Shields.  It  (1955)  Proc.  Natl.  Acad  Sci. 
I  S  A  41.  983  996. 

2.  Patten,  R  A  &  Cfordy.  W.  (1964)  Radial.  Res.  22, 29-44. 

3.  Rexroad.  H.  N.  &  Curdy.  W.  (1959)  I'roc.  Nall.  Acad.  Sci.  t  SA 
45. 256-209. 

4  Ghe,  M.  &  Tench,  A.  J  (1976)  J  Chew.  Phys  64,2370-2374. 

5.  Melainutl.  F...  Schlick,  S.  &  Silver.  B.  t..  (1974)  J.  Magn.  Rcson. 
14.  104  III 

6.  Fessenden,  It  W.  A  Schuler,  It  II  (1966)  J  ('hew.  Phys.  44. 
434-437. 

7.  Adamic,  K  .  lngolil.  K.  l!.  A  Morton.  J.  R.  (1970)  J.  Aw.  (.hew. 
S oc.  92. 922-923. 

8.  Gorily.  \V  (1979)  Theory  and  Applications  of  Electron  Spin 
Resonance  (Wilev-lnterscience.  New  York). 


